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ABSTRACT 


During  the  last  program  period,  September  I960  to  1  July  1961, 
theoretical  and  experimental  investigations  of  high  speed  gas  dynamics 
were  carried  out  under  this  program. 

The  experimental  investigations  included  stagnation  point  heat  trans¬ 
fer  measurements  in  low  Reynolds  number  flows  in  the  existing  helium- 
driven  shock  tunnel.  Electron  beam  densitometer  studies  were  continued. 
Modification  of  the  existing  hypersonic  shock  tunnel  to  incorporate  a 
combustion  driver  and  a  six-foot  diameter  nozzle  was  initiated.  This 
modification  will  allow  an  increased  performance  up  to  a  Mach  number 
equal  to  30  and  density  altitudes  to  300,  000  feet.  This  modification  will 
be  completed  by  January  1962. 

The  theoretical  investigations  included  a  study  of  viscous  and  invis- 
cid  non- equilibrium  hypersonic  flows  about  blunt  and  sharp  ieading  edge 
slender  bodies.  A  numerical  program  to  determine  the  non- equilibrium 
flow  about  a  blunt  body  in  the  vicinity  of  the  stagnation  region  was  also 
initiated.  A  study  of  the  hypersonic  wake  about  blunt  re-entry  bodies 
was  completed. 
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I.  INTRODUCTION 


In  October  of  I960,  the  Aerodynamics  and  Propulsion  Laboratory 
assumed  responsibility  for  conducting  research  operations  in  the  hyper¬ 
sonic  shock  tunnel  and  arc-driven  shock  tube.  The  experimental  program 
being  conducted  at  that  time  under  STL  contract  no.  AF  04(647)- 309  was 
continued  in  an  uninterrupted  fashion  by  Aerospace  personnel,  since  this 
effort  was  applicable  to  the  scope  and  objectives  of  the  High  Mach  Number 
and  Materials  Research  Program,  Phase  I,  planned  for  performance  under 
Aerospace  contract  no.  AF  Q4(647)-930.  * 

Under  Aerospace  management,  the  High  Mach  Number  and  Aerophysics 
Research  Program  was  increased  in  scope  to  include  a  more  extensive 
theoretical  effort  and  advanced  planning  "'"or  future  experimental  investiga¬ 
tions  and  facilities  requirements  to  meet  Air  Force  research  objectives  for 
advanced  missile  and  space  systems. 

Theoretical  work  conducted  under  the  High  Mach  Program  during  this 
reporting  periuu  imluues  extension  ot  the  hypersonic  wake  flow  studies  to 
fulfill  the  requirements  of  the  Penetration  Aids  Program,  viscous  and 
inviscid  non- equilibrium  gas  flows,  and  investigations  of  the  flow  field 
around  a  blunt  body  considering  finite  chemical  reaction  rates.  The  viscous 
and  inviscid  non- equilibr ium  flow  investigations  for  blunt  and  sharp  leading 
edges,  and  slender  body  power  law  geometries  were  extensive  in  nature, 
and  conclusive  results  were  reached  which  are  detailed  in  this  report.  The 
investigations  of  the  non- equilibrium  flow  field  about  a  blunt  body  include  a 
five- component  reacting  gas,  using  the  best  known  values  of  chemical 
reaction  rates.  Ionization  non- equilibrium  effects  are  not  presently  in¬ 
cluded,  but  provision  has  been  made  for  their  integration  into  the  program 
at  a  later  date.  The  program  i:-.  presently  in  the  check-out  process  on  the 
IBM  7090  computer.  A  major  application  of  this  study  will  be  in  the  area 
of  re-entry  plasma  interactions  with  microwave  radiation  in  support  of 
re-entry  communications  problems. 


The  experimental  studies  of  arc  heating  as  a  means  of  Increasing 
shock  tube  performance  were  completed  and  a  final  report  submitted.  The 
limitations  of  this  technique  for  producing  high  enthalpy  flows  were  clearly 
defined  in  the  study  and  detailed  in  the  report.  Future  effort  was  discon¬ 
tinued,  since  this  technique  was  not  optimum  for  providing  the  flow  simu¬ 
lation  capability  required  to  meet  the  objectives  of  the  High  Mach  Program. 
The  one-inch  shock  tube  used  in  these  studies  is  presently  being  used  for 
instrumentation  calibration. 

Performance  analysis  studies  were  performed  on  various  techniques 
for  generating  high- enthalpy,  high-velocity  flows.  These  studies  include 
arc  heated  driver  methods,  spark  discharge  heated  driver  methods,  and 
combustion  driver  techniques.  As  a  result  of  these  studies,  it  was  con¬ 
clusively  shown  that  the  combustion  driver  technique  was  the  most  practical 
method  of  generating  the  required  flow  in  a  research-type  facility.  As  a 
consequence,  a  modified  design  of  the  present  hypersonic  shock  tunnel  was 
completed  and  plans  initiated  to  convert  the  present  facility  to  incorporate 
a  combustion  driver.  The  expected  performance  of  the  modified  facility  is 
enthalpy  simulation  up  to  approximately  20,000  ft/sec.  The  expected  date 
of  completion  of  the  modification  is  December  1961, 

During  the  transition  period  to  combustion-driven  operation,  the  ex¬ 
perimental  investigations  on  stagnation  point  heat  transfer  in  hypersonic 
low  Reynolds  number  flows  have  been  completed  and  investigations  initiated 
on  control  surface  effectiveness  at  high-altitude  hypersonic  speeds. 

Future  plans  include  the  incorporation  of  a  six-foot  diameter  test  sec¬ 
tion  into  the  combustion-driven  facility  in  order  to  permit  larger  model 
size  and  extend  present  instrumentation  capability  for  low  pressure  meas¬ 
urement.  Experimental  investigations  will  be  conducted  in  "controlled 
non- equilibrium  expansions"  correlated  with  free  stream  conditions  existing 
in  flight  at  the  corresponding  Mach  number.  Configuration  tests  will  be 
conducted  in  the  hypersonic  flow  regime  for  lifting  body  configurations 
characteristic  of  orbital  and  super-orbital  re-entry  flight.  Additional 
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theoretical  investigations  will  be  conducted  in  inviscid  three-dimensional 
hypersonic  flows  and  the  Boltzman  integral  method  will  be  investigated  for 
hypersonic  low  density  flows  over  a  flat  plate.  Incorporation  of  a  six- 
component  force  balance  into  the  lifting  body  hypersonic  investigations  is 
planned  for  the  next  six- month  period. 
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II.  THEORETICAL  INVESTIGATIONS 


A.  Viacoua  and  Invlscid  Non-equilibrium  Gas  Flows 
1 .  Introduction 

The  effects  of  leading  edge  or  nose  bluntneBS  and  boundary  layer 
displacement  have  received  considerable  theoretical  and  experimental 
investigation  for  the  case  of  an  equilibrium  dissociating  gas  or  a  perfect 
gas  with  constant  specific  heats.  The  most  comprehensive  theoretical  and 
experimental  treatment  of  this  problem  has  been  made  recently  by  Cheng, 
Hall,  Golian,  and  Hertzberg.  Since  these  effects  are  most  significant  for 
high  altitudes  and  hypersonic  speeds,  the  equilibrium  assumption  must  be 
examined  in  terms  of  the  flight  spectrum  pertaining  to  hypervelocity  vehi¬ 
cles.  While  an  estimate  of  non-equilibrium  flow  regimes  has  been  made 
by  previous  investigators  in  rather  similar  form,  a  cursory  examination 
of  this  problem  will  be  included  here  for  purposes  of  orientation  and  com¬ 
pleteness  . 

The  species  conservation  equation  for  steady  boundary  layer  flow 
may  be  written 


8a,  8a 

pu  W"  +  P  v  W" 


3 —  (p  a.  Y .  )  " 

8y  'r  i 


W. 


(i) 


where  a  -mass  fraction  concentration  of  the  species  i,  "V~  ^  diffusion 
i  y  i 

velocity  of  species  i,  and  w.  -  net  gas-phase  reaction  rate.  An  estimate 
of  the  degree  of  departure  from  chemical  equilibrium  may,  therefore,  be 
determined  by  the  ratio  of  either  the  convection  or  diffusion  rate  to  that 
for  recombination.  This  ratio  ib  written,  for  convection, 


ir 
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Therefore,  within  the  order  of  magnitude  considered  we  can  take  this  esti¬ 
mate  as 
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»  1  {equilibrium) 


0  (1)  (noil- equilibrium) 


«  1 


("frozen" ) 
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As  pointed  out  by  Cheng,  the  diffusion  term  can  also  be  shown  to  be  gov- 
erned  by  C ..  by  noting  that  (6  /x)  —  £/ux  and  "V--  =  ^  8/ 0y  t  n  ttj. 

The  numerical  solutions  of  Fay  and  Riddell,^  for  the  stagnation  point  prob¬ 
lem,  confirm  the  order-of-magnitude  relation  (2), 

an  estimate  of  the  possible  extent  of  the  non-equilibrium  flow  re¬ 
gime  may  be  observed  in  Figure  1 .  Here,  the  recombination  rate  parameter 

for  "air"  has  been  evaluated  at  the  stagnation  point  assuming  a  rate  law  for 
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oxygen,  recombination  as  proposed  by  Byron.  conclusions  regarding  the 
degree  of  departure  from  equilibrium  downstream  of  the  stagnation  point 
may  be  indicated  by  the  ratio  of  the  local  recombination  rate  parameter  to 
that  at  the  stagnation  point,  namely 


w-d 


(d  u  /dx) 

O  8 

"0  '  /x 
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vhere  it  is  assumed  R  ~  R 


ur  «  U  (hypersonic  small -disturbance 
o  on  7  “ 


assumption)  and  w  -  0  for  the  Byron  rate  law.  Equation  (4)  may  be 
written  .  p 
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where  r  =  radius  of  the  leading  edge.  Under  a  predominant  inviscid  leading 
edge  flow  field,  (5)  becomes 


(RsTfl) 


1/2 


TT 
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(6) 


In  general,  the  blast -wave  local  shock  slope  may  be  written 
d6/dx~  (k  d/x)^^,  Since  k  «  1,  we  may  write  (certainly  within  the  order 
of  magnitude  considered  here)  k^^  (d/x)^  ^  «  d6/dx  which,  within  the 
error  of  the  small-disturbance  assumption,  is  estimated  as  d6/dx~ 

oo 

Therefore,  (6)  may  be  written 


(7) 


For  the  inviscid  leading  edge  dominated  flow  field,  it  can  be  readily  shown 
Tg/Tj  <  M  .  Thus,  for  hypersonic  Bpeeds,  criterion  (7)  is  taken  as 


C1 


<  M 


1/2 

oo 


(8) 


Therefore,  according  to  (8),  if  the  flow  in  the  vicinity  of  the  stagnation 
point  is  not  in  chemical  equilibrium,  then  the  flow  field  downstream  of  the 
leading  edge  will  also  exhibit  a  non-equilibrium  characteristic. 


Similarly,  if  the  boundary  layer  interaction  mechanism  is  dominant, 
(4)  may  be  written 

c,  '* 


l  ®  M2  M  i/2  TT 
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where  XQ0  =  value  of  the  interaction  parameter  based  on  a  1-foot  length. 
In  the  altitude  range  of  200,  000  to  300,  000  ft,  40  <  £  400.  Thus,  for 
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hypersonic  speeds,  Ma<  X^S  M^,  criterion  (9)  i«  taken  as 


Therefore,  except  for  extreme  altitudes,*  if  the  flow  in  the  vicinity  of  the 
stagnation  point  is  not  in  chemical  equilibrium,  then  the  flow  field  down¬ 
stream  of  the  leading  edge  will  also  exhibit  a  non -equilibrium  character¬ 
istic.  Therefore,  Figure  1  may  be  taken  as  an  indication  of  the  degree  of 
departure  from  equilibrium  at  any  point  on  a  sharp  or  blunt  leading  edge 
geometry  in  hypersonic  flow. 

Figure  1  presents  the  recombination  rate  parameter  Cj,  based  on 

stagnation  point  conditions,  for  a  1-foot  nose  or  leading  edge  radius.  Thus, 

for  smaller  radii,  the  same  velocity-altitude  conditions  would  be  expected 

to  show  a  larger  departure  from  equilibrium.  This  would  be  anticipated 

since  the  smaller  the  radius,  the  shorter  the  local  residence  time  for  a 

particle.  By  way  of  reference,  a  typical  boost-glide  vehicle  (wing  loading 

=  25psf,  CL  -  0.5)  and  a  skip  re-entry  vehicle  (L/D  a  0.5,  initial  entry 

angle  =  3  degrees)  are  indicated.  According  to  this  figure,  the  boost-glide 

and  skip  re ••  entry  vehicles  indicate  a  considerable  range  for  which  10”^S  C, 

5  10"  .  Further,  it  is  interesting  to  note  the  time  duration  spent  in  the 

"frozen"  regime.  Specifically,  if  the  frozen  composition  model  is  speci- 
-  3 

fied  by  C.  10  then  the  boost-glide  vehicle  would  experience  approxi¬ 
mately  95  minutes  in  which  the  entire  flow  field  is  frozen.  This  represents 
over  80  per  cent  of  the  total  glide  flight  time  of  this  vehicle.  Examination 
of  the  skip  re-entry  flight  path  indicates  similar  conclusions  as  to  the 
importance  of  non-equilibrium  phenomena. 


It  should  be  remembered,  however,  that  the  stagnation  point  value  of  C.  is 
extremely  small  at  these  altitudes  (C.  10"'  at  300,000  ft),  and  its  subsequent 

multiplication  by  M  ^  (in  order  tolconsider  the  state  of  the  flow  field 
downstream  of  the  leading  edge)  would  still  indicate  a  frozen  composition  and 
thus  should  not  introduce  a  large  error  in  the  assumption  that  the  degree  of 
departure  from  equilibrium  at  any  point  on  the  body  is  indicated  by  the  stag¬ 
nation  point  value  of  C, . 
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With  these  considerations  in  mind,  viscous  and  inviscid  frozen  flow 
studies  of  the  "blast-wave"  and  boundary  layer  interaction  mechanism  for 
blunt  and  sharp  leading  edge  slender  bodies  are  made.  These  studies  are 
discussed  in  Sections  II.  A.  2  and  II.  A.  3,  respectively. 

The  viscous  and  inviscid  theoretical  studies  were  generalized  to 
consider  an  arbitrary  (frozen)  degree  of  dissociation  in  the  ambient  free- 
stream.  This  generalization  was  made  to  permit  analysis  of  the  flow  field 
about  geometries  placed  in  wind  tunnel  nozzles  where  dissociation  non¬ 
equilibrium  is  present  within  the  nozzle  itself.  Since  high  "Mach  number" 
gas  flows  may  be  readily  obtained  in  wind  tunnel  nozzles,  where  dissociation 
non-equilibrium  is  present,  the  "simulation"  capabilities,  in  terms  of  the 
previously  discussed  viscous  and  inviscid  hypersonic  flow  phenomena  of 
"controlled  non-equilibrium  expansions"  are  studied.  These  studies  are 
discussed  in  Section  II.  A.  4. 

2.  Inviscid  Hypersonic  Flow 

The  following  assumptions  are  employed  for  the  inviscid  plane  and 
axiBymmetric  analyses: 

a)  The  gas  composition  is  constant  along  a  streamline  and  given 
by  the  equilibrium  value  immediately  behind  the  shock 
envelope. 

b)  The  flow  process  is  particle-isentropic,  except  in  crossing 
the  shock  envelope. 

c)  The  translation,  rotation,  vibrational  degrees  of  freedom  are 
everywhere  in  equilibrium  with  the  local  gas  temperature. 

d)  The  ambient  free -stream  is  uniform  in  velocity  and  homoge¬ 
neous  in  gas  properties, 

e)  The  free-stream  velocity  is  much  greater  than  the  speed  of 
sound  in  the  ambient  free-stream  atmosphere, 

f)  The  local  flow  deflection  angle  (in  radians),  as  measured 
from  the  undisturbed  flow,  is  much  smaller  than  unity 
everywhere  except  in  the  nose  region. 


iuW 


a.  Sharp-Nosed  Slender  Bodies 

The  differential  equations  describing  the  one-dimensional 
particle -isentropic  flow  for  a  plane  or  axisymmetrlc  flow  may  be  written 
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)v  =  ‘  F  3y 


(10) 


(H) 


(12) 


Since  neither  the  momentum  or  continuity  relations  are  influenced 
by  the  gns  model  chosen,  we  need  only  direct  our  attention  to  the  entropy 
transport  equation.  The  specific  entropy,  $  ,  may  be  obtained  from  the 
relation 
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pd  (p  *)  _ 
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+  T 


dP 

dT 


(1  +  o  )  T] 


pR  (1 


o)  Td  (p 
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(13) 


who  re 

p  r  TT^  +  ITt|rt  +  2pa|“ 

and  p  ,  are  the  coefficients  of  energy  content  for  the  molecular  and 

the  atomic  gas  at  the  given  temperature.  For  the  present  hypersonic  flow 
investigations,  we  can  take  PA  =  5/2;  and  pM>  the  vibrational  energy 
content,  is  taken  as  that  for  a  quantized  linear  oscillator. 
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Relation  (13),  for  frozen  flow,  can  be  simplified  by  noting  that, 

for  air,  we  can  take  ( 1  -  a/  1  +  a)(dp^/dT)  «(a/ 1  +  a)(D/RT^)  whichreduces 
(13)  to 


-  a)  PM  +  4a 


a)0 


M 


in  p  +  (constant) 


(14) 


Therefore,  the  entropy  transport  Eq.  (12)  may  be  written 


B  0  \ 

FT  +  v  F7jPP  =  0  05) 

The  exponent  ( 1  -  a)  Pj^  +  5a/(  1  -  a)  p^  +  4a  -  1  can  be  identified  aB  the 
specific  heat  ratio,  y.  for  a  given  streamline.  Thus  for  the  frozen  flow 
over  a  sharp-nosed  slender  body,  we  notice  the  variation  of  the  shock 
layer  prop*  rlies  are  bounded  by  the  7/5  <  y  <  5/3  (a  =  0,  1,  respectively) 
solutions. 

The  solutions  of  Eqs.  (10),  (11),  and  (15)  can  be  readily  obtained 
for  a  =  constant,  since  assumption  (f)  allows  us  to  take  =  7/2,  The 
utility  of  such  a  solution  must,  therefore,  be  examined.  The  variation  of 
a  with  flow  deflection  angle  as  a  function  of  altitude  and  free-stream 
velocity  is  shown  in  Figure  2.  The  results  presented  in  this  figure  are 
for  a  nondissociated  free-stream,  but  an  arbitrary  degree  of  dissociation 
in  the  ambient  irec-stream  will  have  a  negligible  influence  on  these 
results.  Assumption  (f)  and  Figure  1  indicate  that  for  sharp-nosed  slender 
bodies  we  can  take,  to  good  approximation,  (l-a)Pj^  +  5a/  +  4a  -  1 

~7/5.  The  solution  of  Eqs.  (10),  (11),  and  ( 15)  in  terms  of  a  family  of 


11 


8 

|  0.5 

o 

<n 

a 


h-  300,000  FT,^: 28,000  FT/SEC  ■ 
M  200,000  FTUajsZS.OOO  FT/SEC- 


&  04 
Si 


*  n  ■»  I _ 


/ 

/ 


-h  *  300,000  FT,U^  =20,000  FT/SEC 


h  =  ZOO,  000  FT.1^  >20, pOO  FT/SEC 


DEFLECTION  ANGLE,  8  ~  DEGREES 


Figure  Z.  Degree  of  Dissociation  Behind  Oblique  Shocks  in  Air 
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similar  solutions,  is  readily  obtained  by  means  of  the  transformation^ 
*1  =  y/6(t),  via: 


J£_  = 

Poo 


p— ) p  i") 
\P®  / 


+  *  lr^-1  *  (1) 


rOO 


(16) 


00/ 


£  =  i-Vh) 


i 

3- 


t  2 

together  with  the  Hugoniot  relations  for  strong  shocks,  (6/a^)  »  1, 
(6  =  d6/dt). 


Ii  .  /i_\2 

Poo  ^+1\a  00  / 
P1  y  +  1 

r  = 

Koo  1 

V1  2 

T  =  VTT 


(17) 


Relations  (16)  and  (17)  reduce  the  system  of  Eqs.  (10),  (11),  and 
(15)  to  a  system  of  ordinary  differential  equations  involving  P,  R,  V. 

Since  this  system  of  equations  involves  the  ratio  of  specific  heats,  as  a 
result  of  the  boundary  conditions  (17),  the  influence  of  a  dissociated  free- 
stream  must  be  considered  for  the  value  to  be  employed  in  the  boundary 
condition  (17),  Certainly,  for  a  nondissociated  ambient  free-stream,  the 
Y  =  7/5  value  will  be  applicable.  Consideration  of  the  influence  of  a 
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dissociated  free-stream  on  the  Hugoniot  relations,  over  a  wide  range  of 
chemical  and  thermodynamic  variables  in  air,  has  been  made/’  These 
results  indicate,  as  would  be  anticipated,  that  we  can  bound  the  specific 
heat  ratio  by  7/5  Sy  S.5/3,  Numerical  solutions  for  P(t|),  R.(*l)>"V'l,n) 
have  been  obtained  for  various  conditions  by  several  authors,  '  7  and 

some  of  these  results  are  presented  in  Figure  3.  The  <r  =  2/3,  3/4 
represent,  respectively,  the  two-dimenBional  blunt  leading  edge  flat  plate 
and  boundary  layer  interaction  mechanisms.  Since  the  surface  pressure 
is  dependent  on  the  values  of  P(0)  and  "V(0),  the  relatively  small  influence 
of  the  boundary  conditions,  even  for  the  limiting  cases  (y  =  7/5,  5/3),  can 
be  observed. 

With  the  solutions  for  P(0)  andV(O)  determined,  the  local 
pressure  distribution  is  determined  upon  transformation  of  (16)  to  the 
steady  flow  plane.  To  do  this,  we  make  use  of  energy  conservation  in  the 
following  form: 

U  2  2 

a  D  +  -p-  qD  +  (18) 

00  C  u  i 


Therefore,  except  for  a  highly  dissociated  ambient  gas  Btream  (a  —  1), 
the  ratio  of  the  frozen  pressure  distribution  (p^)  on.  a  sharp-nosed  slender 
body  and  the  corresponding  perfect  gas  value  (P7/5)  should  be  given,  for 
the  same  free-stream  velocity,  to  good  approximation  by 


"7/5 
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00 


rEiSi.  (19) 


00, 


where  Tq^  and  Tqoj  are  the  ambient  free-stream  temperature  for  a 
nondissociated  and  dissociated  gas,  respectively.  Correspondingly,  if  we 


Figure  3.  Variation  of  Pressure  and  Velocity  in  Shock  Laver 


consider  wind  tunnel  application,  we  can  examine  the  ratio  (19)  for  the 
same  "Mach  number."  In  this  cgee,  the  ratio  (19)  becomes 


P& 

P?/5 


(20) 


where  =  frozen  free-stream  velocity  associated  with  non-equilibrium 

nozzle  expansion. 


b.  Blunt-Nosed  Slender  Bodies 

The  above  discussion  is  directed  toward  geometries  for  which 
the  nose  drag  may  be  ignored  in  relation  to  the  drag  of  the  afterbody.  We 
will  now  turn  our  attention  toward  blunt-nosed  bodies  for  which  the  nose  drag 
contributes  a  substantial  portion  of  momentum  and  energy  to  the  flow  field 
downstream  of  the  immediate  vicinity  of  the  leading  edge.  ThiB  condition 
is  the  one  associated  with  the  "blast-wave"  analogy  problem, 

Inclusion  of  the  contribution  of  a  finite  leading  edge  necessitates  the  intro¬ 
duction  of  an  energy  and  momentum  balance  in  the  transverse  flow  field. 

The  change  of  momentum  in  the  axial  direction  may  be  written 

rx  drb  v 

DN  f  pw  ^  ,2irrb)  dx  “ 

•  Jo 


p [<p  -  i>  jf  +  t  v2]  (z*y)v  dy  -  pOT 6  (*6>v  [cp.  -  n 

b 

-  (f  f  *  7  “2  *  2 v2  -  7 (2’ylV  4y  (21) 

’  The  "Mach  number"  is  included  in  quotes  because  of  its  definition  in  a  frozen 
ambient  gas  stream.  This  problem  is  discussed  in  detail  in  Section  11.  A.  4. 
For  our  purpose  here,  we  need  only  consider  It  defined  as  "M"  =  V/"a" 
where  "a"Z  =  (0p/8p)  . 
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where 
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Similarly,  mass  conservation  is  written 
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Substitution  of  relation  (18)  into  (21)  yields 
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We  can  ignore  the  last  term  in  Eq.  (23)  since  we  can  take 


for  most  cases  of  interest. 


In  addition,  for  hypersonic  flow  we  can  take 


whereupon  (23)  can  be  written 
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Finally,  substitution  of  the  continuity  relation  yields 


(24) 
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Equation  (25)  ia  observed  to  reduce  to  the  previous  perfect  gas 
results4,  6‘8,  10for  a  =  a^  =  0,  PM  =  7/2,  (3A  =  5/2. 

Before  continuing  any  further,  it  would  appear  desirable  to  obtain 
an  estimate  of  the  parameter  (a  -  a^)  D/iU^/Z),  the  ratio  of  the  "dissocia¬ 
tion  energy"  to  the  free-stream  kinetic  energy,  for  the  range  of  flight 
conditions  applicable  to  hypervelocity  vehicles,  A  quantitative  discussion 
of  this  parameter  can  be  obtained  by  means  of  available  shock  tables.^ 

These  results  would,  therefore,  correspond  to  the  frce-flight  case  where 
aoo"  ^ata  *or  the  shock  case  are  presented  in  Figure  4, 

while  data  for  the  normal  shock  are  presented  in  Figure  5.*  Thus, 
examination  of  Figures  4  and  5  indicates  that  the  numerical  evaluation  of 
the  third  term  in  Eq,  (25)  must  be  made  for  flow  disturbances  which 
violate  assumption  (f).  We  conclude,  therefore,  since  relation  (25)  is 
dependent  on  the  details  in  the  immediate  vicinity  of  the  leading  edge,  the 
previous  similar  solution  techniques  employed  in  the  perfect  gas 
problem4'^"8,  cannot  be  applicable  for  the  frozen  composition  model. 

Since  the  "dissociation  energy"  term  in  (25)  is  dependent  on  the 
detailed  flow  field  in  the  immediate  vicinity  of  the  leading  edge  (in  particu¬ 
lar,  the  shock  wave  envelope),  a  general  solution  for  arbitrary  nose  geome¬ 
tries  is  not  apparent  to  the  author  at  the  present  time.  General  conclusions 
as  to  the  importance  and  interplay  between  a  blunt  leading  edge  and  the 
afterbody  can  be  obtained,  however,  by  examining  special  cases  in  terms 
of  relation  (25),  If,  initially, we  consider  a  blunt  flat  plate,  an  estimate  of 
the  point  on  the  body  for  which  pw  «  is  given  by 
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RATIO  OF  ‘'DISSOCIATION  ENERGY"  TO  TOTAL  ENTHALPY  ~oD 


Figure  4.  Energy  Distribution  Behind  Oblique  Shock 


Figure  5.  Energy  Distribution  Behind  Normal  Shock 

ft 
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The  above  relation  expresses  the  condition  where  the  energy 
introduced  into  the  transverse  flow  field  due  to  the  presence  of  a  finite 
leading  edge  (nose  drag)  is  just  balanced  by  the  dissociation  energy  which 
is  frozen  at  the  shock  wave  envelope.  Indeed,  since  the  nose.drag  term 
is  independent  of  downstream  distance,  relation  (25a)  clearly  indicates  the 
possibility  of  obtaining  values  of  P^/p^  <  I.  Substituting  for  the  nose 
drag  in  (25a),  we  obtain 
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If  we  consider  leading  edge  geometries  (e.g.  ,  flat,  wedge,  cone, 
leading  edges)  for  which  (ct  «  a(j))D/(U^2)  =  constant  for  dimensions  of 
th*>  order  of  d,  then  {<’.5l»)  c<tti  be  written 
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where  a  ;  mass  fraittun  dissociated  behind  shock,  wave  in  region 
bi-twccn  6  -  +_  d/d.  A  plot  of  the  term 
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the  difference  between  the  nose  drag  and  the  "dissociation  energy"  term,  „ 
is  presented  as  a  function  of  nose  geometry,  free-flight  velocity,  and 
altitude  in  Figure  6.  The  results  for  the  circular  nose  are  based  upon 
the  normal  shock  value  of  the  dissociation  energy  term.  According  to 
Figure  6,  for  flight  velocities  greater  than  20,000  fps,  relation  (25a) 
should  be  satisfied  for  6/d  1,  Since  the  "negative  running"  charac¬ 

teristics  from  the  region  6/d  <J  1  will  intersect  the  surface  of  the  body 
very  near  the  leading  edge,  the  condition  pw  «  p ^  should  occur  very 
close  to  the  leading  edge  for  a  flat  afterbody.  Accordingly,  the  "blast- 
wave"  mechanism  usually  associated  with  blunt  nosed  slender  bodieB  in 
hypersonic  flow  will  be  considerably  reduced  for  altitudes  (£200,000  ft) 
where  the  "frozen"  composition  model  is  applicable.  In  thiB  case,  most 
leading  edge  geometries  can  be  considered  "sharp"  nosed  with  the  boundary 
layer  interaction  mechanism  determining  the  downstream  flow  field  for 
flat  afterbody  geometries  (a  detailed'discussion  of  the  frozen  boundary 
layer  interaction  problem  is  presented  in  Section  II.  A.  3). 

In  order  to  obtain  numerical  verification  of  the  preceding  con¬ 
clusions,  a  characteristics  solution  of  the  frozen  hypersonic  flow  about 
a  blunt  leading  edge  slender  body  was  made.  Employing  the  "sonic-wedge" 
model  of  Bertram,^  the  frozen  and  equilibrium  flow  models  have  been 
numerically  examined  for  =  28,000  fps,  h  =  300,000  ft.  To  be  sure, 
viscous  effects  must  be  considered  in  any  accurate  determination  of  the 
flow  field  at  these  flight  conditions.  However,  this  consideration  is  not 
pertinent  to  the  objective  of  this  particular  investigation.  The  characteris¬ 
tics  net  for  the  "sonic-wedge"  problem  is  shown  in  Figure  7.  The  large 
region  downstream  of  the  leading  edge  influenced  by  the  gas  contained 
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FLAT  NOSE  ~  300,000  FT 

■  CIRCULAR  NOSE  ~  200,000  FT 

■  CIRCULAR  NOSE  ~  300,000  FT 
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Figure  7.  Characteristics  Diagram  for  "Frozen"  Flow 
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between  the  centerline  of  the  wedge  and  the  "limiting  characteristic"  can 
be  observed  in  this  figure.  This  iB  a  result  of  the  "sonic-wedge"  geometry 
which  should  provide  a  larger  downstream  influence  than  that  associated 
with  flat  or  circular  noses. 

For  this  particular  case,  . 
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and  relation  (25c)  takes  the  form 
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which,  observing  the  value  of  al)/  U^/ 2  (Figure  4)  associated  with  the 
"sonic-wedge"  leading  edge  region,  requires 


Numerical  integration  of  the  right-hand  term  in  (26)  can  be 
obtained  from  the  wave  diagram  shown  in  Figure  7.  The  results  of  this 
integration  arc  presented  in  Figure  8,  It  can  be  observed,  from  this 
figure,  relation  (27)  iB  satisfied  for  6/ d  =s  0,8,  The  negative  characteris¬ 
tic  at  6/d  a 0.8  is  observed  (wave  diagram  Figure  7)  to  intersect  the  after¬ 
body  of  the  "sonic-wedge"  at  x/ds  5.5.  This  downstream  distance  should, 
therefore,  provide  an  estimate  of  the  point  where  pap. 


'  The  "sonic -wedge"  geometry  was  chosen,  of  course,  in  order  to  avoid  the 
difficult  detached  shock  problem  associated  with  the  flat  or  circular 
nose  geometries. 
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The  surface  pressure  distribution  as  a  function  of  downstream 
distance  in  leading  edge  diameters,  for  both  the  equilibrium  and  frozen 
flow  models,  is  presented  in  Figure  9.  The  agreement  of  the  estimated 
point  at  which  Pw  k  from  relation  (26),  with  the  characteristics  solu¬ 
tion,  can  be  observed  in  this  figure.  In  addition,  the  previously  discussed 
possibility  of  obtaining  values  of  Pw/P00<  1  can  be  observed  in  this 
figure.  This  is  a  result  of  the  condition 


The  large  error  introduced  with  the  vise  of  an  equilibrium 
pressure  distribution  ran  also  be  readily  observed  in  Figure  U  should 
•bo  remembered  at  Ibis  point,  as  previously  discussed,  the  "soni  c -wedge" 
geometry  should  provide  a  maximum  downstream  influence  of  the  loading 
edge.  Detached  shock  leading  edge  geometries,  according  to  (25c),  should 
Pw~  ‘  loser  to  the  lending  edge  than  that  for  the  "soni c -wedge" 
geometry,  thus,  greatly  diminishing  and  eliminating,  for  most  practical 
purposes,  the  "blast -wave"  mechanism. 


The  corresponding  frozen  and  equilibrium  temperature  and 
electron  concentration  distribution  at  the  outer  edge  of  the  boundary  are 
slH>„r>  in  Figures  10  and  11,  respectively.  The  frozen  electron  concentra¬ 
tion  is  determined  assuming  the  electrons  to  be  in  equilibrium  with  the 
local  frozen  gas  temperature.  Finally,  Figure  12  presents  the  frozen  and 
equilibrium  skin  friction  and  heal  transfer  coefficients  as  a  function  of 
downstream  distance.  Both  the  skin  friction  and  heat  transfer  parameters 
are  normalized  with  respect  to  the  wall  velocity  and  enthalpy  gradients, 
respectively,  in  order  to  allow  incorporation  of  the  influence  of  finite 
surface  catalysis. 
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Figure  10.  Temperature  Distribution  About  Blunt  Body 
in  Hypersonic  Flow 
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Figure  12.  Local  Skin  Friction 
About  Blunt  Body 
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The  above  discussion  neglected  the  presence  of  an  afterbody. 

According  to  (25),  the  afterbody  1b  observed  to  govern  the  down-stream 
flow  field  when 

dr 
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The  absolute  value  signs  are  employed  in  (29)  to  include  the  condition 
in  which  rase  the  local  pressure  on  the  surface  of  the  body  will  be  less 
than  that  associated  with  a  "sharp"  nosed  geometry  with  the  same  after¬ 
body.  An  estimate  of  the  inequality  (29)  can  be  obtained  if  we  consider 
blunt  nosed  cones  or  wedges,  namely: 
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which,  for  order-of-magnitude  purposes,  may  be  taken  as 
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Relation  (31)  indicates  that  the  downstream  surface  pressure 
on  a  blunt  nosed  cone  or  wedge  can  be  lower  than  the  corresponding  "sharp" 
nosed  value  for  a  considerable  distance  downstream  of  the  leading  edge. 

The  exact  downstream  distance  for  which  this  is  the  case  is,  according  to 
(30),  dependent  on  the  shock  wave  envelope  and  the  altitude  and  velocity 
Q.e.  ,  (a  -  o^)  D/IT^/Z^  of  the  vehicle. 


3.  Viscous  Hypersonic  Flow 

The  following  assumptions  are  made  for  the  viscous  flow  studies: 

a)  The  flow  field  in  the  shock  layer  is  assumed  to  be  divisible 
into  a  distinct  boundary  layer  and  an  inviscid  flow  external 
to  the  boundary  layer. 

b)  The  inviscid  fluw  field  is  assumed  to  be  governed  by  the 
inviscid  flow  assumptions  listed  in  Section  IT.  A.  2,  (a) 
through  (f). 


With  these  assumptions  in  mind,  we  will  now  direct  our  attention 
to  the  "frozen"  hypersonic  boundary  layer  interaction  problem. 

a .  "Sha r p"  Leading  Edge  Flat  Plate 

Employing  the  Howarth  transformation,  the  boundary  layer 
displacement  thickness,  f>  ,  may  be  written 

,  Vp  "  U6, 
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where  s'  =  J  p^p^u^ds.  Since  the  gas  composition  in  the  boundary  layer 
is  frozen  at  the  value  at  the  outer  edge  of  the  boundary  layer,  we  can  write 
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where  H/H  =«,£'*>  u/u..  We  can  further  write 
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Downstream  of  the  leading  edge,  the  temperature  is  sufficiently  low  so  that 
we  may  take 
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Therefore,  the  integral  in  (32)  may  be  written 
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For  the  lower  temperatures  associated  with  downstream  of  the  leading 
edge,  we  can  take 


Therefore,  relation  (32)  can  be  written 
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If  we  assume  inviscid  pressure  distributions  of  the  form  of  relations  (16) 
and  (17)  and  note  that  for  the  boundary  layer  interaction  problem,  we  can 


write 


w  -  3t 


bH 


n-Bt11'1  =  VjV(O) 


and,  substituting  for  Vj,  we  obtain 


y  +  1  B 


where  we  have  taken  6  =  At^  and  =  Bt^ .  Identifying  the  boundary  layer 
displacement  thickness  with  r^t)  we  can  write  (40)  as 
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(43) 


The  above  equation  can  be  treated  analytically  by  considering  the  limiting 
cases  of  large  and  small  Og  (i . e . ,  in  order  to  let  1  +  ag  ^  constant).  In 
this  case,  we  may  remove  the  term  (1  +  Qg)  from  the  integral.  In  lieu  of 
this  relation,  (43)  must  be  solved  numerically  for  a  prescribed  Og  =  Ug(s) 
which  then  introduces  additional  complication  into  I  (P,  Sw,  a)  since 
OgD/H^  4  constant.  Therefore,  for  constant,  we  can  write 
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For  either  the  hypersonic  nozzle  problem  or  the  free-flight  case,  we  can 
take 
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which  simplifies  (44)  to 
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Making  use  of  relation  (42)  and  noting  &  -  3/4  for  the  boundary  layer 
interaction  problem,  we  can  solve  for  A: 
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Substitution  of  (47)  into  (16)  yields  the  surface  pressure  in  the  form 
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where  X.  5  (X./Re  }^2,  Re  =  p,U  s/p  .  Since  we  are  concerned 
with  the  downstream  problem,  we  can  take  y  =  7/5,  as  in  the  previous 
study,  for  all  cases  except  those  where  o  "*  1 ,  The  only  other  term 
which  would  provide  any  numerical  difficulty  would  be  /V(0). 

Since  the  previous  remarks  bound  the  value  of  this  parameter  between 
the  7/5  i  v  £5/3  value,  an  inspection  of  Figure  3  indicates 
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for  y  =  5/3 


(49) 


Therefore,  the  y  =  7/5  value  should  numerically  introduce  a  small  error 
into  relation  (48).  Further,  we  should  also  note  we  can  take 
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Numerical  evaluation  of  I((3,  S  ,  n)  will  complete  the  nutner - 
ical  requirements  for  relation  (48).  In  general,  numerical  evaluation  of 
I(fi,  S  ,  n)  can  only  be  made  for  a  given  set  of  altitude -velocity  conditions 
since  g  and  f,  among  other  things,  are  functions  of  Prandll  number  and 
thus  vary  due  to  the  dissipation  term  in  the  energy  equation  according  to 
N.  Kemp.  If,  however,  we  ignore  this  variation  for  the  range  of  altitude 
and  velocity  we  are  concerned  with  and,  also,  consider  g  and  f  independ¬ 
ent  of  n^D/H^  (an  assumption  which  is  certainly  open  to  question  for  large 
values  of  q^D/H^),  we  can  ttien  numerically  evaluate  I(p,  S^,  n)  as  a 
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function  of  agD/H^.  Therefore,  l(p,  Sw,  a)  is  evaluated  numerically 
subject  to  the  following  boundary  conditions: 


f<0)  =  0,f'(0)  =  0,  f'(oo)  =  1 
cig  D 

S(0}  =  jq —  .  g(oo)  =  1 


q'{0)  =  0 

where  the  primed  condition  denotes  d/dtp 

The  above  conditions  imply  a  highly  cooled,  "noncatalytic"  wall. 
The  numerical  results  are  presented  in  Figure  13.  As  would  be  anticipated, 
for  OgD/H  =  1,  I  ((3,  Sw,  n)  ~  0.  Namely,  if  all  the  energy  is  frozen  in 
the  dissociation  mode  n  viscous  dissipation  mechanism,  the  boundary  layer, 
is  not  required.  To  be  sure,  as  «gD/H^-*l,  tho  assumptions  (i.e., 
hypersonic  small  disturbance,  invariance  of  g  and  f,  neglect  of  ionization) 
employed  in  this  analysis  are  no  longer  valid;  and  therefore,  the  region  in 
the  vicinity  of  I((l,  S^,  Q)  ~  1  should  not  be  considered  too  seriously. 

Before  any  numerical  evaluation  of  (48)  is  made,  one  must  nat¬ 
urally  inquire  as  to  the  range  of  values  of  tig  for  sharp  leading  edges  which 
are  associated  with  realistic  flight  conditions.  This  consideration  can  be 
made  in  a  simple  manner;  namely,  consider  the  flow  field  geometry  in 
Figure  14.  Mass  conservation  requires 

p  U  h .  =  pudy  (52) 
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which,  under  the  Howarth  transformation,  can  be  written 
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Figure  13.  "Displacement  Thickness  Integral"  for  a 
"Non-Catalytic"  Surface 
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Substitution  of  relation  (48)  into  (52a)  and  differentiating  we  obtain 
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If  we  substitute  for  a  =  0,  v  =  ?/5,  we  obtain 

cn 


For  practical  flight  conditions,  for  which  the  division  of  a  distinct  boundary 
layer  and  inviscid  flow  is  possible,  ^  ^  *  ].  According  to  (54),  this 

condition  implies  flow  deflections  of  less  than  25  degrees  which,  upon  exami¬ 
nation  of  Figure  4,  indicate  values  of  alJ/lI^  <  0.1.  Therefore,  "freezing" 
of  the  dissociation  energy  should  have  a  small  influence  on  the  boundary 
layer  interaction  problem  for  a  "sharp"  flat  plate  Specifically,  the  frozen 
boundary  layer  interaction  on  a  "sharp"  flat  plate  should  be  given,  to  good 
approximation,  by  the  perfect  gas  y  =  7/5  prediction. 

b.  "Blunt"  Leading  Edge  Flat  Plate 

The  boundary  layer  interaction  problem  for  the  blunt  leading 
edge  flat  plate  is  defined  as  the  condition  when  the  streamlines  at  the 


outer  edge  of  the  boundary  layer  are  those  which  passed  through  the  near¬ 
normal  shock  in  the  vicinity  of  the  leading  edge.  In  contrast  to  the  investi¬ 
gations  in  Reference  li  the  decreased  importance  of  the  frozen  blast-wave 
mechanism,  discussed  in  II.  A.  2.  b,  allows  examination  of  this  problem  with 
neglect  of  the  downstream  influence  of  the  leading  edge.  In  addition,  since 

we  are  net  concerned  here  with  the  low-Reynolds  number  problem,  we  will 

1  3 

ignore  vorticity  corrections  such  as  those  considered  by  Lees. 

With  these  remarks  in  mind,  we  need  simply  evaluate  the  con¬ 
ditions  at  the  outer  edge  of  the  boundary  layer  for  incorporation  into  rela¬ 
tion  (48).  In  this  case  ,  is  equal  to  the  value  behind  a  normal  shock  for 
the  given  free-3tream  velocity  and  altitude.  The  value  of  I(p,  S^,  a)  is 
equal  to  the  value  obtained  from  Figure  13  with  a^D/H^  equal  to  that  behind 
a  normal  shock  for  the  given  free -stream  velocity  and  altitude.  Finally,  the 
velocity  at  the  outer  edge  of  the  boundary  layer  is  obtained  by 


with  u  =  value  behind  normal  shock.  With  these  considerations  in  mind, 
the  frozen  boundary  layer-induced  pressure  for  a  "blunt"  and  "sharp"  flat 
plate  is  presented  in  Figure  15.  The  "catalytic"  curve  is  taken  as  that  for 
which  1(P,  S^,  a)  is  equal  to  the  value  for  a^D/H^  =  0.  The  large  differ¬ 
ences  which  are  observed  in  thiB  figure  are  due  mainly  to  variations  of 
I(P,  S  ,  a)  with  a  d/H  The  relatively  weak  dependence  on  free-stream 
velocity  can  also  be  observed  in  this  figure.  As  would  be  anticipated,  the 
"sharp"  leading  edge-induced  pressure  is  given  closely  by  the  perfect  gas 
value  (i.e.,  p/p^a  0.16 x)-  The  "blunt"  leading  edge  "catalytic"  wall 
value  is  slightly  higher  as  a  result  of  the  term  [(1  +  a^)//l  +  “^l]  in 
relation  (48). 


i 
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Figure  I  5.  F’rozen  Boundary  Layer  Interaction  on  a 
"Blunt"  and  "Sharp"  Flat  Plate 
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4.  "Controlled  Non -Equilibrium  Expansions" 

One  of  the  most  serious  problems  associated  with  the  development 
of  high  enthalpy-low  density  hypervelocity  facilities  is  the  presence  of  dis¬ 
sociation  non-equilibrium  within  the  expansion  nozzle.  Indeed,  the  stag¬ 
nation  conditions  required  for  experimentally  investigating  dissociative 
non-equilibrium  phenomena  can  give  rise  to  large  departures  from  equi¬ 
librium  v/ithln  the  nozzle  itself.  This  problem  has  received  considerable 
theoretical  investigation  for  single  component  systems. In  addition, 

a  recent  experimental  investigation  of  this  problem,  for  a  multi-component 

17 

system  (air),  has  been  made  by  Nagamatsu,  et  al.  Since  the  analysis  of 
the  non-equilibrium  nozzle  flow  problem  has  been  discussed  in  detail  in 
References  5,  14-16,  we  shall  not  discuss  it  further.  Rather,  we  Bhall 
employ  the  results  in  the  present  discussion. 

If  we  define  a  Mach  number  "M"  such  that 


"M" 


V^p/9p,s 


then,  for  a  dissociating  gas  we  can  write 
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It  can  be  easily  shown  that  (57),  for  a  constant  gas  composition,  can  be 
written 
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where  we  have  neglected  dfJj^/dT  in  comparison  to  the  remaining  terms. 
With  "Mach  number"  in  a  frozen  gas  stream  defined  by  (56)  and  (58),  we 


I*  f » '  *1  <"'<  mnKHnwmtmmm 


have  calculated  "M"  as  a  function  of  a  nozzle  expansion  ratio  assuming  a 

3 

rate  law  for  "air"  recombination  to  be  given  by  that  proposed  by  Byron  for 
oxygen  recombination,  namely  k^(T)  =  2.  16  x  10^  cm^/mole^  sec. 

(T  in  °K). 

To  be  sure,  the  complex  kinetics  associated  with  the  coupled  reac¬ 
tions  in  air  cannot  be  represented  by  a  simple  temperature  dependence; 
however,  the  relatively  weak  influence  of  k^(T)  on  the  non -equilibrium  nozzle 
expansion  process  will  lend  some  degree  of  generality  to  the  results.  In 
addition,  the  present  results  can  easily  be  modified  to  incorporate  more 
accurate  rate  data  when  it  becomes  available. 

Figure  16  presents  "M"  as  a  function  of  nozzle  expansion  ratio  for 
a  stagnation  temperature  of  8000°K  and  for  stagnation  pressures  (PQ)  of  10, 

100,  and  1000  atmospheres.  The  calculations  were  made  for  a  conical  noz¬ 
zle  with  a  7. 5-degree  half-angle  and  a  throat  diameter  of  0,1  inch.  How¬ 
ever,  since  the  nozzle  dimensions  and  geometry  have  a  weak  influence  on 

5,15-17,  .  ,  .  ,  . 

the  non-cquuihrium  process,  the  results  can  be  conuiueieu  quite 

general.  These  stagnation  pressures  would  correspond,  respectively,  to 
arc-jet,  shock  tunnel,  spark-heated  (hot  shot)  operation.  The  difference 
between  the  equilibrium  and  non-equilibrium  "Mach  number"  iB  due  to  the 
large  difference  in  static  temperature  between  the  equilibrium  and  non- 
equilibrium  expansion  process.  Figure  16,  therefore,  indicates  the  attrac¬ 
tive  possibilities  for  generating  High  "Mach  number"  gas  flows  by  means  of 
"controlled  non-equilibrium  expansions"  (using  stagnation  pressure  varia¬ 
tions).  The  main  consideration  here,  of  course,  is  the  utility  (in  terms  of 
aerodynamic  simulation)  of  this  frozen  composition  high  "Mach  number"  gas 
stream.  A  measure  of  the  dissociation  energy  frozen  in  the  nozzle  is  shown 
In  Figure  17.  In  this  figure  the  ratio  of  the  frozen  dissociation  energy  in  the 
nozzle  to  the  stagnation  enthalpy  Is  presented  as  a  function  of  stagnation  pres¬ 
sure.  The  trade -offs  between  "Mach  number"  and  frozen  dissociation  energy 
can,  therefore,  be  observed  in  Figures  16  and  17.  The  simulation  capabili¬ 
ties,  in  terms  of  "M"  and  x  >  for  a  stagnation  pressure  of  10  atmospheres 
are  shown  in  Figure  18.  Indeed,  the  possibility  of  concurrent  "Mach  number" 
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and  Reynolds  number  simulation  can  be  observed  from  this  figure,  These 
parameters  are  indicated  for  the  previously  discussed  Byron  rate  law  and 
a  rate  law  two  orders  of  magnitude  "faster"  (i.c.,  Byron  rate  law  x  102). 
The  previously  discussed  weak  dependence  on  kj^(T)  can  be  observed  in 
this  figure.  With  these  "performance  charts"  for  orientation,  we  shall 
now  turn  our  attention  toward  a  discussion  of  the  simulation  capabilities  of 
high  "Mach  number"  frozen  composition  gas  flows. 

If  first  we  consider  the  simulation  capabilities  in  terms  of  the 
similitude  for  slender  sharp-nosed  bodies  establi  shed  by  Cheng,  we  can 
write  the  surface  pressure  similitude  as 


rr  {  M  r, 

oc 


P  _ .  P  ,  Q  ) 
1  oo  '  oo  on 


(59) 


where  p  ,  p  ,  a  expresses  the  free-stream  atmosphere  simulation  re- 

quirement.  Since  the  surface  pressure  has  previously  indicated1^  a  weak 

dependence  on  p  ,  p  ,  a  ,  we  can  examine  this  similitude,  in  terms  of 

the  Hugoniot  relations,  for  the  non-equilibrium  nozzle  flows  indicated  in 

Figures  16  to  18  where,  indeed,  p  ,  Pu),  n(jj  are  not  properly  simulated. 

These  results  are  presented  in  Figure  19.  In  this  figure  the  pressure  ratio 

across  an  oblique  shock  is  plotted  against  the  similarity  parameter  "M"  t 

(t  wedge  angle)  for  the  "Mat  h  number"  range  6.7  "M"  62.5.  The 

agreement  with  the  similitude  curve  for  a  nondisHoclated  free-slrcanvat- 

2  2 

mo  s  phere,  given  to  go  od  approximation  by  p/p^  -  y(y+1)/2"M"  t  ,  can 

he  observed  in  this  figure.  Similar  agreement  with  the  lift  curve  similitude 

can  be  observed  in  Figure  20.  Numerical  results  for  the  skin  friction  and 

heat  transfer  coefficient  '.re  presented  in  Figure  21.  The  solid  curve  labeled 

"free-flight"  represents  a  correlation  of  the  Hugoniot  relations  for  the  range 

10,  000  <  V  <  28,000  fps  at  h  -  300,000  ft,  The  agreement  of  the  frozen 

skin  friction  and  heat  transfer  coefficients  with  the  "free-flight"  values 

from  10  £  "M"  25  can  be  observed  in  this  figure,  The  results  are  pre- 

2 

sented  for  both  the  Byron  rate  law  and  the  Byron  rate  x  10  .  The  corres¬ 
ponding  lack  of  simulation  for  " M"  =  35.5,  -15,  62.5  can  also  be  observed. 
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Thus,  for  the  simplified  flow  mech&niem  considered,  meaningful  pressure, 
lift,  skin  friction,  heat  transfer  simulation  for  "M"  5  25  appears  possible. 

We  will  now  turn  our  attention  toward  consideration  of  more  "complicated" 
flow  mechanisms. 

In  order  to  examine  the  blunt  leading  edge  slender  body  problem,  nu¬ 
merical  solutions  of  the  flow  about  a  "sonic-wedge"  model  at  "M"  =  30.  3 
were  made.  The  value  of  "M"  =  30.  3  was  selected  in  order  to  allow  a 
comparison  with  the  previous  results  presented  in  Figures  9  and  10  for 
V  =  28,  000  fps  and  h  =  300,000  ft,  corresponding  to  "M"  =  30.3.  Numer¬ 
ical  results  were  obtained  for  P  =  10  and  100  atmospheres,  T  =  8000°K. 

0  4  ®  3 

"M"  =  30.  3  corresponds  to  an  expansion  ratio  of  4.  4  x  10  and  6  x  10  for 

P  -  10  and  100  atmospheres,  respectively.  The  wedge  angle  corresponding 
to  the  "  sonic-wedge"  model  for  100  and  10  atmospheres  were,  respectively, 

52  degrees  and  57  degrees.  The  numerical  results  are  presented  in  Fig¬ 
ure  22.  The  calculations  were  only  carried  out  for  the  wedge  shock  and  the 
expansion  around  the  wedge  shoulder  because  of  the  computational  time  in¬ 
volved  in  continuing  the  calculation  further  by  characteristics.  A  compari¬ 
son  of  Figure  22  and  Figures  9  and  10  indicates  the  simulation  for  the 
stagnation  conditions  chosen.  It  is  felt  the  major  difference  can  be  ascribed 
to  the  lower  enthalpy  (i.  e.,  lower  than  that  corresponding  to  28,  000  fps)  for 
the  given  stagnation  conditions  with  the  associated  decrease  in  the  parameter 
oD/H  behind  the  shock  wave  envelope.  This  conjecture  is  currently,  being 
evaluated  numerically.  While  this  simulation  is  far  from  adequate,  it  should 
be  examined  In  terms  of  that  obtained  by  meanB  of  existing  facilities.  Namely, 
the  differences  in  pressure  and  temperature  indicated  by  a  comparison  of 
Figure  22  and  Figures  9  and  10  should  be  compared  to  those  differences  asso¬ 
ciated  with  the  corresponding  perfect  gas  values. 

The  simulation  capabilities  of  high  "Mach  number"  frozen  gas  flows, 
in  terms  of  the  boundary  layer  interaction  problem,  are  presented  in  Fig¬ 
ure  23.  The  free-flight  curves  are  obtained  with  =  0  and  employing 
relation  (48). 
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Figure  22.  Controlled  Non-Equilibrium 
Expansion  "Simulation" 


Finally,  it  should  be  noted  that  this  investigation  is  based  upon 
employing  a  "froaen"  composition  model  for  the  region  immediately  behind 
the  shock  wave  envelope.  It  has  already  been  indicated  that  this  model  is 
bounded  at  lower  altitudes  by  those  densities  which  are  high  enough  so  as 
to  yield  particle  rosidence  times  of  the  same  order  as  those  associated 
with  the  chemical  rates.  Certainly,  this  model  is  bounded  at  the  higher 
altitudes  by  the  forward  problem,  namely,  the  time  to  dissociate  in  the  re¬ 
action  zone  behind  the  shock  envelope  (cf.  Reference  5).  In  addition,  low- 
Reynolde  number  effects  will  .also  serve  to  introduce  an  upper  limit 
in  altitude. 

5 .  Concluding  Remarks 

A  theoretical  study  has  been  made  of  the  "frozen"  viscous  and 
inviscid  hypersonic  flow  about  blunt  and  sharp-nosed  slender  bodies. 

The  inviscid  investigation  indicates,  as  would  be  anticipated, 
"freezing"  of  the  dissociation  energy  in  the  region  Immediately  behind  the 
shock  wave  envelope  has  a  small  influence  on  the  sharp-noaed  slender  body 
problem.  In  this  case,  the  flow  field  should  be  given,  to  good  approximation, 
by  the  perfect  gas  (y  =  7/5)  prediction. 

Conversely,  the  frozen  flow  over  blunt-noBed  slender  bodies  exhibits 
large  departures  from  either  the  equilibrium  or  perfect  gas  predictions. 

An  energy  balance  in  the  transverse  flow  field  indicates  this  large  differ¬ 
ence  is  primarily  due  to  the  condition  where  the  energy  introduced  into 
the  transverse  How  field  by  the  presence  of  a  finite  leading  edge  (nose 
drag)  is  offset  by  the  corresponding  dissociation  energy  which  is  frozen 
at  the  shock  wave  envelope,  In  addition,  since  the  nose  drag  term  is 
independent  of  downstream  distance,  the  energy  distribution,  (i.  e.  ,  the 
large  amount  of  energy  frozen  in  dissociation  indicates)  the  possibility  of 
obtaining  surface  pressures  lower  than  the  ambient  free-stream  value. 

These  conclusions  are  verified  numerically  by  a  characteristics  solution 
for  the  frozen  and  equilibrium  hypersonic  flow  about  a  "sonic-wedge”  Blab. 
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The  viscous  investigations  take  the  form  of  an  investigation  of  the 
frozen  boundary  layer  interaction  problem.  As  in  the  sharp-nosed 
inviscid  studies,  the  frozen  composition  model  has  a  small  influence  on 
the  sharp  leading  edge  problem  for  values  of  the  hypersonic  interaction 
parameter,  y  ,  less  than  1000.  The  large  differences,  compared  to  the 
-y  =  7  /5  value,  for  the  blunt  leading  edge  case  are  primarily  due  to  the 
decrease  in  the  ’’boundary  layer  displacement  integral”  for  a  noncatalytic 
wall. 

The  viscous  and  inviscid  theoretical  studies  were  generalized  to 
consider  an  arbitrary  (frozen)  degree  of  dissociation  in  the  ambient  free- 
stream.  This  generalization  was  made  to  permit  analysis  of  the  flow  field 
about  geometries  placed  in  wind  tunnel  nozzles  where  dissociation  non- 
equilibrium  is  present  within  the  nozzle  itself.  Since  high  "Mach  number" 
gas  flows  may  be  readily  obtained  in  wind  tunnel  nozzles  where  dissocia¬ 
tion  non-equilibrium  is  present,  the  simulation  capabilities  in  terms  of 
the  previously  discussed  viscous  and  inviscid  hypersonic  flow  phenomena, 
by  means  of  "controlled  non-equilibrium  expansions"  were  studied.  The 
results  indicate  attractive  simulation  possibilities  which  are  currently 
under  further  theoretical  and  experimental  investigation. 
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SYMBOLS  (Section  II.  A) 

a  speed  of  sound 

A  area 

Cf  skin  friction  coefficient 

C  recombination  rate  parameter 

drag  coefficient 
lift  coefficient 

c  heat  transfer  coefficient 

q 

d  leading  edge  diameter 

D  dissociation  energy 

Dj^  nose  drag 

e  internal  energy 

f'  boundary  layer  velocity  profile,  f'  =  u/u^ 

g  boundary  layer  enthalpy  profile,  g  -  H/H^ 

h  Btatic  enthalpy 

H!  total  enthalpy 

k  nose  drag  coefficient,  k  s  D^/l/2  d/2  (ir  d/2)V 

specific  rate  constant,  recombination 
M  molecular  weight 

p  static  pressure 

stagnation  pressure 
Pr  Prandtl  number 

q  heat  transfer  rate 

R  universal  gas  constant 
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SYMBOLS  (continued) 


local  radius  of  the  body 

Reynolds  number  per  length,  Re  =  ,">u/(ji 

f* 

defined  by,  s'  =  J  pgPgUgds 
entropy 


static  temperature 

axial  velocity  component 

free-stream  velocity 

normal  velocity  component 

velocity,  V  =  \/i+v^ 

diffusion  velocity 

axial  distance 

distance  normal  to  surface 

mass  fraction  dissociated 

dimensionless  enthalpy  defined  by  h  =  Pp/p 

3  1/2 

hypersonic  interaction  parameter,  X  =  M‘  (X/Re^) 

nondimensional  distance  between  shock  wave  envelope  and  body 
surface,  r|  =  y/6 

shock  wave  envelope 

boundary  layer  displacement  thickness 

d6  /dt 

ratio  of  specific  heats 

constant  in  linear  viscosity-temperature  law 

viscos  ity 
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SYMBOLS  (continued) 

v  0,  1,  respectively,  for  two-dimensional  and  axisymmetric  flow 
«r  shock  wave  power  law,  6<—>tff 
p  density 

t  fineness  ratio 

6  wedge  angle 

w  reaction  rate 

Subscripts 

A  atomic  species  translational  energy  term 

E  conditions  at  equilibrium 

f  frozen  composition 

i  refers  to  i1^  atomic  chemical  species 

M  molecular  species  energy  term 
8  evaluated  at  stagnation  point 

w  evaluated  at  the  surface 

6  evaluated  at  the  outer  edge  of  the  boundary  layer 
oo  ambient  free-stream 
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B.  Non-EquUlbrtum  Flow  About  a  Blunt  Body 
1.  Introduction 

An  inverse  method  for  analyzing  compressible  flow  about  a  blunt 

body  was  first  suggested  by  Lin  and  Rubinov.  ^  This  method  starts  with 

an  assumed  form  of  shock  shape.  The  flow  field  behind  the  shock  and  the 

body  shape  are  determined  by  numerical  analyses.  Further  investigations 

w*>re  carried  out  by  Lin  and  Shen,^*  Cabannes,^  Van  Dyke,^  Zlotnick 
24  25 

and  Newman,  Mangier  and  Evans.  Other  refined  analyses  classified 

26 

as  inverse  methods  were  also  given  by  Garabedian  and  Lieberstein,  and 

27 

Vaglio- Laurin  and  Ferri.  A  detailed  review  of  these  analyses  may  be 

found  in  Ref.  28.  More  recently  a  comprehensive  review  of  advances  in 

29 

non-equilibrium  dissociating  gas  dynamics  was  given  by  Li.  The  case 
of  a  reacting  and  relaxing  steady  flow  around  a  blunt  body  was  treated  by 
Lick^  who  used  the  inverse  method  of  analysis.  Some  rather  interesting 
results  were  presented.  However,  these  results  are  only  of  limited  quali¬ 
tative  value  because  of  the  assumptions  used  in  the  analysis.  The  two 
important  assumptions  employed  are  listed  as  follows: 

a)  Air  is  considered  as  a  binary  mixture  composed  of  atoms 
and  molecules.  Therefore,  the  distinction  among  NO,  O,, 
N.,,  and  the  distinction  between  O  and  N,  are  not  taken 
into  account. 

b)  A  parameter,  \  =  (R  Ip  V  )(ui7/p),  is  introduced  in  the 

analysis  where  Rq  is  the  radius  of  curvature  of  the  nose  at 
the  stagnation  point:  represents  the  mass  rate  of  produc¬ 

tion  of  the  molecular  species  by  chemical  reactions. 
Numerical  computations  were  carried  out  for  the  flow  fields 
of  constant  values  of  x-  When  the  value  of  \  is  large,  for 
example  \  =  100,  the  solution  corresponds  closely  to  that 
of  equilibrium  dissociation.  When  the  value  of  \  1®  zero, 

it  corresponds  to  that  of  frozen  dissociation.  Case  of  inter¬ 
mediate  values  of  x  will  then  provide  qualitative  information 
on  the  non- equilibrium  flow  phenomena. 
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The  first  assumption  mentioned  above  represents  a  hypothesis  of 
using  an  idealised  chemical  system  in  place  of  the  actual  one,  while  the 
second  one  is  employed  mainly  for  mathematical  expedience.  It  is  easy  to 
see  that  the  value  of  x  should  vary  instead  of  remaining  constant  in  the 
flow  field. 

In  the  present  analysis,  these  two  restrictive  conditions  are  re¬ 
moved,  The  method  of  solution  follows  the  inverse  approach  previously 
explored  in  the  literature.  It  may  be  argued  that  since  the  predicted  re¬ 
action  rates  are  far  from  accurate,  the  first  assumption  used  in  Lick's 
analysis  is  thus  well  justified.  However,  in  the  continual  search  for  a 
better  solution  for  the  flow  field,  a  more  realistic  representation  of  the 
chemical  system  will  certainly  provide  additional  information  of  interest. 
Especially,  the  relative  difference  of  rates  among  oxygen,  nitrogen,  and 
nitric  oxide  may  have  significant  effects  on  the  final  solution  of  the 
flow  field. 

2.  Analysis 

For  a  prescribed  shock  shape,  the  conditions  immediately  behind 
the  shock  may  be  directly  determined  from  the  oblique  shock  relations 
by  assuming  that  immediately  behind  the  shock,  the  vibrational  degrees 
of  freedom  have  been  fully  excited  while  dissociation  has  not  yet 
begun.  Next,  a  smoothing  technique  such  as  a  polynomial  fit  may 
be  used  to  establish  the  partial  derivatives  of  these  conditions  with 
respect  to  a  direction  tangential  to  the  shock  surface.  Knowing 
the  values  of  the  variables  and  their  partial  derivatives  along  the  shock 
surface,  the  values  of  their  partial  derivatives  normal  to  the  shock 
surface  may  be  determined  from  the  governing  equations  of  motion  and 
the  equation  of  state  when  a  shock-oriented  curvalinear  coordinate 
system  is  chosen.  Linear  extrapolation  with  respect  to  the  direction 
normal  to  the  shock  surface  based  on  difference  formula  may  now  be 
used  to  calculate  the  values  of  the  variables  along  a  new  boundary  located 
a  short  distance  away  from  the  shock  surface.  The  same  smoothing 
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technique  establishes  the  partial  derivatives  along  this  new  boundary.  In 
turn,  the  partial  derivatives  normal  to  this  boundary  are  determined  from 
the  governing  flow  equations.  Successive  applications  of  this  scheme  al¬ 
lows  the  computation  to  be  extended  to  the  unknown  body  surface  which  will 
be  established  as  a  part  of  the  solutions  based  upon  the  consideration  of 
flow  continuity. 

The  equations  of  motion,  in  terms  of  the  shock-oriented  curvilin¬ 
ear  coordinates  may  be  written  (Fig.  24), 


8a.  8a.  cj. 

u  1  1  1 

"PrKyT5T  +  v-87  =  T 


W  (Pur)  +  ^  [pv  r  (1  +  Ky)]  =  0 


u|j  +  Kuv  +  (1  +  Ky)v§^  =  -  I|g 


By 


8y 


d(Eh,d.)  +  udu+vdv  =  0 


i  i 


The  equation  of  state  is  written  as 

p  =  (R PT 

Where  R  is  the  universal  gas  constant,  and  is  the  molecular  weight, 


(60) 

(61) 

(62) 

(63) 

(64) 

(65) 
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n  denotes  the  total  number  of  chemical  components  in  the  system.  The 
above  system  of  (n  plus  5)  equations  will  be  used  for  the  determination  of 
a^ ,  T,  o,  p,  u,  and  v,  which  completely  describe  the  flow  field. 
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In  Eq.  (60)  the  mass  rate  of  production  of  the  component  by 
chemical  reactions,  ,  is  given  by 


u 


i 


Mi  dvi 
"  X  dT 


(66) 


-a 


where  A  is  the  Avogadro'a  number,  v  is  the  number  of  particlee  of  the 

til  *  ^ 

iin  component  per  unit  volume,  (particle /cm J), 

In  general,  chemical  reactions  can  proceed  in  both  forward  and 
backward  directions.  A  set  of  opposing  chemical  reactions  may  be  repre¬ 
sented  by  (Ref.  31), 


II  X\£  H 

i.  +  )  a!  v.  <■ )  a!'  v,  +  u. 

j  L.  i  \  Z-  i  i 

i»l  b  i=l 


(67) 


where  aj  and  a!'  are  stoichiometric  constants. 

In  these  representations,  any  "third  body"  chemical  species,  p., 

J 

which  acta  as  an  energy  source  or  sink  in  a  chemical  reaction  but  is  other¬ 
wise  unchanged,  is  included. 

The  net  rate  of  production  of  ith  species  is  then  given  by 


dv. 

i 

dt 


(a!'  - 


ap 


n 

n 

i=l 


<vi> 


a! 

l 


V  |a 


i 


‘i)kb 


n 

n 

i  =  l 


(vp 


(68) 


At  thermodynamic  equilibrium  there  is  no  net  change  in  composition, 

therefore,  the  rate  constants  k.  and  k,  must  be  related  through  the  equi- 
•V  D  I 

librium  constant,  K.  That  is, 

^  n  a"  -  a ' 

(kj/Kh)  =  K  s  n  (v1>e)  .  (69) 
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When  multiple  chemical  reactions  are  involved,  the  net  rate  production  of 
ith  species  can  be  obtained  by  summing  Eq.  (68)  over  all  the  chemical  re¬ 
actions  in  the  system.  Where  Ifc  is  the  equilibrium  constant  in  the  unit  of 
(particle/cm^)  while  the  rate  constant,  k,  has  the  unit  of  [cm^/(sec  par- 
tide  )] .  In  the  case  of  atmospheric  air  the  dominant  chemical  reactions 
with  their  respective  rate  constants  are  listed  in  Table  1. 

The  rate  constants  for  all  the  chemical  reactions  with  various 

32 

"third  bodies"  used  in  this  analysis  are  those  given  by  AVGO.  It  should 
be  understood  that  these  numerical  values  of  the  reaction  constants  are  to 
be  considered  as  recommended  estimates  only.  Modifications  can  readily 
be  introduced  when  better  experimental  data  become  available.  Related 
discussion  on  chemical  kinetics  may  be  found  in  Ref.  33. 

Combining  EqB.  (60),  (66),  (68),  and  (69)  on  the  basis  of  chemical 
reactions  and  their  rate  constants  given  in  Table  1,  the  following  set  of 
five  equations  may  be  obtained  in  place  of  Eq.  (60): 
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where 


v  = 
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(75) 

In  obtaining  the  above  equations,  proper  care  has  been  taken  to 

3  i 

convert  the  rate  constants  kj  ,  k^  ,  in  the  unit  of  (cm  /particle  sec)  to  kjj  , 
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,  kj  in  the  units  given  in  Table  1,  Furthermore,  the  following 
non-dimensional  quantities  have  been  introduced: 


y  ■  * 


JL  5  =  -£- 


^i  - 


if  K 

r 


where  L,  p  ,  p  ,  T  ,  yft  ,  u  are  chosen  reference  quantities 
r  rr  r  (r  r 


u  =  *Jp  /  p 
r  rr  rr 


Pr  =  Pr  ^5")Tr 


In  terms  of  these  same  non-dimensional  quantities,  Eq.  (62)  may  be 


rewritten  as: 


-1  f  1  Sp  ,  v —  -  9u"l 

- - — —  K  u  v  +  u  ■— 

v(l  +  Ky)  p  8x  9xJ 


Eq.  (63)  may  be  rewritten  as: 
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Eq.  (61)  may  be  rewritten  as: 


=  __^£_f(l  +Ky)  9|+®|  +  |®| 
v(l  +  Ry)  9y  9x  p  8x 


+  (1  + 
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Differentiating  the  equation  of  state,  Eq.  (65),  with  respect  to  y,  we  have, 


where 


l  8T  _  1  Bp  1  8%  1  Bp 
T  By  p  By  7ft  8y  P  8y 


(87*/By>  =  -  SU/^jH 8a./8y) 


Dividing  the  partial  differential  form  of  the  energy  equation,  Eq.  (64),  by 
the  absolute  temperature,  we  obtain, 

isfh.  !!iU  *  0  (81 


T  By  i  \  1  dT/  T  i  \  1  By /  T  0y  T  By 
Substitution  of  (81)  into  (82)  yields, 

E  ia,  I  2  (jL  -  I  i£  +  i  S  (h.  — +  H.  ££  +  -  —  =  0  (83) 

i  1  dT/|_p  8y  i  ^  By  f  p  ByJ  T  i  \  1  By  /  T  By  T  By 

Eliminating  (0p/8y)  and  (Bp/8y)  from  Eq.  (83)  by  using  Eq.  (63)  and  Eq.  (61) 
respectively,  we  obtain 
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The  non-dimensional  form  of  Eq.  (84)  reads, 
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The  equation  of  state  may  be  rewritten  as, 


that  is, 


+  +  i22  +  Ifi  +  ; 

^NO  *0,  J 


PT 


(86) 


Now  the  reacting  flow  field  inside  the  shock  layer  is  governed  by  the  differ¬ 
ential  system  consisting  of  Eqs.  (70,  71,  72,  73,  74,  78,  79,  80,  85,  81, 
or  86).  There  are  ten  equations  for  the  ten  unknown  quantities 


P.  P  ■  T 
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In  this  inverse  method,  the  shape  of  the  shock  wave  is  prescribed. 

The  flow  conditions  immediately  behind  the  shock  wave  can  be  determined 
from  the  shock  wave  relations.  In  turn,  the  partial  derivatives  of  the  flow 
conditions  with  respect  to  x  can  be  numerically  calculated  because  of  the 
shock-oriented  coordinates  chosen  in  the  analysis.  From  the  governing 
differential  system,  we  then  calculate  successively  the  partial  derivatives 
of  the  unknown  quantities  with  respect  to  y  in  the  following  order: 

1)  Determine  {8a. /8y)  individually  from  Eqs.  (70)  through  (74). 

2)  Determine  (8u/By)  from  Eq.  (78). 

3)  Determine  (8v/8y)  from  Eq.  (85). 

4)  Determine  (9p/8y)  from  Eq.  (79). 

5)  Determine  (8p/8y)  from  Eq.  (80). 

6)  Determine  (8T/8y)  from  Eq.  (81). 

Knowing  the  partial  derivatives  normal  to  the  shock  surface,  linear  extra¬ 
polation  based  on  difference  formula  will  be  used  to  calculate  the  value  of 
variables  along  a  new  boundary  located  a  small  y  distance  away  from  the 
shock  surface.  Step  (6)  may  be  deleted,  and  the  temperature  distribution 
along  the  new  boundary  may  be  determined  directly  from  the  equation  of 
state,  Eq.  (86),  after  all  the  other  variables  along  tin  new  boundary  become 
known.  The  detail  of  the  computation  rnay  beBt  be  illustrated  by  the  fol¬ 
lowing  block  diagram. 


* 


75 


Where  0  denotes  any  of  the  dependent  variables  0m  is  the  value 
of  <p  at  the  mth  step.  The  values  of  (d<p/dx)Q  are  initially  obtained 

from  the  oblique  shock  calculations.  The  computations  are  carried  out  by 
applying  this  scheme  successively  until  the  body  streamline  1b  established 
from  continuity  consideration. 

Consider  a  streamline  crossing  the  curved  shock  at  an  arbitrary 
point.  B  (Fig.  25).  The  amount  of  fluid  enclosed  between  this  streamline 
AB  and  the  streamline  along  the  axis  of  symmetry  must  be  equal  to  the 
amount  that  crosses  the  section  BC  inside  the  shock  layer.  It  should  be 
understood  that  for  the  axisymmetric  flow,  the  control  surfaces  will  be 
obtained  by  revolving  AB  and  BC  about  the  axis  of  symmetry.  The  above 
mentioned  physical  condition  may  be  represented  by  the  following  equation: 


To  find  tiie  shape  of  the  body  surface,  we  numerically  integrate  the  right 
hand  side  of  the  equation  above- -a  process  which  can  be  carried  out  once 
the  distributions  of  p  and  u  are  determined  inside  the  shock  layer--and 
determine  the  yw  coordinate  at  which  this  equality  is  satisfied. 

After  the  body  streamline  shape  is  established,  the  mass  flow  rate 
is  then  integrated  from  the  body  surface  toward  the  shock  to  determine 
various  streamline  positions. 

Numerical  computations  (IBM7090)  for  selected  geometries  and 
flight  conditions  are  currently  being  determined. 


«f  * 
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SYMBOLS  (Section  II.  B) 


^o2,n 


N2  ° 


Avogadro  number 

specific  enthalpy 

h/hr  ,  where  hy  = 

reaction  rate  constant,  cm^/fsec  (particles)^] 
dissociation  rate  constant,  [cm^/(mole  sec)],  see  Table  I 

recombination  rate  constant,  [cm^/(mole^  sec)],  see  Table  I 

3 

reaction  rate  constant,  [cm  /(mole  sec)],  see  Table  I 

curvature 
KL  =  K  £2. 

3 

Equilibrium  constant,  particlcB/cm 

2 

Equilibrium  constant,  s  (Pq^Pq  ) 

2 

2 

Equilibrium  constant,  5  (p^  /p^  ) 

2 

Equilibrium  constant,  5  [  (Pj^o  Pq^Po  ,  PN^ 

2 

Equilibrium  constant,  =  [  (pN  p0)/(pN0  PN}1 

2 

reference  length,  =  fi. 

Mach  number 
molecular  weight 
pressure 


(p/pr)  ,  where  pr  *  pa 

radial  distance  from  the  axis  of  symmetry 

(r /a. ) 
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SYMBOLS  (Continued) 


Universal  gas  constant 

time,  sec 

temperature,  °K 

=  (T/T  )  ,  where  T  =  T 

r  r  oo 

velocity  component  along  the  x  direction 
=  u/ur  ,  where  ur  =  *'/pQO/pQO 

free  stream  velocity 

velocity  component  along  the  y-direction 

v/u  ,  where  u  =  sp  Jp 

r  r  rco 

coordinate  distance  along  the  shock  surface 
a  x/n  . 

coordinate  distance  normal  to  the  shock  surface 

=  y/tt . 

=  P-JP 

=  specific  heats  ratio  for  perfect  gaB 

3 

=  number  density  (particle/cm  ) 

3 

=  number  density  of  "third  body"  species  (particle/cm  ) 

=  density  ' 

=  shock  angle 
=  arbitrary  function 
=  radius  of  curvature 

=  radius  of  curvature  of  shock  wave  at  stagnation  streamline 


SYMBOLS  (Continued)- 

* 

Subscripts: 

i  ith  species  or  arbitrary  index  on  rate  constants 

j  jth  species  or  arbitrary  index  on  rate  constants 

0  quantities  atx  =  o  or  sea  level  standard  conditions 

r  reference 

s  conditions  immediately  behind  the  shock 

w  wall  conditions 


oo 


free  stream 


1  'ij  n  w  'WBlliWilff  mp-m 
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C.  Hypersonic  Wake  Flow  Studies 


1. 


Introduction 


As  a  body,  such  as  a  nose  cone,  enters  the  earth's  atmosphere, 

e 

it  leaves  a  deposit  of  energy  in  the  form  of  a  trail.  The  rate  at  which 
energy  is  deposited  is  equal  to  the  drag  energy  DU  where  D  is  the  aerody¬ 
namic  drag  and  U  is  the  flight  velocity.  Part  of  this  energy  goes  into 
heating  the  air  and  part  goes  into  giving  the  air  directed  motion  or  momen¬ 
tum  in  the  flight  direction.  It  is  the  heated  air  that  gives  rise  to  observ¬ 
ables  in  the  form  of  an  ionized  trail  and  gas  radiation. 

After  the  air  has  passed  through  the  bow  shock,  around  the  body, 
and  through  the  trailing  shock  system,  as  shown  in  Fig.  26,  it  returns  to 
ambient  pressure  at  some  distance  xq  behind  the  body.  At  thiB  station  a 
certain  velocity  and  temperature  profile  would  persist  downstream  if  it 
were  not  for  laminar  diffusion,  turbulent  mixing,  and  radiation. 

The  amount  of  energy  radiated  from  the  wake  is  extremely  small 
compared  iu  the  uthci  controlling  processes  of  energy  transfer,  and,  hence, 
can  be  neglected  in  any  wake  analysis. 

Figures  26a  and  26b  are  schematics  of  what  is  qualitatively  ex¬ 
pected  in  the  trail  behind  a  blunt-nosed  body  and  a  slender-pointed  body, 
respectively.  Both  bodies  are  assumed  to  have  about  the  same  base  area 
and  volume.  The  skin  friction  or  viscous  dragB  of  the  two  bodies  will  be 
roughly  the  same  at  the  same  velocity  and  altitude.  This  fact  is  reflected 
in  the  temperature  and  velocity  profiles  shown  in  Figures  26a  and  26b,  and 
is  represented  by  the  small  bumps  at  the  center  of  the  profiles.  This  is  a 
schematic  representation  only.  In  reality  the  momentum  defect  associated 
with  the  bump  should  include  not  only  skin  friction  drag,  but  nose  drag  in 
the  case  of  a  sharp  body  with  a  small  nose  radii,  and  the  drag  associated 
with  the  base  recompression  at  the  neck.  This  bump  in  the  velocity  profile 
constitutes  the  inner  core  and  its  spread  will  be  dominated  by  viscous  pro¬ 
cesses,  either  laminar  or  turbulent,  if  transition  has  occurred. 
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Figure  26b.  Schematic  of  Flow  Past  a  Sharp  Body 


:  W*  wmw-x  M*I 


0 

The  associated  bump  in  the  temperature  or  enthalpy  profile  will 
reflect,  in  addition,  the  heat  taken  out  of  the  flow  by  the  body,  and  its  peak 
will,  therefore,  be  much  lower  than  in  the  case  of  an  insulated  body.  As 
will  become  apparent  in  the  analysis,  heat  transfer  to  the  body  plays  a 
dominant  role  in  erasing  observables  in  the  wake  of  a  sharp  body  which  are 
dependent  on  the  maintenance  of  a  high-temperature  environment.  Although 
the  heat  transfer  to  the  blunt  body  is  of  the  same  order  as  the  heat  transfer 
to  the  sharp  body  at  the  same  velocity  and  altitude,  a  high-temperature 
environment  is  maintained  by  the  high  temperature  air  that  has  passed 
■  through  the  steep  portion  of  the  bow  shock  and  surrounds  the  viscous  core 
as  shown  in  Fig.  26a,  so  that  in  the  blunt-body  case  heat  transfer  to  the 
body  does  not  play  an  important  role  in  erasing  observables. 

2.  The  Spread  of  Turbulence  in  Axisymmetric  Wakes 

This  section  presents  a  simplified  analysis  of  the  spread  of  tur¬ 
bulence  in  the  trail  wherein  certain  simplifying  assumptions  are  made. 

These  assumptions  may  be  summarized: 

a.  Equilibrium  flow, 

b.  Negligible  radiation, 

c.  Negligible  laminar  diffusion  compared  to  turbulent  mixing, 
and 

d.  Turbulence  originates  in  the  viscous  core  or  on  the  body  and 

is  fully  developed  and  in  equilibrium  by  the  time  it  has  reached 
the  50  radii  station. 

Turbulence  is  thus  confined  to  a  circular  core,  the  front  of  which 
moves  out  and  "swallowB  up"  mass  from  the  outer  entropy  trail.  Eventually 
the  front  will  have  moved  out  into  the  outer  cool  air  at  the  edge  of  the  en¬ 
tropy  trail,  and  the  entire  trail  begins  to  cool  by  turbulent  mixing.  As  to 
how  long  this  process  takes  (or  how  far  downstream)  before  the  entropy 
trail  produced  by  the  normal  portion  uf  the  shock  of  a  blunt  body  is  cooled 
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sufficiently  to  erase  any  observables  each  as  electron  density,  will  depend 
upon  certain  initial  parameters  at  the  station  xq  .  These  parameters  are: 

a.  Total  drag  coefficient,  * 

b.  Ratio  of  drag  associated  with  the  turbulent  inner  core  to  the 
total  drag, 

c.  Width  of  the  entropy  trail,  and 

d.  Width  of  turbulent  inner  core. 

3.  Basic  Equations 

We  will  use  the  stream  function  4*  and  x  aa  independent  variables 
in  this  analysis.  The  stream  function  is  related  to  the  physical  radial  co¬ 
ordinate  r  through  the  formula 

r 

4*  =  2  tr  pur  dr  (88) 

"o 

ana  r  can  be  found  from  the  inverse  relation 


(89) 


Thus  the  radius  R  of  the  stream  tube  4^  in  the  free  stream  up  to  the  bow 
ehock  is 


R 


IT  O  U 
'  GO  00 


(90) 


Lot  4<{  be  the  stream  tube  which  enters  the  effective  front  of  the 

inner  turbulent  core  at  x  and  let  4*r  be  the  stream  tube  that  enters  at  x  , 

f  o 

the  origin  of  calculation.  Then  from  conservation  of  momentum  in  the 
x-direction  we  have 


(91) 
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since  there  is  no  x- momentum  transported  through  stream  tube  boundaries 
in  the  outer  inviscid  flow  and  the  pressure  is  assumed  to  be  consts&t. 

A  second  relation  is  found,  the  integral  relation  between  the  kinetic 
energy  and  the  dissipation  in  the  turbulent  core,  viz. 


d  f £  1  2  .  ,  1  2d^f  ,  2  f 

as  J  Iut  d*  =  ruf  -dx  -  4tt  J 


ft  2  2ut‘ 

*tr  ptut  ~r  d*  <92> 


By  substituting  the  profile 


U  +  (u,  -  u  )  -r- 
a  f  a  $7 


into  the  integral  relations  (91)  and  (92)  one  obtains 


I(uf  +  ua)4,£  =  I(uf  +ua  **f  +/  uid  + 

o  o  n  J 


e  (u* '  uaJ  -asr  ft 


jt'f 

=  4  ir  ZJ  Sr2Ptut(uf-ua)2^d  + 


Where  both  Eq.  (94)  and  its  differential  form  were  used  in  reducing  the  le£‘ 

hand  side  in  Eq.  (92)  to  that  in  Eq.  (95).  4^ ,  ua  >  and  uf  are  assumed  to  be 

purely  functions  nf  *  with  initial  values  of  ij'f  •  ua  >  and  u#  at  x  =  . 

Ao  o  40  ° 

Equation  (95)  can  then  be  further  reduced  to 


2  “f  •  "a  f f  2 

”  T 3-J  *«'  h-t 


In  order  to  obtain  a  simple  evaluation  of  the  dissipation  integral  on 
the  right  hand  side  of  Eq.  (96)  we  restrict  ourselves  to  the  case  where  the 


i^lf ■  tv  •-  . .  »tt-  <i8Ht^iiypi>EiqKW ■VFTifitwwnifii wi— n^t  r 


variations  of  velocity  and  density  are  small  across  the  core.  For  this  case 
we  make  the  "incompressible"  assumption 


T,  pt  Ut  Tr  pf  Uf 


and  approximate  «  by 


€  t  =  Kp.  r,  u.  ■  u 
t  f  f 1  i  a 1 


Substituting  these  into  (96)  we  obtain 


K,/+rpi  | 
nThV  uf  f  a' 


TVJ  =  12  KVSt;  luf  -  «alt«f  -  "a* 


How  from  the  relation  (94)  we  obtain  u  in  terms  of  u.  and  an 

a  f 


exorcssion  for  rl  u  /d:: 
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(uf  +  Ua  Hi 

(.•  DO 


2u,  duf 

— —  -  — —  -  (uf  +  u  )  ip 
+<•  d  k|>f  f  a 


+f 

f  -  A/  uicH 

°  J 


du  d  d>, 
a  f 

cT^  tHT 


Since  m  is  assumed  to  be  known  as  a  function  of  41 ,  it  is  a  simple  matter 
to  obtain  x  as  a  function  of  ^  by  numerical  integration  of  Eq.  (99). 

We  have 


r  £  d  „  /S',  u, 

- — - L -  \  J—L  di|n  =  12KVtt  (x  -  j:  ) 

J  Tuf  -  -al V  Pf  f 
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Prom  the  graph  of  x  vi  ^  we  can  obtain  4*£  vs  x  and,  hence,  the  flow 
quantities  and  equilibrium  electron  densities  as  a  function  of  x  and  t|i  and 
ultimately  as  a  function  of  x  and  r  through  the  relation  (89). 

This  formulation  is  for  an  arbitrary  as  a  function  of  +  .  We 
shall  look  at  the  special  case  where  u.  can  be  approximated  by 


u.  (*)  =  (1  -  ae’b+) 


(103) 


This  distribution,  by  the  proper  choice  of  a  and  b,  can  be  made 
to  fit  the  velocity  defect  produced  by  a  blunt  body  fairly  well.  It  gives  the 
bell-shaped  nearly  Gaussian  curve  for  the  outer  entropy  trail  shown  in  Fig.  26a 
when  plotted  in  the  physical  coordinate  r  .  (Note  that  ^  from  Eq.  (88.) 

At  the  initial  station  x  =  x  let  the  drag  associated  with  the  inner 

o 

turbulent  core  be  D  and  let  D  be  the  total  drag.  We  then  have 
w 


rf 

I  '  1 


\  I  » 

/  u  y 


I  (uf  *  ua  > 
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(104) 


OJ 

D  -  DW  =  (Uoo  •  uf  }  +f  +  /  (Uoo  '  ui}  dl|j  =  <Uco  •  uf  >  ^ 
n  i )  r%  r 


+  ~V  e'b+f  =  i  +  (U  -  uf  ) 
b  oo  o  [b  \  l  co  ij 


(105) 


from  which  one  readily  obtains  the  relation 


ri  u,  -  u 
0  fo  ao 


(106) 
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We  sh&ll  introduce  at  this  point  the  dimensionless  quantities: 


* 

u 

= 

u/u 

CD 

(107) 

p 

= 

p/pm 

OO 

(108) 

T 

= 

T/T^  P^/p  =  1/p  (since  p  =  const.  =  p^) 

(109) 

CD 

= 

D4  "«  Uoo2  A 

(110) 

CD  = 
\v 

°w4',a»',»  A 

dll) 

k 

= 

CD  /CD 

(112) 

w 

= 

^b 

(113) 

+ 

= 

l/h\ 

(114) 

>>  — # 

n 

= 

b^f  =  /y 

(115) 

c  * 

CD 

= 

V* 

(116) 

Uco 

flight  speed 

= 

atmospheric  density  at  flight  altitude 

T 

oo 

= 

absolute  temperature  at  flight  altitude 

Pm 

= 

ambient  pressure  at  flight  altitude 

A 

= 

re-entry  body  maximum  x- sectional  area  =  it  r^2 

rb 

= 

maximum  body  radius 

ik 

0 

= 

p  U  A  -  stream  function  value  of  free- stream 

00  CO 

tube  equal  to 

body  x-eectional  area 

D 

= 

total  drag 

D  =  drag  contained  in  turbulent  core  at  the  start  of  calculation 
Q  =  total  heat  transfer  rate  to  the  body 

Q*  -  Q/i  ^  =  dimensionless  heat  transfer  rate  coefficient. 
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Equation  (101)  can  than  be  written  in  non-dimensional  form,  viz 


!w »-r  •  m  •f5www*"t«aw 


ir; 


It  can  be  seen  immediately  that,  If  Eqs.  (120)  and  (121)  are 

substituted  into  Eq.  (117),  a  factor  l/C—  can  be  pulled  out.  Thus  multi- 

$  U 

plying  Eq.  (117)  by  we  obtain 


Y  •  k  ,  n 


;  n)dn 


12KC*  X  -  X 
D  o 


=  g 


(125) 


with 


1-k  'o 

TTTe 
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.-n 


F  = 


'••nr' 


1  - 


1  -k 

1+  T) 


e  °  (1  +ti)  e 


-*1 


I  * 

VI  uf  Tf 


(126) 


4.  Discussion  and  Numerical  Examples 

Examination  of  Eq.  (125)  and  (126)  reveals  that  the  major  param¬ 
eters  that  control  the  spreading  of  the  turbulent  core  occur  only  in  the 
combinations 


1  -  k 
1  +,1o 


1_!m2 

2  oo 


1 2  K  C 


D 


kC 


D 


«/j* 


=  1 


-  u 


and  that  (x  -  x^/r^  can  be  expressed  entirely  in  terms  of  these  parameters. 
It  can  further  be  shown  that  these  parameters  can  all  be  expressed  in  terms 
of 


$  $ 

v  ,  M  ,  K  ,  ,  k  ,  T,  ,  u, 

OO  U  i  I 

o  o 


Although  the  exponential  variation  of  velocity  given  by  Eq.  (103) 
was  chosen  on  the  basis  of  a  good  fit  for  a  blunt  body  with  a  spherical  nose, 
the  approximation  can  still  be  used  for  sharp  bodies,  such  as  a  cone,  with¬ 
out  loosing  any  of  the  salient  features  in  describing  the  turbulent  core 
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spreading  and  temperature  on  the  axis  as  a  function  of  (x  -  xQ)/r^  .  This 
is  due  to  the  fact  that,  for  a  sharp  body,  the  radial  distribution,  of  velocity 
and  temperature  in  the  inviscid  flow  is  fairly  flat  and  the  gradients  at  the 
edges  are  not  important.  This  is  demonstrated  in  the  two  numerical  cal¬ 
culations  discussed  in  the  ensuing  paragraphs. 

Two  numerical  calculations  at  150,  000  ft  and  Mach  20  are  pre¬ 
sented:  one  for  a  typical  blunt  body  such  as  that  shown  in  Fig.  26a  and  one 
for  a  typical  sharp  body  shown  In  Fig.  26b.  The  drag  contained  in  the  tur¬ 
bulent  core  bucket  is  assumed  to  be  the  same,  that  is,  =  .06.  This, 

uw 

in  both  cases,  includes  the  skin  friction  and  base  drag,  plus  nose  drag  in 

the  case  of  the  sharp  body.  At  150,  000  ft,  the  skin  friction  coefficient  is 

about  0,03,  and  using  Reynolds  analogy,  the  heal  transfer  coefficient  is 
*  13. 

OAr  !■»  '>  A  V.Oii  .  These  values  were  assumed  lor  both  cases.  Both 

i.  '  UL.I  OU 

calculations  were  started  with  the  assumption  that,  the  velocity  on  the  axis 
was  zero  at  the  beginning  which  would  be  the  case  at  the  effective  rear 
stagnation  point,  Each  case  will  now  be  discussed  separately. 

a.  Blunt  Body 

A  typical  pressure  drag  coefficient  for  a  blunt  body  is  0.54, 
giving  a  total  drag  coefficient  of  Cn  =  0.6  with  k  =  0.1  .  The  results  are 

u  *■’  >’f 

shown  in  Fig.  27  where  *| ,  T,  and  T  with  and  without  heat  transfer  are 
plotted  against  (x  -  xQ)/r^  .  A  measure  of  how  far  thefront  has  moved 
laterally  in  terms  of  Caussion  widths  of  the  entropy  trail  is  given  by  r| , 

It  can  be  seen  that  in  the  blunt  body  case,  the  front  must  move  many 
Gaussian  widths  before  the  center  of  the  turbulent  core  begins  to  cool. 

This  distance  is  of  the  order  of  a  hundred  body  radii.  The  solid  curve  for 
T&  is  the  temperature  ratio  on  the  axis  for  an  insulated  hody,  while  the 
dashed  curve  is  for  the  assumed  value  of  the  heat  transfer  rate.  It  can  be 
n-aen  that  the  heat  transfer  rate  to  the  body  has  little  effect  in  depressing 

$  Jjt 

the  temperature  except  at  the  very  beginning  where,  for  u  =  0,  T  is 

Si  3. 

about  0.56  times  the  free  stream  stagnation  temperature  (strictly  speaking 
$ 

T  is  the  enthalpy  ratio  and  not  the  temperature  ratio). 

Hf* 
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It  can  be  concluded  that  for  a  typical  blunt  body  with  a  two-foot 
diameter  the  temperatures  will  be  sufficiently  high  to  produce  large  amounts 
of  radiation  and  sustain  high  electron  concentrations  visable  to  radar.  Fur¬ 
thermore,  the  width  of  the  turbulent  core  becomes  very  large  for  the  blunt 

body  case,  as  evidenced  by  the  large  values  of  *1  in  terms  of  body  radii. 

$  $ 

When  T  has  reached  roughly  half  the  value  of  T,  ,  the  initial  front  tem- 
a  ta 

perature  in  the  entrophy  trail,  q  has  the  value  5  corresponding  to 

rf^rb  =  aa  ca^cu^atec^  irom  Eq.  (89),  and  (x  -  xQ)/r^  is  117, 

b.  Sharp  Body 

A  typical  pressure  drag  coefficient  for  a  sharp  body,  say  a  10° 
half  angle  cone  or  an  equivalent  ogive,  ii  0.0b  giving  a  total  drag  coefficient 
of  Cp  =  0.12  with  k  =  0.5  .  The  results  are  shown  in  Fig.  28  where  the 
same  quantities  are  plotted  as  in  the  blunt  body  case.  It  can  be  seen  that 
the  wake  is  cooled  in  a  much  shorter  length  than  in  the  blunt  body  case. 

>}c 

The  temperature  on  the  axis  has  reached  the  value  T&  a  9  at  a  value  of 
(x  -  XQ)/rh  =  33,  with  n  a  2.1  for  a  value  of  rf/rb  =  5.4. 

It  can  be  concluded  that  for  a  typical  sharp  body  the  dimensions 
of  the  high  temperature  region  in  the  wake  are  about  a  factor  of  three  less 
than  for  the  blunt  body  case.  The  heat  transfer  rate  to  the  body  has  a  more 
significant  effect  on  the  reduction  of  these  dimensions  than  in  the  blunt  body 
case. 

The  values  of  the  parameters  chosen  for  the  sharp  body  are 
conservative  and  they  tend  to  make  the  wake  larger  than  that  which  might 
actually  be  expected  at  this  altitude. 


IS.  EXPERIMENTAL  INVESTIGATIONS 


A.  Stagnation  Point  Heat  Transfer  Measurements  in 
Hypersonic  Low  Reynolds  Number  Flows 

1.  Introduction 

In  determining  the  point  at  which  protective  fairings  for  delicate 
payloads  should  be  jettisoned  during  the  launch  trajectory  of  satellites  ve¬ 
hicles,  free  molecule  flow  heat  transfer  rates  have  been  assumed.  Calcu¬ 
lations  based  on  this  assumption  provide  the  maximum  heat  transfer  rate 
predictions,  and  hence,  the  maximum  payload  protective  fairing  require¬ 
ments.  As  a  result,  the  tendency  has  been  to  carry  the  fairings  to  altitudes 
which  are  much  beyond  optimum  because  of  the  uncertainty  of  the  heat 
transfer  rates  in  the  very  high  atmosphere.  These  rates  were  not  signifi¬ 
cant  relative  to  exit  and  re-entry  heat  rates  and  these  unknowns  were 

tolerable.  Now  situations  are  arising  where  it  is  no  longer  advisable  to 

34 

ignore  these  uncertainties.  Also  theoretical  analysis  in  the  low  density 
regime  is,  for  the  most  part,  unsubstantiated  by  experiment.  It  is,  there¬ 
fore,  the  purpose  of  this  study  to  attempt  to  measure  the  stagnation  point 
heat  transfer  rates  to  cylinders  and  spheres  in  the  hypersonic  shock  tunnel 
at  flow  Mach  numbers  from  13  to  23.  The  Reynolds  number  range  to  be 
covered  is  from  5  to  500,  based  on  leading  edge  diameter,  corresponding 
to  the  regime  between  free  molecule  flow  and  the  continuum  regime. 

Because  of  its  high  reservoir' stagnation  temperature,  the  shock 
tunnel  is  an  instrument  peculiarly  suited  to  making  measurements  at  high 
Mach  numbers  and  low  Reynolds  numbers  typical  of  the  high  altitude  por¬ 
tions  of  the  launch  trajectory  of  satellite  vehicles.  Furthermore,  the  shock 
tunnel  can  yield  a  variety  of  stagnation  pressures,  and  hence,  Reynolds 
numbers  for  a  given  stagnation  temperature.  The  chief  drawback  inherent 
In  shock  tunnels  is  the  comparatively  short  time  interval  during  which 
measurements  can  be  made.  As  a  result,  equilibrium  conditions  between 
the  flow  and  the  measuring  instrument  usually  cannot  be  obtained  and  must 
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be  infer. red  from  the  transient  response  of  the  measuring  device.  In  thiB 
report  heat  transfer  rate  measurements  determined  from  th?  transient 
response  of  thin  film  resistance  thermometers  mounted  at  the  stagnation 
point  of  models  will  be  described  and  the  results  of  the  measurements  pre-; 
sented,  The  shock  tunnel  facility  in  which  these  tests  were  made  is  shown 
in  Fig.  29. 

2,  Theoretical  Predictions 

_  til 

Jt-'robstein  and  Kemp'  have  considered  the  problem  of  calculating 
the  aerodynamic  heating  to  be  expected  in  hypersonic,  rarified  gas  flows, 
where  the  mean  free  path  becomes  loo  large  for  the  use  of  the  inviscid 
f  low-boundary  lay  or  approach  and  too  small  for  the  tree  molecule  flow 
regime  to  apply.  They  have  solved  the  Navier  stokes  equations  for  two  seta 
of  boundary  conditions;  (1)  a  "viscous  layer”  class  in  which  a  thin  shock 
wave  exiBts,  but  the  region  between  the  shock  and  the  blunt  body  is  fully 
viscous,  and  (2)  a  "merged  layer"  class  representing  a  more  rarified  con¬ 
dition.  In  the  merged  layer  the  shock  wave  is  no  longer  thin  and  the  Navier 
stokes  equations  are  used  to  provide  a  solution  which  includes  the  shock 
structure  and  free  stream  conditions  as  outet  boundary  conditions 

From  this  analysis  a  deviation  is  predicted  in  the  stagnation  point 
heat  transfer  rate  from  that  expected  on  the  basis  of  extrapolated  boundary 
layer  theory.  The  viscous  layer  theory  predicts  higher  heat  transfer  to 
both  spheres  and  cylinders  normal  to  the  flow  as  the  Reynolds  number 
(based  on  nose  diameter)  is  decreased  to  the  order  of  10. 

The  "incipient  merged  layer"  theory  modifies  the  heat  transfer 
rates  to  values  which  are  less  than  in  the  viscous  layer,  but  greater  than 
in  the  boundary  layer  theory.  The  predicted  results  for  spheres  and  cyl¬ 
inders  are  shown  in  Fig.  30a  and  30b.  Although  an  incipient  merged  layer 
prediction  for  the  cylinder  is  not  made,  it  is  easy  to  suppose  that  it  lies 
between  the  viscous  layer  and  boundary  layer  value.  These  curves  were 
computed  on  the  basis  of  a  density  ratio  across  the  shock  of  0.1  and  a  y 
for  the  gas  of  11/9. 
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igure  30a.  Heat  Transfer  Rate  to  Cylinders  Normal  to  the  Flow  as 
Function  of  Stagnation  Point  Reynolds  Number 


For  shock  tunnel  reservoir  conditions  of  about  240  atmospheres 
and  2640°K,  and  test  section  Mach  numbers  from  10  to  20,  the  heat  trans¬ 
fer  rates,  according  to  Fay  and  Riddell,^  were  calculated.  These  values 
are  plotted  in  Fig,  30  for  comparison  with  the  results  of  Probstein  and 
Kemp.  A  free  molecule  flow  value  for  an  accommodation  coefficient  of 
unity  is  plotted  also. 

3.  Experiments 

In  the  shock  tunnel  for  shock  Mach  numbers  of  4.5  to  4.7,  reser¬ 
voir  test  conditions  (which  remain  steady  for  about  4  millisec/  r  V-*. 

220  and  240  atm  and  2S40  and  3000  K  are  obtained.  Using  the  twelve-inch 
diameter  test  section  and  varying  the  throat  diarn.  ?.<•/  u»cr  a  range  of  .04 
to  .375  inches,  Mach  numbers  from  10  to  23  i.re  achieved.  For  each  test, 

the  shock  Mach  number  (M  )  in  the  shock  tube,  and  the  pitot  pressure  (P  ) 

8  t 

in  the  test  section,  are  measured.  From  the  shock  Mach  number,  the 
reservoir  conditions  (that  is,  the  conditions  behind  the  reflected  shock) 
are  ascertained.  These  conditions,  along  with  the  pitot  pressure  measure¬ 
ment,  determine  the  flow  Mach  number  in  the  test  section 

Spherical  and  cylindrical  shells,  of  Pyrex  7710  glass  ranging  in 
diameter  from  3  to  20  mm,  were  employed  as  models  Two  thin  film  re¬ 
sistance  thermometers  were  applied  at  the  stagnation  point,  of  the  cylinder 
and  half  an  inch  on  either  side  of  the  mid-point  of  the  cylinder.  The  films 
were  applied  by  airbrushing  Hanovia  No.  05X  liquid  bright  platinum  paint 
in  a  strip  0.5  mm  wide  and  half  an  inch  long  and  th  in  firing  it  at  i250°F. 
Firing  silver  micro  paint  in  a  heavy  layer  provides  low  resistance  contacts 
with  the  measuring  circuit.  Fig.  31  is  a  photograph  of  the  cylindrical 
model.  This  model  was  mounted  in  the  test  section  normal  to  the  flow  so 
that  the  rrid-point  of  the  cylinder  corresponded  with  the  axis  of  the  tunnel 
teBt  section.  The  cylinder  was  supported  from  the  wall  in  cantilever  fash¬ 
ion.  The  spherical  model  shown  in  Fig.  32  has  a  shorter  thin  film  at  its 
stagnation  point  and  is  supported  from  the  model  stand  as  illustrated.  The 
pitot  pressure  gauge  location  is  indicated  also  in  this  photograph 
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‘igure  32.  Mounting  of  the  Spherical  Model  in  the  Tunne.  Showing  also 
the  Pitot  Tube  for  Measuring  Test  Section  Total  Pressure 


The  technique  of  thin  film  resistance  thermometry  is  quite  well 
known  and  very  adequately  covered  in  the  literature.  In  essence,  the  ther¬ 
mometer  measures  the  surface  temperature  history  of  the  material  on 
which  it  resides.  This  measured  temperature  history  is  converted  to  a 
heat  transfer  rate  by  the  use  of  one- dimensional  heat  conduction  theory  and 
one  basic  calibration  of  the  thermometer.  The  basic  calibration  is  gener¬ 
ally  done  by  the  pulsing  technique  in  which  a  known  heat  flux,  produced  by 
a  current  step  pulse  through  the  resistance  thermometer,  is  applied  to  the 
gauge.  The  calibration  is  used  to  determine  the  lumped  thermal 
of  the  substrate  of  the  thin  fiim,  \'tr  pCk,  where  p  ,  C  ,  and  K  are  the  den¬ 
sity,  specific  heat,  and  conductivity  of  the  substrate.  The  thermal  coeffi¬ 
cient  of  resistance  of  the  film  is  determined  by  measuring  the  resistance 
at  the  freezing  and  boiling  points  of  water. 


The  equation  relating  the  heat  transfer  rate  (Q)  to  the  model  along 
with  the  temperature  history  of  the  model  surface  (ATM)  is  given  as 
r  'I 

AT 


(127) 


T  is  the  temperature  change  of  the  surfai e  and  I  is  the  initial  tempera- 

r 

ture  of  the  surface.  The  quantity  in  the  braces  J  represents  the  correc¬ 
tion  to  the  heal  transfer  rate  which  must  be  included  in  order  to  account  for 
the  change  in  the  thermal  properties  oi  the  substrate  when  the  surface 
temperature  becomes  large  in  comparison  with  the  calibration  tempera¬ 
ture.35  b/k  is  a  function  of  the  material,  which  for  Pyrex  is  4.73 
o 

Measurements  were  made  in  the  shock  tunnel  of  the  heat  transfer 
rates  to  cylinders  3mm.  10mm,  and  20mm  in  diameter  and  to  spheres  10 
and  20mm  in  diameter.  Reservoir  stagnation  conditions  ranged  from  2800 
to  3000°K  with  pressures  of  28  to  29  atm  and  220  to  240  atm.  Mach  num¬ 
bers  in  the  test  section  varied  from  10  to  23,  but  most  of  the  data  was  in 
the  range  of  13  to  16. 


4.  Results 


In  calculating  the  value  of  Ret  ,  the  Reynolds  number  based  on 

leading  edge  diameter,  the  value  of  the  viscosity  is  determined  according 

to  the  results  of  the  study  of  oxygen  recently  reported  by  Hartunian  and 
36 

Marrone.  This  study  indicates  that  the  Sutherland  viscosity  values  based 

on  kinetic  theory  and  hard  spherical  molecules  are  as  much  as  50%  too  low 

for  temperatures  up  to  3000°K.  This  modified  value  of  u  also  was  used 

1  Z 

to  determine  the  theoretical  value  of  ,  according  to  Fay  and  Riddell, 

Mach  numbers  and  the  heat  transfer  rates  are  tabulated  in  Table  I.  Values 
of  Q^Re^/ipU*  as  a  function  of  Ret  are  plotted  in  Fig.  30a  and  30b,  along 
with  the  theoretical  curves.  The  results  for  the  cylinder  seem  to  follow 
the  boundary  layer  theory  right  on  down  to  the  free  molecule  flow  regime 
with  only  a  slight  increase  at  Reynolds  numbers  less  than  20.  The  consid¬ 
erable  scatter  in  the  data  makes  a  definitive  determination  of  the  best 
theory  impossible. 

The  experimental  data  for  the  axisymmetric  probe  (sphere)  are 
even  more  scattered,  and  therefore,  inconclusive. 

Some  remarks  on  the  accuracy  of  these  expci  iments  should  be 
made.  The  shock  Mach  number  can  bo  measured  to  ±  1%  and  this  is  about 
the  error  expected  in  the  reservoir  conditions  The  total  pressure  can  be 
measured  to  within  ±  10%  and  this  should  be  the  expected  error  in  p  .  As 
a  result,  the  theoretical  curve  which  is  plotted  from  the  experimentally 
obtained  valueB  for  the  tunnel  conditions  could  have  a  width  of  ±  10%.  The 
heat  transfer  measurements  have  a  possible  error  of  between  ±  5%  and 
±  18%,  depending  on  the  smoothness  and  clarity  of  the  heat  transfer  record. 
An  example  of  a  typically  adequate  record  is  shown  in  Fig.  33  for  a  flow 
Mach  number  of  16,5  and  reservoir  conditions  of  243  atm  and  3100°K. 

No  explanation  is  immediately  available  for  the  greater  amount  of 
scatter  in  the  sphere  data.  Since  the  spheres  were  blown  from  Pyrex  tub¬ 
ing  some  irregularities  existed  in  their  radii.  However,  the  heat  transfer 


Figure  33.  Typical  Temperature  History  for  the  Stagnation  Point  of  a 

Cylindrical  Model 

P  =  222  atm 

T  =  2980°K 

M  =16.0 
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rate  varies  inversely  as  the  square  root  of  the  body  radius,  and  hence, 
irregularities  in  the  radii  should  not  give  significant  effects.  The  much 
smaller  gauge  on  the  spheres  could  limit  the  accuracy  of  the  results. 

5.  Conclusions 


The  heat  transfer  rates  to  cylinders  and  spheres,  as  measured  in 

p 

the  shock  tunnel,  are  adequately  predicted  by  the  theory  of  Fay  and  Riddell 
down  to  Re  £  20  .  At  this  point,  the  heat  transfer  rate  increases  some- 
what,  as  predicted  by  Probstein  and  Kemp,  ’  reachinp  thp  fro*  mnUHiU 
licw  Jieat  transter  rate  at  a  value  nf  between  2  and  3, 


107 


B .  Electron  Beam  Densitometer  Studies 
1 .  General 

When  an  electron  beam  is  projected  through  a  space  containing 
molecules  of  a  gas,  the  electrons  are  scattered  due  to  both  elastic  and 
inelastic  collisions  with  gas  molecules.  The  incident  electron  loses 
energy  and  causes  the  parent  beam  to  become  diffuse  and  inhomogeneous. 
Electron  scattering  in  a  gas  is  a  function  of  electron  energy,  particle  den- 
pitv  molecular  croce  enrH-iwr  •»*--«*  1  -  !  L  T'  *  ’  *" 

to  measure  tmj  molecular  number  density  tor  a  known  gas  by  electron 
scattering. 

A  high  energy  beam  of  electrons  was  utilized  to  measure  gas 
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density  in  the  flow  field  of  the  low  density  wind  tunnel  at  Berkeley.  This 
is  a  continuous -run  tunnel  operating  at  Mach  2  with  free  stream  pressures 
of  about  0.1  millimeters  mercury,  The  Berkeley  electron  beam  apparatus 
was  entirely  contained  within  the  large  diameter  loot  sei  ilun  of  the  tunnel. 


Perturbation  ol  the  flow  and  limitations  of  response  time  and 
sensitivity  preclude  the  use  of  more  conventional  pressure  sensors  in  a 
hypersonic  low  density  shock  tunnel.  The  electron  beam  probe  is  thus  a 
theoretically  ideal  sensor  sinc  e  cle'  trons  are  effectively  massless  c  om¬ 
pared  to  the  free  stream  flow  and  have  a  high  order  of  sensitivity  and 
response  speed. 


2 .  Electron  Beam  Att e nuation  in  a  Gas 

When  an  electron  beam  is  projected  into  a  gas,  a  fraction  of  the 
initial  beam  will  survive  the  passage  without  deflection  and  its  path  is  a 
continuation  of  the  original  ballistic  trajectory.  If  only  the  unscattered 
electrons  arc  accepted  by  a  suitable  means  of  detection-  Lambert's  law  of 
absorption  will  hold  for  the  loss  of  beam  intensity  along  its  path:  * 


dl 

dx 


U1  T 

=  V-  P1 


(128) 


I 
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Where  I  is  the  beam  intensity,  x  is  the  path  length,  p  is  the  gas  density, 
and  p  is  an  absorption  coefficient. 

The  absorption  coefficient,  p,  is  related  to  the  effective  scatter¬ 
ing  cross  section  of  the  gas.  The  cross  section  Q  ($)  is  defined  as  the 
probability  per  scattering  center  that  an  electron  will  be  deflected  by  an 
interaction  with  a  gas  molecule  into  a  cone  of  half  angle  which  is  greater 
than  the  effective  angle  subtended  by  the  detector  aperture  measured  from 
the  point  of  scattering,  If  there  are  n  molecules  oi  aas  in  the  vr.himp 
■„-_v  upiod  by  the  electron  beam,  i'ne  totai  probability  of  an  electron  inter¬ 
action  with  a  gas  molecule  is  nQ  ,  or  over  an  increment  of  path  length  dx, 
is  nQdx,  and  the  loss  of  beam  intensity  becomes: 


dl_ 

dx 


nQT 


(129) 


Integration  of  Eq.  (129)  yields  I  (x)  =  Iq  E‘nQx  (130) 
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Markevitch  and  Huribut  hnvp  shown  that  the  experimental  vaiue 
of  Q  (rr)  for  high  energy  electrons  in  air  is  close  to  the  theoretical  value 
based  on  clastic  collisions.  This  value  is  energy  dependent  and  has  been 
found  to  have  the  relation 


Q  iir) 


3.6  x  10'13 
V 


2 

cm 


(111) 


where  V  is  in  volts. 

Taking  Avogadro's  number,  N  -  6.025  x  10^3/mol.,  and  the 
standard  volume,  V,  =  22.4146  x  103  cm3/mol.,  then  at  a  pressure  of  one 

°  13 

micron  there  are  3.537  x  10  molecules  per  cubic  centimeter  of  a  perfect 
gas.  Then  n  =  3.537  x  1013  p,  and  pp  =  nQ  =  3.537  x  1013  pQ,  from 
which 

Q  =  2.827  x  10"14  (132) 

P  v 


U0 


I 


t 


i 
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This  quantity  Q  is  the  total  collision  cross  section  of  the  gas,  and  has  the 
same  significance  as  the  sense  of  Q  (it)  in  Eq.  (131);  that  is,  it  represents 
the  total  cross  section  which  was  determined  from  scattering  into  a  cone 
or  half  angle  ir  radians. 

We  mav  now  obtain  the  absorption  coefficient  p  by  combining 
Eq.  (131)  and  (132): 

HP  =  12.73  £  (133) 

ior  air,  where  p  is  in  microns  and  V  is  in  volts. 

The  equation  for  the  intensity  of  an  electron  beam  in  air  at  a 
constant  pressure  p  microns  is  then 


1  (x)  =  Iq  E' 12-73  2Z  (134) 

where  I  (x)  is  measured  in  centimeters  from  the  point  of  origin  ana  Iq  is 
the  incident  beam  current  of  an  energy  V  electron-volts. 

Equation  (134)  is  plotted  on  logarithmic  coordinates  in  Figs.  34 

and  35. 

3,  Electron  Beam  Absorption  in  an  Axially  Symmetric  Flow  Field 

The  electron  beam  densitometer  will  be  used  to  obtain  free  stream 
density  measurements  in  ihe  flow  field  of  the  hypersonic  shock  tunnel  and 
the  variation  of  free  stream  density  in  the  direction  normal  to  the  flow. 
Consider  a  uniform  beam  of  electrons  injected  into  an  axially  symmetric 
flow  field  where  the  gas  density  is  a  function  of  r  ,  the  radial  distance 
from  the  central  flow  axis,  as  in  Fig.  36.  It  will  be  required  to  find  the 
gas  density,  p  ,  as  a  function  of  r  .  Consider  a  uniform  beam  of  electrons 
ol'  intensity  Iq  originating  outside  the  flow  region  and  directed  in  the  posi¬ 
tive  x  direction. 
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Lambert's  law  of  absorption  holds  as  a  partial  derivative  Bince 
the  beam  attenuation  is  a  function  of  x  and  y  . 


which  has  a  solution 


(135) 


1  n  T/tr  \  ..  In  T 
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where  p  «  f (r )  and  is  a  constant. 

2  2  2 

From  the  geometry  of  Fig.  36,  r  =  x  +  y  ,  and  for  constant  y  , 
rdr  =  xdx  =  <Jr^  .  dx-  Also,  when  x  =  o,  r  =  y,  and  when  x  =  x, 
r  -  R  then  Eq.  J136)  becomes 

R 

F(y)  =  -2p  f  p(r)  ■■  (L  -  2x)  (137) 

y  \ ]  r“  -  y** 

2  2 

A  second  change  of  variables  is  made  by  putting  S  =  R  -  r  and 
g(s)  =  p(r).  Now  dS  =  -2rdr,  and  at  r  =  y ,  s  -  R2  -  y2  and  at  r  =  R 
S  =  0.  Then  Eq.  (137)  becomes 


F(y)  =  (i 


l  , 

R  -y 


g(s) 


ds 


J~RZ  -  y2  -  S 


-  ppj  (L  -  2x) 


(138) 


2  2 

This  may  be  simplified  by  substitution  of  t  =  R  -  y  ,  then 


c 

F(y)  =  -  , -  V-Py  (L  -  2*) 


(139) 


\ 
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and  define 


f(y)  =  G(t)  =  -  (j. 


f  g(s)ds 

J  JT- 


(140) 


f(y)  =  F(y)  +  jlPi  (L  -  2x) 


(141) 


Equation  (140)  in  cu  me  A&ei  type  with  unknown  g(s)  and  has  the 


solution 


g(8)  =  ll  -iL  f  GjOdt, 


(142) 


To  integrate  by  parts,  choose  U  =  G  (t)  and 


dV  = 


Equation  (142)  becomes 


g(s)  =  +  _L  f  o’  (t>dt 

VT  ^  J  ^/sT7 


From  Eq.  (140)  G(o)  =  o  and 


G'(t)  =  f’(y)  *  --  f(y) 
dy 


so  that 


5 

JL  [  f'(y)dt 
T T\i  J  S  -  t 


(144) 
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2  2 

We  now  change  variables  by  the  substitutions  t  =  R  -  y  and 
_  2  2 

S  =  R  -  r  .  Note  that  when  y  =  R  ,  t  =  o  and  when  y  =  r  ,  t  =  a  ,  then 


g(s)  =  p  (r) 


-  -  f~ 

*i*  j  r 


(y)dy 


J  JIT  2 

r  Vy  -  r 


(145) 


By  differentiation  of  Eq.  (141), 


f'  (y)  =  ~  lnl(y)  -  2  pp 


f~l  2 

vR  -  y 


(146) 


and  substitution  into  Eq.  145, 


p(r) 


■ih . 'A?  / 


ydy 


J  f  i  2  JZz  Z 

r  V  y  -  r  VR  -  y 


(147) 


p(rl  =  i?  /£1"My)Xi7*  +  p' 

r  y 


(148) 


I (y )  may  be  obtained  by  measuring  segments  of  the  core  flow 
under  identical  tunnel  conditions  and  p(r)  is  evaluated  numerically.  A 
preferred  technique  is  to  incorporate  a  series  of  detectors  and  a  rapid 
sweep  on  the  electron  beam  to  measure  all  segments  during  a  single  run. 

4.  Instrumentation 

The  electron  beam  densitometer  consists  of  two  major  systems: 
the  electron  source  or  gun,  with  the  associated  high  voltage  accelerator, 
heaters,  bias  and  deflection  power  supplies,  and  the  detector,  with  its 
associated  power,  bias,  and  deflection  supplies.  Oil  diffusion  pumps  sup¬ 
ply  the  high  vacuum  environment  for  the  electron  gun  and  the  detector. 


ft 
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A  stage  of  differential  pumping  separates  the  test  chamber  from  the 

electron  gun.  The  electron  gun  generates  a  narrow,  slightly  convergent 

beam  of  electrons  in  the  energy  range  from  15  to  50  Kev,  at  a  beam  cur- 
.  3 

rent  up  to  10  amperes.  The  electron  source  is  a  triode  gun  using  a  10 
mil  diameter  tungsten  hairpin  filament,  emission  aperture,  Wehnelt  aper¬ 
ture  and  anode  aperture.  It  was  adapted  from  the  gun  used  at  Berkeley. 

A  similar  gun,  probably  the  original  of  this  type,  was  designed  by 
39 

K.  H.  Steigci-wcilu  for  use  with  the  Tie1""  EM8  electron  microscope, 
designated  the  Kd  far-focus  cathode  system.  A  similar  gun  is  being  used 
by  Hamilton-Electrona  in  their  electron  beam  welding,  milling,  and  hole 
drilling  machines.  The  Aerospace  electron  gun  is  shown  in  Figs.  37,  38, 
and  39. 

Collimatiori  of  the  beam  is  done  with  apertures  of  small  diameter 
These  apertures  also  separate  the  scavenging  chambers  required  to  pro¬ 
vide  differential  pumping  (gas  molecules  escaping  from  the  teBt  chamber 
are  pumped  out  to  prevent,  degrading  the  vacuum  in  the  electron  gun). 

The  detector  uses  a  Faraday  cup  electron  collector  which  may  be 
biased  at  a  positive  voltage  with  respect  to  ground.  An  internal  shield 
biased  negatively  with  respect  to  the  collector  insures  suppression  of  sec¬ 
ondary  electrons. 

A  general  view  of  the  densitometer  is  shown  in  Fig.  40. 

Preliminary  tests  have  been  made  ol  the  apparatus  with  satisfac¬ 
tory  detection  of  the  scattered  beam  achieved  in  a  test  gas.  Evaluation  of 
the  preliminary  data  is  now  underway  to  establish  feasibility  and  imple¬ 
mentation  of  this  technique  to  the  combustion-rlr  iv  en  shock  tunnel. 
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'igure  39.  Electron  Bean  Gun  and  Tilting  Mechanic 


C.  High  Enthalphy  Combustion- Driven  Hypersonic.  Shock  Tunnel  Design 


1,  General  Requirements 


The  requirements  for  extending  the  performance  capability  of  the 
existing  helium- driven  hypersonic  shock  tunnel  may  be  most  easily  ob¬ 
served  by  examining  re-entry  vehicle  simulation  requirements  by  means 
of  similitude  considerations.  If  the  flight  envelope  for  future  lifting  and 
non-lifting  re-entry  vehicles  is  taken  as  that  shown  in  Fig.  41,  then  the 
simulation  requirements  for  each  regime  may  be  examined.  These  require¬ 
ments  will  be  examined  for  both  slendei  bodies,  typical  of  high  lift-drag 


ratio  vchic 


i*iuii>  bodies,  typical  <_>f  ballistic  rc-  entrv  velric.  leu 


a.  Slender-Body  Simulation  Considerations 


Regime  I: 


In  this  regime  as  a  result  of  viscous  non- equi hbi  iurn  con¬ 
siderations,  complete  thermodynamic  velocity,  and  model  suite  simulation 
are  required.  It  is  possible  in  this  regime  to  relax  the  above  requirements 
for  the  inv  iav  id  How.  i.  c.  the  actual  model  length  need  not  be  duplicated. 
Here  model  stilling  can  bo  made  by  simulating  the  loial  residence  time  of 
a  partic  le  U  /I.  In  this  (.me  the  lnvtsc’d  similitude  is  given  by 
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Since  this  regime  occurs  below  an  altitude  of  approximately  200,  000  ft 
(Fig.  41),  the  low  Reynolds  number  phenomena  associated  with  slender 
bodies  may  not  be  too  significant,  and  therefore,  the  influence  on  the 
inviscid  external  flow  can  bo  considered  small.  Reference  5  quantitatively 
indicates  the  region,  in  terms  of  the  hypersonic  interaction  parameter  y  , 
where  the  interaction  of  the  boundary  layer  with  the  external  inviscid  flow 
field  may  be  ignored.  Therefore,  if  these  effects  were  to  be  considered 
small,  the  only  need  for  reproducing  model  length  in  this  case  would  be 
for  skin  friction  and  heat  transfer  measurements.  If  the  surface  of  the 
vehicle  were  to  be  catalytic,  p-esent  theoretical  considerations  indicate 
that  non- equilibrium  phenomena  would  not  significantly  influence  either 
the  skin  friction  or  heat  transfer,  thus  offering  the  possibility  of  relaxing 
model  length  simulation  requirements. 

All  of  the  above  remarks  for  the  altitude  bounds  indicated  in 
Fig.  41  necessarily  must  be  made  on  the  basis  of  assumed  air  dissociation- 
recombination  rates  and.  therefore,  only  can  be  considered  quantitative 
within  the  accuracy  with  which  these  rates  arc  known.  There  is  also  some 
question  as  to  the  importance  of  simulating  the  ambient  free-stream  ther¬ 
modynamic  environment  (P  ,  p  ,  a  )  .  Since  the  attainment  in  practice 
of  exact  simulation  of  these  thermodynamic  parameters  represents  one  of 
the  most  difficult  simulation  problems  for  any  facility,  it  would  certainly 
be  desirable  to  relax  these  requirements  if  possible.  If  it  was  found  that, 
for  example,  the  free-stream  temperature  simulation  (within  certain  limits) 
did  not  significantly  influence  the  aerodynamic  or  heating  characteristics 
measured  on  a  model,  then  it  might  be  possible  to  simulate  the  non¬ 
equilibrium  effects  by  means  of  simulation  of  the  recombination  rate 
5 

parameter,  C.  .  Unfortunately,  this  type  of  information  is  not  available 
either  theoretically  or  experimentally  and  one  must,  at  this  time,  be 
directed  by  the  rigorous  mathematical  similitude  requirements;  \.e.,  com¬ 
plete  thermodynamic  velocity,  and  model  length  simulation. 


Regime  II  -  Frozen  Composition 

In  this  regime  the  chemical  reaction  rates  art:  considered 
slow  compared  to  either  the  rate  of  convection  along  stream  lines  or  the 
rate  of  diffusion  along  stream  lines.  A  quantitative  delineation  of  this 
regime  as  well  as  the  dissociation  non- equilibrium  regime  is  shown  in 
Fig.  41.  Mathematical  1  y ,  the  similitude  requirements  in  this  case  may  be 
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It  can  be  noted,  according  to  the  above  similitude,  that  freedom  in  the 
choice  of  both  model  length  and  possibly  free-stream  velocity  exists.  If, 
however,  we  deviate  from  the  mathemalica  1  formulation  and  examine  Ihe 
practical  aspect  of  the  simulation  requirements,  we  can  then  observe  :hat 
we  indeed  do  not  have  freedom  in  the  choice  of  free-stream  velocity  or 
model  length.  In  this  regi  :e,  we  are  interested  in  determining  the  influ¬ 
ence  of  large  amounts  of  energy  distributed  in  the  inert  degrees  of  freedom 
(particularly  dissociation).  With  this  consideration  in  mind,  coupled  with 
the  free- stream  thermodynamic  simulation  requirement,  it  is  felt  that  the 
proper  total  specific  enthalpy  level  must  be  simulated.  This  requirement 
di<  fates  simulation  of  free-stream  velocity.  We  would,  therefore,  conclude 
that,  we  have  freedom  only  in  model  length,  but  this  is  not  so  since  concurrent 

si  mu  hit  ion  of  P  ,  n  ,  u  ,  y  ,  dictates  reproduction  of  model  length. 
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Therefore,  while  the  frozen  composition  model  initially  indicates  some 
simplifications  in  terms  of  facility  simulation  requirements,  these  simpli- 
cations  cannot  he  considered  realistic  which  means  we  mutt  revert  back 
to  complete  thermodynamic,  velocity,  and  model  simulation. 

In  addition,  the  above  remarks  are  based  on  the  assumption 
that  the  dissociation  process  comes  to  equilibrium  immediately  behind  the 
shock  envelope.  If  the  rate  of  equilibration  of  the  forward  dissociation 
process  is  considered  then  at  some  point,  in  the  frozen  flow  regime,  we 
again  encounter  non- equilibrium  phenomena.  Furthermore,  since  these 
non- equilibrium  eilects  occur  at  high  altitudes,  where  the  viscous  effecis 
die  important  rnese  considerations  (viscous  non- equilibrium)  further  dic¬ 
tate  the  need  for  complete  thermodynamic  velocity  and  model  sizr 
simulation. 

It  should  be  emphasized  again  that  the  importance  of  complete 
thermodynamic  simulation  must  be  examined;  however  our  present  lack 
of  any  experimental  or  theoretical  understanding  of  this  flow  regime  does 
not  permit  <uiy  relaxation  ol  this  requirement  at  the  present  time. 

Regime  1ft  -  Low  Reynolds  Number 

'line  regime  also  requires  complete  thermodynamic,  velocity 
and  model  length  simulation  but  this  particular  item  is  separated  from 
Regime  I  for  two  reasons:  (a)  For  the  most  optimistic  case  this  regime  is 
observed  to  be  governed  as  a  minimum  by  the  complete  Navior-Stokes 
equations;  and  (b)  as  in  the  previous  regime  the  problem  of  the  time  to 
dissociate  is  present  thus  yielding  viscous  non- equilibrium  simulation 
requirements.  As  a  result  of  these  considerations  it  is  difficult  to  see 
how  in  the  foreseeable  future,  quantitative  theoretical  or  experimental 
data  will  be  available  to  relax  any  of  the  <  ompiete  simulation  requirements 
Certainly,  relaxation  of  the  simulation  requirements,  as  a  result  of  a 
more  quantitative  understanding  of  the  flow  phenomena  will  have  to  occur 
first  within  either  Regime  I  or  Regime  II  before  any  such  understanding 
can  be  possible  in  Regime  III, 
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Regime  IV  -  Perfect  Gas 
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This  regime  is  delimited  by  consirlprinn  that  if  we  tie-up  a 
relatively  small  amount  of  energy  in  dissociation  behind  the  vehicle  shock- 
envelope,  them  the  influence  of  any  "real  gas"  phenomena  on  cither  the 
aerodynamic  or  heat  transfer  problem  should  be  relatively  small.  This 
consideration  was  made  quantitatively  for  the  slender  body  problem  by 
considering  the  ratio  of  the  dissociation  energy  to  the  free-stream  kinetic 
energy,  which  is  taken  as  the  total  gas  enthalpy  lor  hypersonic,  flow.  The 
numerical  work  was  carried  out  for  a  disturbance  strength  of  30°,  This 
would  correspond  to  a  local  inclination  of  the  veh,<  in  tn.rfn.rc  of  30°  with 
respect  to  tin-  free-stream  velocity  vector.  With  those  considerations  in 
mind,  this  estimate  is  certainly  not  valid  for  the  stagnation  region;  how¬ 
ever,  it  should  provide  a  good  estimate  for  the  down- sti  earn  region  and 
tints,  for  tin*  over-all  vehicle  ae mdyna mic  and  heat  transfer  characteristics. 
In  this  region  lh<*  aerodynamic  and  he. it  transfer  characteristics  should  be 
given  to  a  pond  approximation  by 
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As  would  be  anticipated,  the  above  similitude  is  recognised  as  essentially 
that  h  r  a  perfect  gas  with  constant  specific  heat.  The  simplification  asso¬ 
ciated  with  simulation  of  the  above  compared  to  the  previous  simulation 
requirements  outlined  in  Regimes  I,  If,  and  III  can  be  readily  observed. 
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Here  we  can  relax  the  requirements  on  particle  velocity,  model  length, 
and  ambient  thermodynamic  properties.  It  is  felt  the  above  conditions  can 
be  satisfactorily  simulated  with  existing  hypervelocity  facilities  (e.  g., 
shock-tunnels,  Ohio  State,  and  University  of  California)  and,  therefore, 
should  not  be  considered  too  critical  at  the  present  time.  In  fact,  except 
for  possibly  the  very  high  altitude  conditions,  they  can  be  considered,  for 
slender  body  application,  fairly  well  understood  within  the  framework  of 
continuum  fluid  mechanics. 

Regime  V  *  Equilibrium 

This  regime  is  considered  to  be  given  by  high  Reynolds  num¬ 
ber  Reynolds  number,  chemical  equilibrium  flow  phenomena.  It  is  there¬ 
fore  felt  one  should  be  able  to  employ  the  information  obtained  from 
Regime  I  to  easily  extrapolate  to  Regime  V  for  aerodynamic  and  heat  trans¬ 
fer  data. 

Finally,  it  should  be  noted  the  above  remarks  for  Regimes  I 

through  V  arc  applicable  for  slender  bodies  down- stream  of  the  immediate 

« 

vicinity  of  the  leading  edge.  Therefore,  they  should  be  considered  perti¬ 
nent  for  slender  body  aerodynamu  and  heat  transfer  characteristics 
exclusive  of  the  stagnation  region.  The  bounds  indicated  cannot  be  con¬ 
sidered  valid  for  the  stagnation  region,  particularly  when  the  low  Reynolds 
number  and  dissociative  non- equilibrium  phenomena  are  important.  This 
problem  will  now  be  discussed  in  the  ensuing  blunt  body  similitude  consid¬ 
erations,  where  the  general  remarks  will  apply  for  the  stagnation  region 
of  a  slender  body  with  the  appropriate  dimensional  changes. 

b.  Blunt-Body  Simulation  Considerations 

The  remarks  in  this  section  will  be  directed  toward  low  fine¬ 
ness  ratio  (<  5)  lifting  and  non-lifting  re-entry  vehicles.  As  remarked 
previously,  these  considerations  will  apply  to  the  slender  body  stagnation 
regions  also. 

It  should  be  remembered  this  bound  will  be  different  for  the  blunt-body 
problem 
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Regime  I  -  Non- Equilibrium 

The  previous  remarks  made  for  the  slender  body  problem  are 
applicable  here.  However,  the  viscouB  interaction  problem  rnnsideraticns 
have  decreased  importance  for  low  fineness  configurations.  Since  we  are 
concerned  with  blunt  inviscid  dominated  flow  fields,  one  is  tempted  to 
examine  the  possibility  of  relaxing  the  velocity  simulation  requirement  by 
means  of  particle  residence  time  simulation  (i.  e.,  U^/L).  Unfortunately, 
as  in  the  previous  slender  body  case,  the  necessity  of  simulating  the  re¬ 
combination  rate  parameter  does  not  permit  us  to  take  advantage  of  this 
possibility.  Therefore,  we  conclude  that  we  must  provide  complete  simu¬ 
lation  as  in  the  slender  body  problem.  This  regime  is  particularly  impor¬ 
tant  for  advanced  re-entry  vehicles  because  maneuverability  and  the 
associated  aerodynamic  problems  must  be  considered  more  carefully  than 
those  associated  with  present  non-lifting  vehicles. 

Regime  II  -  Frozen  Composition 

Again,  the  remarks  made  for  the  slender  body  problem  apply 
here.  In  addition,  this  regime  will  be  particularly  important  for  advanced 
re-entry  vehicle  aerodynamic  (maneuvering)  and  communications  black-out 
phenomena.  Decoy  aerodynamic  design,  in  order  to  simulate  re-entry 
observables,  will  also  be  important  in  this  regime. 

Regime  III  -  Low  Reynolds  Number 

All  of  the  remarks  made  for  the  slender  body  problem  in  this 
regime  can  be  applied  here. 

Regime  IV  -  Perfect  Gas 

This  regime  is  delimited  by  considerations  similar  to  those 
previously  discussed  for  the  slender  body  problem.  Since  the  shock-layer 
associated  with  these  vehicles  will  be  dominated  by  streamlines  which  pass 
through  the  near  normal  shock,  this  regime  must  apply  to  lower  velocities 
than  those  associated  with  the  slender  body  case.  The  bound  for  the  blunt- 
body  problem  is  indicated  in  Fig.  42. 
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When  the  use  of  high-termperature,  low-density  model  testing 
facilities  for  non- equilibrium  gas  flow  problems  are  considered,  the  nature 
of  the  nozzle  expansion  process  must  be  carefully  examined.  Indeed,  the 
stagnation  conditions  required  for  low-density  model  testing  can  give  rise 
to  large  departures  from  equilibrium  in  the  nozzle  itself.  The  existence 
of  a  departure  from  equilibrium  flow  in  the  nozzle  expansion  process,  and 
the  resultant  influence  on  test  section  thermodynamic  and  flow  properties, 
must  be  quantitatively  considered  before  any  meaningful  discussion  of  aero- 
thermodynamic  simulation  can  be  made.  The  influence  of  dissociative 
non- equilibrium  flow  phenomena  in  a  hypersonic  nozzle  for  a  stagnation 
temperature  of  8000°K  and  stagnation  pressures  between  1  and  1,000  atmos¬ 
pheres  was  previously  discussed  and  numerical  results  presented  in  Figs.  16, 
17,  and  18.  For  the  more  general  simulation  considerations  discussed  here, 
a  carpet  of  the  limiting  value  of  the  frozen  degree  of  dissociation  in  the  noz¬ 
zle  (a j.)  as  a  function  of  stagnation  temperature  and  pressure,  is  shown  in 
Fig.  43. 5  The  range  4000°K  <  Tq  <  12,  000°K,  1  <  Pq  <  l,  000  atmospheres, 
is  presented.  Therefore.  Fig.  43  provides  an  estimate  of  the  required 
stagnation  conditions  to  achieve  complete  equilibrium  expansion  in  the  nozzle 
or  to  provide  a  specific  frozen  degree  dissociation  for  (controlled  non¬ 
equilibrium  expansion)  experiments.  It  can  be  readily  observed  from  this 
figure,  that  the  achievement  of  particle  velocity  simulation  with  an  equilib¬ 
rium  expansion  requires  stagnation  pressures  of  the  order  of  1,000  atmos¬ 
pheres  or  greater. 

With  the  above  simulation  and  stagnation  condition  requirements 
in  mind,  we  will  now  turn  our  attention  toward  consideration  of  the  perfor¬ 
mance  capabilities  of  the  proposed  facility. 

2.  Combustion-Driven  Shock  Tunnel  Performance 

The  objectives  in  the  proposed  design  are  to  attain  stagnation  pres¬ 
sures  greater  than  1,000  atmospheres  with  corresponding  stagnation  tem¬ 
peratures  of  8,  000°K  or  higher.  ThiB  implies  that  a  strong  primary  shock 
in  the  shock  tube  is  required,  which  can  be  obtained  by  using  high  pressure 
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helium  in  the  driver,  and  a  reasonable  dea~  '  ariver  gas  heating.  The 
'j  Ul  AV  WJi.  gas  win  tnen  be  helium  with  a  stoichiometric  mixture  of 
hydrogen  and  oxygen  to  supply  the  energy  for  heating  and  to  increase  the  • 
driver  gas  pressure  through  a  combustion  process.  In  the  operation  of  the 
combustion  driver,  however,  the  helium  is  regarded  as  the  diluent  in  sup¬ 
pressing  detonation  during  the  combustion  process.  Thus,  a  combustion 

mixture  is  usually  referred  to  by  its  helium  content  in  percent  of  the  total 
40 

volume.  Nagametsu  has  performed  experiments  for  hydrogen-oxygen- 
helium  mixtures  of  various  ratios.  Detonation  was  experienced  when  a 
65%  helium  diluent  was  used  with  a  stoichiometric  mixture  of  hydrogen  and 
oxygen.  Difficulties  in  igniting  the  ga.s  were  encountered  when  an  80% 
helium  dilution  was  used. 

In  calculating  the  performance  of  the  shock  tunnel,  helium  is  re¬ 
garded  as  the  drivei  gas.  The  product'of  combustion.,  vapor  of  H^O,  is 
then  regarded  as  the  contaminant  in  the  driver  gas.  Since  the  water  vapor 
contributes  to  an  appreciable  increase  in  molecular  weight  oi  the  final  gas 
mixture  over  that  of  the  pure  helium  its  effect^  cannot  be  ignored  com¬ 
pletely  in  the  calculations.  A  detailed  study  of  the  combustion  process 
would  give  an  accurate  prediction  oi  the  driver  gas  conditions  with  regard 
to  temperature,  pressure,  speed  oi  sound  ape*  if ic  heat  ratio,  etc  A 
detailed  study  of  this  nature  in  order  to  be  adequate  for  our  present  pur¬ 
poses.  should  trace  the  expansion  process  of  the  driver  gas  after  the  shock 
formation  has  been  initiated.  This  implies  the  necessity  for  the  construc¬ 
tion  of  portions  of  a  M.olher  diagram  for  each  mixture  ratio  under  consid¬ 
eration,  in  order  that  the  regions  of  interest  can  be  adequately  covered. 

For  our  present  purpose,  the  equilibrium  gas  temperature,  pressure,  and 
specific  heat  ratio  y  were  determined  at  the  time  when  the  combustion 
process  was  completed.  Isentropic  expansion  of  the  driver  gas  was  assumed 
in  order  to  calculate  the  conditions  behind  the  interface.  Four  mixture 
ratios  were  considered  in  calculating  the  driver  gas  conditions,  and  subse¬ 
quently.  the  shock  tunnel  performance  with  regard  to  stagnation  conditions 
in  the  region  after  the  reflected  primary  shock.  It  is  on  this  stagnation 
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A  brief  description  of  the  process  involved  is  in  order  when 
formulating  the  problem.  Four  mixture  ratios  were  included  to  cover  the 
possible  operating  range  between  no  ignition  and  detonation.  A  constant 
volume  combustion  process  was  considered.  In  each  case,  the  driver  gas 
was  taken  as  the  gas  mixture  after  complete  combustion,  assumed  homo¬ 
geneous  and  at  equilibrium.  Instant  diaphragm  removal  was  considered  at 
a  time  when  the  combustion  pressure  would  be  at  its  peak.  It  was  assumed 
further  that  the  expansion  process  in  the  driver  gas  was  isentropic.  Using 
the  expression  relating  the  shock  tube  pressure  ratio  to  the  primary  shock 
Mach  number  we  have 


P4  P2(P_4\ 
Pi  '  Pi  IP 3/ 


(152) 


whcr  e 


^4 


by  taking  y^  equal  to  some  constant,  the  pressure  ratio  can  be  readily 

obtained  for  a  given  shock  Mach  number.  Figure  44  is  a  plot  of  Eq.  (152) 

for  the  four  mixture  ratios.  During  this  analysis,  extensive  use  was  made 

41 

of  the  hypersonic  gas  dynamic  charts  for  equilibrium  air  by  Feldman. 

After  selecting  the  shock  Mach  number,  Ms,  various  parameters 
with  regard  to  pressure  and  temperature,  etc.,  can  be  obtained  for  differ¬ 
ent  regions  in  the  shock  tube  A  plot  of  the  reservoir  conditions  versus 
shock  Mach  number  is  given  in  Fig.  43.  Figure  46  is  a  graph  of  the  tunnel 
reservoir  conditions  for  different  shock  tube  initial  conditions. 

By  superimposing  the  shock  tube  conditions  onto  the  chart  for  the 
non- equilibrium  regime,  the  region  of  coverage  can  be  obtained.  This 
coverage  is  presented  in  Figs.  4?  and  48. 
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■'igure  44,  Combustion -Driven  Shock  Tube  Performance  as  a 
Function  of  Percent  Helium  Diluent  in  the  Stoichiometric 
Mixture  of  Hydrogen  and  Oxygen 
(plot  of  Eq,  152) 
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Figure  47,  Variation  of  with  Stagnation  Conditions  for  Air 
Showing  Region  of  Coverage  with  Shock  Mach  Numbers 
Attainable  in  the  Shock  Tube 
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Calculations  for  a  shock  tunnel,  using  pure  hydrogen  as  the  driver 
gas,  is  given  in  Fig.  49  as  a  comparison.  Hydrogen  is  very  efficient  for 
use  as  a  driver  gas;  however,  its  inherent  dangerous  handling  character¬ 
istics  make  it  unattractive  for  use  in  the  present  facility. 

a.  Area  Ratio  Effect  in  Shock  Tunnels 

The  preceding  analysis  was  appropriate  for  an  ideal  shock 
tube  with  constant  area  throughout.  For  shock  tunnels  having  unequal 
driver  and  driven  tube  diameters  and  a  throat  section  at  the  diaphragm 
location,  the  performance  calculated  previously  would  be  modified  to  com¬ 
pensate  for  area  ratio  effect  and  the  restriction  at  the  throat.  A  shock  tube 
with  such  a  restriction  is  given  schematically  (see  Fig.  50).  The  effective 

4j> 

pressure  ratio  for  such  a  configuration  would  be 


where 
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is  the  gain  factor  modifying  the  actual  pressure  ratio  to  an  effective  pres¬ 
sure  ratio.  With  the  boundary  conditions  such  that  *  1  and  if 

A.,  /A,*  is  given,  then  M,  can  be  obtained  by  solving  the  equation 
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Figure  49.  Variation  of  Qf  with  Stagnation  Conditions  for  Air. 
Region  of  Coverage  Using  Pure  Hydrogen  as  the 
Driver  Gas  is  Indicated 
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Again,  can  be  obtained  by  solving  the  following  equation: 


~2  ♦  - 11  Tngm 

_2  +  (,4  .  1|M|J 
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3b 


(156) 


A  plot  of  the  gain  factor  g  versus  diameter  d^  is  given  in  Fig.  51  where  a 
fixed  ratio  d^/d^  =  1.5  is  used  (existing  configuration  for  the  Aerospace 
Corp.  shock  tunnel). 

b.  Shock  Tunnel  Testing  Time 

The  testing  time  for  a  hypersonic  shock  tunnel  is  determined 
by  the  length  of  time  that  the  reservoir  pressure  remains  constant:  after 
the  primary  shock  has  been  reflected  from  the  end  of  the  shock  tube. 
Theoretically,  this  time  interval  is  the  time  difference  between  the  pri¬ 
mary  shock  arrival  at  the  end  of  the  shock  tube  and  the  reflection  of  the 
first  reflected  shock  wave  at  the  interface. 

The  ideal  testing  time  for  a  non-tailored  shock  tunnel  is  a 
function  of  the  shuck  wave  speed,  shock  tube  length,  and  given  initial  con¬ 
ditions.  Figure  52  shows  a  plot  of  the  testing  time  as  a  function  of  shock 
Mach  number  with  driven  tube  length  as  a  parameter.  The  limit  for  a  dif¬ 
ferent  driver  length  is  also  indicated.  The  significance  of  this  limit  is  that 
for  a  given  driver  section  tube  length  the  test  time  will  decrease  progres¬ 
sively  as  the  shock  Mach  number  is  decreased  below  the  limit.  Since  ideal 
conditions  will  not  be  encountered  in  practice,  the  actual  testing  time  will 
be  some  fraction  of  the  ideal  testing  time.  A  rough  rule  of  thumb  gives 
this  fraction  to  be  approximately  0.5. 

43 

Figure  53  shows  the  regime  for  tailored  operation  at  differ¬ 
ent  reservoir  conditions  with  shock  tube  initial  conditions  at  parameters. 
Figure  54  shows  the  driver  pressure  versus  shock  Mach  number  with  the 
shock  tube  initial  conditions  as  parameters. 
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The  calculation  so  far  has  been  applied  to  a  stoichiometric 

mixture  of  hydrogen  and  oxygen  with  various  percentages  of  helium  diluent. 

It  is  theoretically  feasible  to  extend  the  tailoring  regime  by  adding  one  or 
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more  gases  of  different  molecular  weight  to  the  combustible  mixture. 

Such  a  possibility  can  be  explored  for  combustion-driven  operation  should 
it  be  desired. 

1  Figure  55  shows  the  tailoring  regime  superimposed  on  a  chart 
of  the  reservoir  conditions.  The  upper  limits  of  a  non-tailored  tunnel  op¬ 
eration,  having  a  1-3/4  millisecond  testing  time,  are  indicated  tor  various 
driven  section  lengths.  Figure  56  is  a  plot  of  the  theoretical  ideal  testing 
time  possible  for  tailored  operation,  as  a  function  of  primary  shock  Mach 
number,  with  driven  section  tube  length  aa  a  parameter,  The  limiting  line 
of  zero  advantage  of  tailoring  is  indicated  for  a  20-foot  driver.  This  means 
that  tailoring  will  be  useful  only  above  certain  shock  Mach  numbers  if  the 
tube  length  is  given.  It  is  also  noted  that  the  usefulness  of  tailoring  will  be 
more  pronounced  for  shock  tunnels  with  short  driven  section  tube  lengths. 

c .  Teat  Section  Conditions 

Once  the  conditions  behind  a  reflected  shock  have  been  estab¬ 
lished,  the  flow  parameters,  as  a  function  of  nozzle  expansion  ratio,  may 
be  determined  by  means  of  techniques  discussed  in  II.  A.  These  calculations 
have  beon  made  for  reflected  shock  conditions  which  may  be  attained  by  the 
combustion  driver.  All  the  test  conditions  possible  will  not  be  presented 
here,  but  rather,  the  results  for  high  enthalpy  reservoir  conditions  are 
presented  (Fig.  57).  In  this  figuro,  the  Mach  number  is  presented  as  a 
function  of  nozzle  expansion  ratio  for  a  stagnation  pressure  of  1000  atmos¬ 
pheres  and  a  stagnation  temperature  of  8720°K.  The  results  are  presented 
for  both  an  assumed  equilibrium  expansion  and  the  non- equilibrium  expan¬ 
sion  as  predicted  by  the  analysis  outlined  in  11.  A.  The  non- equilibrium 
expansion  will  freeze  out  approximately  11  perdent  of  stagnation  enthalpy 
in  the  nozzle  expansion  process.  From  this  figure,  it  can  be  observed  that 
the  Mach  number  for  the  assumed  equilibrium  expansion  is  approximately 
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27  at  an  expansion  ratio  of  10^,  while  the  corresponding  non -equilibrium 
Mach  number  is  approximately  40.  The  expansion  ratio  of  10^  will  be 
attainable  in  the  modified  facility  as  a  result  of  the  inclusion  of  a  six-foot 
diameter  nozzle  (the  nozzle  design  and  detailed  modification  is  presented 
in  the  next  section,  "Description  of  Modified  Facility").  The  performance 
shown  in  Fig.  57  will  allow  experimental  investigation  of  the  use  of  con¬ 
trolled  non- equilibrium  expansions  discussed  in  II.  A.  The  particle  velocity 
at  this  condition  will  be  approximately  19,  300  feet  per  second.  The 
Reynolds  number  and  density  altitude  simulation  at  this  expansion  ratio  (10^) 
corresponds  to  flight  at  approximately  300,000  ft.  Thie  condition  is,  there¬ 
fore,  ideal  for  Investigating  both  low  Reynolds  number  and  dissociative 
non- equilibrium  flow  phenomena. 

d .  Design  Considerations  for  the  Combuatlon-Orivcn  Shock  Tunnel 

In  designing  the  Aerospace  Corporation  combustion-driven 
shock  tunnel,  the  foregoing  analysis  was  applied  to  determine  the  final  con¬ 
figuration.  Peripheral  factors,  howovei,  such  as  available  laboratory 
apace,  played  important  roles  in  the  ovor-all  consideration.  Compromises 
have  had  to  be  made  with  regard  to  the  ultimate  in  performance  as  against 
a  more  practical  design  limited  by  mechanical  or  economical  considerations. 

The  configuration  arrived  at  is  a  shock  tunnel  having  a  driver 
section  length  of  20  feut  and  a  drivon  section  length  of  30  feet  While  both 
sections  have  an  inside  diameter  of  three  inches,  the  major  diaphragm  lo¬ 
cation  has  a  neck  down  region  with  a  diameter  of  two  inches. 

The  advantage  in  using  a  smaller  diameter  on  the  diaphragm 
station  is  apparent  when  mechanical  properties  of  the  diaphragm  are  taken 
into  consideration,  since  a  larger  diameter  diaphragm  is  relatively  unre¬ 
liable,  especially  at  the  hish  pressure?  under  consideration.  The  present 
diaphragm,  which  has  been  well  developed  in  use,  will  be  retained  in  the 
new  configuration.  From  a  gasdynamics  standpoint,  a  throat  would  intro¬ 
duce  a  local  steady  flow  process  which  would  alleviate,  to  some  extent,  the 
expansion  wave  propagation  from  the  driver  towards  the  end  of  the  shock 
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tube.  This  effect  would  tend  to  increase  the  test  time  of  the  shock  tunnel. 

The  designed  driver  operating  pressure  is  2000  atmospheres,  with  allow¬ 
ance  for  a  faetor  of  five  increase  in  the  operating  pressure  in  the  event  of 
detonation.  An  over-all  safety  factor  of  1.1  is  then  used  to  calculate  the 
ultimate  yield  strength  of  the  driver.  The  driven  section  was  designed  to 
withstand  an  internal  pressure  of  20,000  psi.  Since  laboratory  space  is 
restricted,  a  30-foot  driven  section  length  was  chosen.  Tailored  operation 
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Maximum  use  is  made  of  the  existing  shock  tunnel  hardware.  The  intended 
driver  gas  mixture  ratio  to  be  used  is  between  70  and  75  percent  helium 
Within  the  tailored  operation  regime  tho  maximum  stagnation  pressure, 
after  the  reflected  shock,  is  approximately  1250  atm.  The  maximum  stag¬ 
nation  temperature  is  approximately  9800(,K,  and  a  theoretical  testing  time 
of  about  four  milliseconds  is  anticipated. 

3 .  Description  of  Mi >d i f i ed  Shock  Tunnel  Facility 
a.  Introduction 

The  new  tailored  combustion-driven  shuck  tunnel  is  being 
constructed  in  order  to  achieve  the  high  stagnation  temperatures  and  pres¬ 
sures  and  the  required  lest  tunes  spec  tfiod  in  Sections  III  C  1  and  III  C  2 
of  this  report 

According  to  the  length  requirements  spec  ified  in  Sec  lion  III  C  .2, 
a  new  driver,  20  feet  long  and  three  inches  in  diameter,  is  to  be  constructed. 
The  driven  section  will  cuiisist  of  the  three  sections  of  the  cold  helium 
driven  shock  tube  (30  feet).  A  new  nozzle  and  a  six-foot  test  section  is  to 
be  constructed,  utilizing  the  existing  vacuum  dump  tank.  Two  cylindrical 
extensions  will  be  added  to  the  vacuum  tank.  The  nozzle  will  be  fabricated 
of  plastic  sections  of  varying  lengths  in  order  to  permit  flow  Mach  numbci 
variation.  The  nozzle  and  test  section  are  enclosed  in  the  dump  tank.  The 
design  Mach  number  range  is  from  16  to  32.  A  schematic  of  the  shoc  k 
tunnel  facility,  as  it  will  be  installed  in  the  laboratory  building,  is  shown 
in  Fig,  58.  A  detailed  description  of  the  driver  section  and  gas  handling 
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lew  of  Shock  Tunnel  Showing  12-inch  Test  Section 


system,  the  driven  section,  and  the  noszle  test  section  dump  tank  is  given 
below: 


% 


b.  Driver  Section  and  Gas  Charging  System 

Due  to  the  possibility  of. detonation  of  the  combustion  mixture, 
the  driver  section  is  designed  for  an  ultimate  pressure  of  165,  000  psi.  The 
working  pressure,  in  general,  will  not  exceed  30,000  psi.  The  temperature 
is  expected  to  reach  between  2000°  and  30QQ°K,  but  it  is  transient,  as  is 
the  high  combustion  pressure. 

In  order  to  meet  these  requirements,  the  driver  section  and 
its  closely  associated  fittings  will  be  fabricated  from  forged  300M  steel  and 
heat  treated  after  rough  machining  to  maximize  the  effectiveness  of  the 
neat  treating.  The  vessel,  rap  hatching  mechanism,  and  diaphragm  pack¬ 
age  are  detailed  iu  Fig.  5‘). 

The  gas  handling  system  for  loading  the  drive.r  with  a  com¬ 
bustible  mixture  of  oxygen  and  hydrogen,  diluted  with  helium,  is  shown 
schematically  in  Fig.  60.  A  gas  accumulator,  with  a  volume  equal  to  that 
of  the  driver  section,  acts  as  a  fast  high-pressure  hydrogen  pump.  A  stoi¬ 
chiometric  mixture  of  hydrogen  and  oxygen  diluted  with  75  percent  helium, 
at  a  total  presume  of  4000  psi  before  combustion,  is  used,  Oxygen  is 
pumped  into  the  driver  sect  ion  from  the  storage  marked  O.,  to  a  pressure 
of  330  psi.  Helium  is  then  added  until  a  pressure  of  3130  psi  is  reached. 
Then  hydrogen  is  pumped  into  the  hydrogen  accumulator  to  a  pressure  of 
660  psi.  Helium  gaB  then  passes  into  the  accumulator  behind  the  piston. 

As  the  helium  is  stored  at  very  high  pressures  (15.000  -  20,000  psi),  it 
forces  the  hydrogen  into  the  driver.  This  sequence  of  loading  never  per¬ 
mits  a  mixture  of  hydrogen  and  oxygen  to  exist  in  richer  concentrations 
than  the  mixture  to  be  employed  in  the  test.  Hence,  the  probability  of  deto¬ 
nation  occuring  during  the  filling  process  is  low. 

Proper  distribution  and  circulation  of  the  driver  charge  is 
obtained  by  using  a  long  ported  stainless  steel  tube  inside  the  driver  to 
serve  as  the  fill  tube.  The  port  area  varies  exponentially  along  the  tube. 
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iguro  60.  Gas  Handling  System  for  Charging  the  Driver  Section 


coAnjj  greatest  at  the  end  o.;  fill  tube.  In  this  way,  the  gas  enters  the 
driver  section  uniformly  a.ong  the  entire  length.  After  the  driver  is 
charged,  the  gas  is  ignited  by  exploding  a  long,  thin  (0.0015  in.  }  copper 
wire  alloy  parallel  to  and  above  the  fill  tube.  The  insulated  high  voltage 
terminal  is  in  the  end  cap  and  the  Hii  tube  acts  as  a  ground  return.  The 
wire  is  exploded  by  the  discharge  ux  a  5  mid  capacitor  charged  to  30,  000 
volts.  The  exploding  wire  mothod  was  chosen  as  the  most  reliable  ignition 
technique.  Pilot  experiments  show  that  spark  plugs  fail  to  operate  satis  ¬ 
factorily  at  the  high  initial  pressure  of  the  driver  charge.  The  driver 
vessel  thus  requires  no  spark  plug  ports  In  the  side  walls,  reducing  design 
complexity. 

The  diaphragm  shown  in  Fig.  61  is  inserted  into  its  package 
and  fitted  in  place  against  the  driver.  The  locking  nut  with  interrupted 
threads  slides  into  position  on  the  driver  and  is  rotated  into  the  closed  po¬ 
sition.  The  "O"  rings  around  the  diaphragm  package  insure  innate  sealing. 
After  the  diaphragm  and  ignition  wire  are  installed  and  the  driver  charged 
with  gas,  the  shmk  lube  is  ready  to  fire.  When  the  gas  mixture  has  burned, 
a  mixture  of  82  percent  helium  and  18  perc  ent  water  vapor  forms  the 
driver  gf.s, 

c.  Driven  Sec  turn 

The  new  driven  section  consists  of  the  pxistinn  driven  section 
and  the  existing  driver,  resulting  in  a  total  driven  section  30  feet  in  length, 
Ihis  unit  is  moved  forward  bIx  feel  to  a  new  position  and  the  main  buffer  or 
anchor  to  fix  the  position  of  the  shock  tube  is  fastened  as  shown  in  Fig.  56. 
All  pressure  ports  used  in  the  existing  shock  tube  are  available  for  the 
instrumentation  of  the  drive  section  with  the  exception  of  the  last  down¬ 
stream  position.  This  provides  the  vacuum  pump  connection  for  operating 
the  driven  section  at  initial  pressures  of  less  than  one  atmosphere. 

d .  Tunnel  and  Test  Section 

Beyond  the  buffer  plate  which  immobilizes  the  driven  section, 
a  straight  length  of  shock  tube  arts  as  a  mechanical  transition  section  fo* 
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connecting  to  the  tunnel  (Fig.  5S).  The  nozzle  consists  of  a  two-inch  thick 
plate  containing  the  throat  block  which  is  bolted  directly  to  the  dump  tank. 
The  conical  expansion  section  and  the  test  section  are  contained  entirely 
in  the  vacuum  dump  tank  (Fig.  62).  The  nozzle  sections  do  not  encounter 
the  pressure  load  of  an  evacuated  chamber. 

The  30°  conical  nozzle  is  formed  by  truncated  cones  of  vary¬ 
ing  lengths,  The  flow  Mach  number  is  thus  changed  by  adding  or  removing 
sections,  The  "test  section"  is  thus  the  exit  of  the  largest  truncated  coni 
cal  section,  A  reinforced  plastic  laminate  is  used  to  fabricate  the  nozzle 
sections.  In  order  to  accommodate  the  nozzle  and  test  section,  the  dump 
tank  is  enlarged  by  adding  two  midsections  to  the  dump  tank. 

The  end  section  of  the  dump  tank  is  mounted  on  tracks  and 
permitted  to  slide  open,  allowing  easy  arress  to  the  model  and  internal 
instrumentation.  Ten-hub  openings  in  the  dump  tank,  at  points  corres¬ 
ponding  to  the  test  section  of  each  truncated  cone,  arc  provided  for 
schlieren  or  other  observation.  The  model  support  enters  the  dump  tank 
at  the  fi*nt  opening  downBtrcnm  of  the  test  section  window.  The  model 
support  ifl  scpaiuted  from  the  tunnel  striu  turc  and  sealed  by  a  rubber  dia¬ 
phragm  where  it  enters  the  dump  tank  as  in  the  existing  shock  tunnel. 

e.  Com  lusion 

The  new  hypet  sonic  shock  tunnel  with  the  combustion  driver 
will  provide  the  ui<  reused  capability  required  to  support  present  and  future 
AF  systems.  The  reservoir  stagnation  condition  will  be  increased  from  a 
maximum  of  J2U0°K  and  250  atm  to  980Q°K  and  1600  atm  This  corres¬ 
ponds  to  an  increase  in  the  flow  velocity  at  the  test  section  of  from  9000 
feet  per  second  to  21,000  feet  per  second.  'Ihe  Mach  number  which  can  be 
achieved  in  the  test  section  oefore  condensation  occurs  is  ini  teased  from 
20  to  35.  The  new  six-foot  test  section  will  make  possible  an  im  tease  in 
the  model  size  and  the  instrumentation  capacity  of  each  model 


4.  Combustion  Experiments 


The  primary  modification  of  the  current  facility  consists  of  the 
design  and  Incorporation  of  a  combustion  driver.  A  stoichiometric  mix¬ 
ture  of  oxygen  and  hydrogen  is  ignited  at  high  pressure  within  the  driver 
to  heat  the  primary  driver  gas  (helium).  The  initial  loading  pressure  of 
the  driver  is  approximately  4000  psi  and  the  pressure  after  ignition  reaches 
values  as  high  as  'JO,  000  pas.  In  order  to  provide  for  safe  operation  of  the 
facility,  however,  the  design  point  of  the  driver  must  be  over  150,000  psi 
to  pruvide  sufficient  margin  in  the  event  of  detonation  within  the  driver. 
Since  all  of  the  current  facilities  employing  combustion  drivers  operate  at 
considerably  lower  pressures  it  was  necessary  to  undertake  an  experimen¬ 
tal  investigation  of  the  combustion  processes  of  such  a  mixture  at  the  pres¬ 
sures  of  interest  as  well  as  methods  of  uniform  ignition  of  the  mixture 

The  experimental  program  undertaken  im  ludes  the  following 
investigations: 

a  Method  of  ignition 

1.  sparkplug 
2  exploding  wire 

b  Combustion  proi  esses  with  regard  to 

l .  nuxtur e  ratio 
l.  effoit  of  loading  pressure 
i.  i  ombuslitm  pressur  e 
4,  detonation 

t 

Spark  plugs  have  been  routinely  operated  sun  essful.ly  for  sho  k  tunnels 
with  a  loading  pressure  in  the  range  600  to  1200  psi  and  a  i  or  i  esponding 
pressure  rise  after  combustion  to  6000  -  12,000  psi.  This  pressute  is 
within  the  design  limitations  of  commercially  available  spark  plugs.  With 
an  increase  in  loading  pressure,  the  pressure  after  ignition  may  ext  cod  the 
mechanical  stress  limit  of  commercial  spark  plugB. 
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Another  major  consideration  is  the  breakdown  potential  required 
for  the  spark  to  bridge  the  electrode  gap  at  high  pressure  in  a  dominantly 
helium  atmosphere.  For  minimum  achievable  electrode  spacinga,  the  log¬ 
arithmic  rise  of  breakdown  potential  with  pressure  requires  voltages  of  the 
order  of  1Q0KV  at  4000  psi.  Extremely  small  values  of  electrode  spacing 
limits  spark  plug  life  to  a  single  run. 

The  first  method  of  igniting  the  gas  mixture  investigated  was  a 
spark  plug  ignition  system,  in  view  of  the  low  pressure  operational  success 
with  this  technique  in  other  facilities.  Specially  designed  plugs  were  con¬ 
structed,  since  commercially  available  spark  plugs  are  not  designed  for 
4000  psi  operation.  An  assembly  drawing  of  the  high-pressure  plug  is 
shown  in  Fig.  63.  Figure  64  is  a  photograph  of  the  finished  spark  plug. 

The  valve  body  and  the  central  electrode  were  made  of  stainless  steel  md 
the  ends  of  the  electrode  were  machined  flat.  The  filler  material  is  teflon. 
A  gas-tight  assembly  is  achieved  when  the  retaining  nut  is  scr  ewed  on 
tightly  from  the  top.  compressing  the  insulating  teflon.  The  0.010  inch 
elertrnH»  gap  is  the  internal  annular  area.  A  ring  insert  at  the  tip  iorms 
the  ground  electrode  and  permits  adjustment  of  the  gap  distant  e.  The  shar 
radii  of  the  flat  electrodes  provide  high  field  strength  for  initiation  of 
breakdown.  The  plug  was  hydraulically  tested  to  a  statu  ptessurc  ol 
10  000  psi  without  failure 

Initial  combustion  tests  were  made  with  this  spark  plug  configura¬ 
tion  in  a  s pet  i  dly  designed  test  chamber.  Gas  charge  pressures  up  to  800 
psi  were  successfully  ignited  by  the  spark  plug  with  a  potential  of  15-20  KV 
across  the  electrode  gap.  Mechanical  failures  occasionally  dec  eloped  in 
the  electrode  insert.  Experiences  encountered  elsewhere  show  electron e 
pitting  during  repeated  firing  at  high  pressure,  resulting  in  an  im  reased 
effective  electrode  spacing.  Gap  distance  increase  by  pitting  can  be  sulfi- 
da*  to  cause  erratic  or  incomplete  breakdown  at  pressures  onsiderably 
below  4000  psi. 
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Figure  63.  Combustion  Driver  Spark  Plug 
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For  the  combustion  driver,  18  spark  plugs  must  be  used  along  the 
length  of  the  driver  for  uniform  combustion.  The  probability  of  one  or 
more  plugs  malfunctioning  during  any  given  run  is  high.  This  is  obviously 
undesirable  as  the  possibility  of  induced  detonation  and  erratic  performance 
increases,  The  problems  involved  in  removing  18  spark  plugs  with  high- 
pressure  fittings  for  routine  chocks  at  regular  intervals  do  not  present  an 
attractive  picture  from  an  operational  standpoint. 

Another  important  design  factor  is  the  acute  stress  problem  aris¬ 
ing  in  the  mechanical  design  of  the  high-pressure  driver  due  to  the  spark 
plug  holes,  These  factors  presented  such  difficulties  that  other  ignition 
techniques  were  considered  for  the  hlgh-proasuro  combustion  driver. 

Ah  a  result  of  these  considerations,  investigations  were  undertaken 
on  exploding  wire  techniques,  Successful  operation  has  been  obtained  in 
many  laboratories  using  this  approach  but  at  lower  gas  charge  pressures 
than  4000  psi.  The  exploding  wire  considered  for  use  in  the  Aerospace 
combustion  driver  consists  of  a  series  ol  wires  spanning  the  entiup  length 
of  the  driver  in  order  to  obtain  uniform  combustion  in  both  the  axial  and 
radial  directions.  Suite  'nitintiuu  ol  the  explosion  is  a  prn< ess  involving 
oniy  the  wire,  there  are  tu>  problems  associated  with  electrical  breakdown 
of  the  gas  and  cleetrode  erosion.  1  he  potential  required  to  initiate  wire 
explosion  is.  tor  all  prni  Meal  purposes,  independent  of  gas  prossui  e. 

The  equipment  used  in  the  experiment  consisted  simply  of  two  3pf 
capacitors  connected  In  series  with  the  wire,  ’the  capacitors  ate  rated  at 
15KV  each,  and  a  total  enpat  itance  of  l,5pf  is  obtained  at  10  KV.  The 
schematic  is  shown  in  Fig.  65. 
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Figure  65.  Schematic  of  Exploding  W ire  System 

With  the  capacitors  charged  t.i  the  desired  potential,  argon  was 
allowed  to  flow  into  the  switch.  A  minimum  potential  of  12  KV  is  required 
to  cause  electrical  breakdown  in  the  argon  for  the  gap  spacing  used  in  the 
switch  (0.5  inch  between  flat  faced  rods).  The  ionized  path  formed  in  the 
argon  completes  the  circuit,  dincharging  the  stored  energy  into  the  wire, 
resulting  in  a  violent  explosion. 

Tile  circuit  discharge  characteristic,  derived  from  the  lmpcdnnt  c 
of  the  circuit  component*  (im  hiding  the  exploding  wire)  is  i  ritu  al  for  a 
short  duration  discharge  resulting  in  a  strong  explosion.  When  the  circuit 
is  completed  by  ionization  o'  the  argon  switch,  current  begins  to  ilow 
through  the  wire.  Due  to  the  wire  resistance,  energy  dissipation  within 
the,  wire  occurs  in  the  form  ot  joule  heating.  The  amount  of  current  pass¬ 
ing  through  the  wire  is  determined  by  the  total  resistance  of  the  wire  and 
the  charge  potential.  Thus,  a  longer  wire  requires  a  higher  threshold  p0‘ 
tential  to  obtain  the  required  current  for  a  rapid  explosion.  Changing  the 
wire  gauge  and  material  influences  the  threshold  potential.  As  the  wire  ie 
raised  rapidly  to  high  temperatures,  vaporization  is  initiated.  Present 
theory  and  experiments  show  that  metal  droplets  are  formed  initially,  and 
as  further  heating  occurs,  the  droplets  eventually  form  vapors  and  expand 
radially  in  the  form  of  a  cylindrical  shock  wave.  During  the  droplet 


formation  period,  the  electrical  circuit  can  be  interrupted  due  to  physical 
separation  of  the  conducting  material.  If,  however,  the  separation  of  the 
droplets  and  sub-droplets  is  below  a  critical  value,  space  charge  effects 
will  maintain  the  electrical  current,  provided  the  residual  potential  is  high 
enough.  Further  heating  creates  a  plasma  which  continues  to  conduct  cur¬ 
rent  until  dissipation  due  to  the  radial  expansion  process  quenches  the 
discharge.  Each  step  of  these  processes  occurs  in  a  very  short  time  inter¬ 
val,  determined  by  the  material,  the  length,  and  the  uniformity  of  the 
cross  section  of  the  wire.  The  discharge  time  constant  is  important  in  that 
it  must  be  short  or  comparable  in  magnitude  to  the  effective  separation  time 
of  the  droplets  or  the  dissipation  time  of  the  plasma. 

Experiments  have  been  performed  to  seek  the  threshold  charge 
potential  required  to  produce  a  strong  explosion  for  varying  lengths  of  wire. 
In  addition,  different  wire  gauges  and  materials  were  investigated  These 
included  0.005  inch  tungsten,  0.001  inch  copper,  and  0.0005  inch  platinum. 
The  most  successful  wire  found  to  date  is  0.0015  inch  copper  alloy.  Wires 
oi  this  material  are  easy  to  handle  and  have  high  mechanical  strength.  A 
single  wire  up  to  seven  feet  in  length  has  been  successfully  exploded  with 
the  present  power  supply.  The  results  of  this  series  of  experiments  are 
presented  in  Fig. 66.  A  linear  relationship  is  noted  between  charge  poten¬ 
tial  and  wire  length.  Extrapolating  the  results  to  a  20-foot  wire  indicates 
a  charge  potential  of  about  75KV  is  required  for  a  strong  explosion.  This 
creates  a  problem  in  the  driver  design,  since  the  close  spacing  within  the 
driver  means  that  arcing  to  the  driver  wall  is  difficult  to  suppress  at  this 
high  poten.ial.  In  addition,  the  lead  wire  insulation  must  withstand  high 
voltage  and  high  pressure  simultaneously. 

A  possible  solution  is  to  divide  the  wire  in  three  sections  and  con¬ 
nect  them  in  parallel.  In  this  way,  the  breakdown  potential  is  considerably 
reduced.  The  problems  associated  with  simultaneous  ignition  of  multiple 
wires  is  now  being  investigated.  The  experimental  apparatus  is  shown  in 
Figs.  6  7  and  68.  Preliminary  combustion  experiments  with  the  exploding 
wire  were  conducted  by  installing  a  wire  on  a  spark  plug  modified  for  the 
purpose  (Fig  69). 
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Figure  67.  Explosion  of  -’ll ree  Wires  in  Parallel 
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A  strong  explosion  provides  adequate  margin  for  igniting  the  gas 
mixture  at  high  pressure  and  eases  ignition  problems  associated  with  mix¬ 
ing  non-uniformities.  A  closer  study  df  the  ignition  problem  indicates  that 
complete  burning  (non-explosive)  of  the  wire  may  be  adequate  to  ignite  a 
high-pressure  gas  mixture.  Recent  reports  from  other  laboratories  show 

that  burning  of  aluminum  tape  has  been  used  successfully  to  ignite  gas  mix- 

■ 

tures  with  pressure  levels  comparable  to  the  combustion  driver. 

A  specially  designed  small  scale  pilot  combustion  chamber  was 
constructed.  The  chamber  volume  is  1.5  cubic  inches  with  a  minimum  wall 
thickness  of  one  inch.  The  cylindrical  chamber  has  a  spark  plug  at  one 
end,  and  a  diaphragm  seal  at  the  opposite  end  of  the  chamber.  The  dia¬ 
phragm  is  calibrated  to  break  at  about  5500  psi.  This  enables  a  mixture 
gas  at  4000  psi  to  be  loaded  into  the  chamber.  The  pressure  increase  after 
ignition  is  such  that  the  diaphragm  rupture  pressure  is  exceeded,  resulting 
in  diaphragm  break  and  chamber  pressure  drop. 

The  combustion  chamber  is  connected  to  bottles  of  hydrogen  and 
oxygen.  Helium  supply  is  routed  from  the  main  reservoir  of  the  present 
shock  tunnel.  A  schematic  of  the  installation  is  shown  in  Fig.  70.  The 
sequence  of  operations  begins  with  flushing  the  chamber  with  oxygen,  fol¬ 
lowed  by  loading  to  the  desired  pressure.  Hydrogen  is  injected  next  by 
setting  the  regulator  to  the  correct  pressure  and  then  opening  the  hand 
valve.  Each  hand  valve  must  be  closed  immediately  after  loading  each 
component  gas  in  order  to  prevent  feed  back  of  gas  mixtures  into  the  com¬ 
ponent  gas  lin^s.  Helium  from  the  shock  tunnel  reservoir  is  injected  last. 
Difficulty  in  obtaining  the  correct  mixture  ratio  in  the  pilot  combustion 
chamber  has  been  encountered  due  to  the  excessive  volume  in  the  connect¬ 
ing  gas  lines.  Since  a  finite  length  is  required  to  connect  the  test  chamber 
and  the  end  valves,  the  connecting  line  volume  must  be  eliminated  or 
accounted  for. 

Solenoid  valves  were  installed  directly  onto  the  test  chamber  to 
eliminate  this  excess  volume.  Good  results  were  obtained  in  this  manner 
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Figure  7  0,  Schematic  of  Pilot  Combustion  Test  Facility 
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but  the  solenoid  valves  were  damaged  repeatedly  at  pressures  much  below 
the  rated  operating  pressure.  The  test  combustion  chamber  volume  was 
enlarged  and  a  well-tested  hand  valve  incorporated  to  permit  operation  at 
the  required  pressure  level.  A  redesigned  combustion  chamber  with  an 
internal  volume  of  20  cubic  inches  was  constructed  with  hand  valve  seats 
machined  into  the  wall  of  the  chamber.  Combustion  experiments  are  now 
being  initiated  in  this  new  test  cell. 
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D.  Advanced  Instrumentation 
1.  General 

Although  any  one  of  a  large  number  of  tranaduction  methods  could 
be  adapted  to  shock  tunnel  testing,  initial  consideration  was  restricted  to 
those  classed  as  self-generating;  i.  e.,  no  external  source  of  electric  power 
is  required.  This  approach  was  chosen  to  simplify  instrumentation  and  to 
eliminate  electronic  development  to  accompany  the  transducer  development. 
Of  the  self-generating  systems,  the  piezoelectric  ceramic  transducer  ap¬ 
peared  to  be  the  most  attractive  for  a  reasonably  rapid  developmental 
program, 

A  variety  of  piezoelectric  transducers  are  commercially  available. 
Most  of  these  are  either  high-sensitivity  types  classed  as  microphones  or 
low-senaitivlty  devices  Intended  for  high-pressure  instrumentation.  Micro¬ 
phones  are  generally  too  large  for  model  instrumentation  and  not  suited  to 
shock  tunnel  testing  for  a  number  of  reasons.  Some  of  the  high-pressure 
transducers  would  be  suitable  in  many  respects  except  for  inadequate  sen¬ 
sitivity.  These  transducers  usually  have  a  circular  cylindrical  configura¬ 
tion  which  is  easy  to  install  in  models.  This  shape  is  a  logical  consequence 
of  the  type  of  piexoeloctric  elements  usually  employed;  hollow  or  solid 
cylinders  operating  In  compression.  Such  elements  provide  a  mechanically 
■tiff  system  with  high  resonance  frequency  which  limits  the  ultimate  sensi¬ 
tivity  attainable  with  elements  of  practical  dimensions.  The  higher  sensitiv¬ 
ity  is  achieved  in  microphones  in  part  by  the  larger  size  and  to  a  greater 
degree  by  the  type  of  pifezoolectric  element  employed.  ThiB  Ib  usually  a 
composite  bar  or  plate  with  the  output  resulting  from  flexure  or  torsion  of 
the  element.  Only  flexural  (bonder)  elements  can  he  conveniently  made 
with  ceramics.  The  resonance  frequency  1b  usually  much  lower  >(than  that 
of  the  cylindrical  types,  which  limits  the  frequency  response  of  tho  trans¬ 
ducer.  Thus  the  principle  objective  of  the  design  is  to  achieve  the  desired 
sensitivity  with  a  sufficiently  high  resonance  frequency. 
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2.  Bender  Element  Transducers 


A  functional  drawing  of  the  conventional  bender  element  transducer 
is  shown  in  Fig. 71.  A  positive  external  pressure  causes  a  downward  deflec¬ 
tion  of  the  element.  This  produces  tensile  stresses  in  the  upper  portion  of 
the  element  and  compressive  stresses  in  the  lower  portion.  The  element 
consists  of  a  "sandwich"  of  two  ceramic  strips  and  a  shim  in  the  neutral 
plane  of  the  composite  bar.  For  the  most  common  applications,  polarisa¬ 
tion  is  such  that  conducting  electrodes  on  the  upper  and  lower  surfaces 
provide  the  proper  output.  Analysis  of  such  a  structure  is  usually  made 
by  making  use  of  the  electrical  analog  of  the  mechanical  vibrating  system, 
including  the  mechanical-to- electrical  transduction  characteristics.  ,  1  ^ 
Figure  72  shows  the  equivalent  circuit  of  the  arrangement  shown  in  Fig.  VI. 
A  simplified  equivalent  circuit  ior  frequencies  well  below  the  system  reso¬ 
nance  is  shown  in  Fig  73.  A  cantilovei  element  is  shown  m  Fig.  7  1  although 
other  crystal  support  and  drive  systems  are  possible  without  altering  the 
configuration  of  the  equivalent  circuit.  The  values,  of  the  components  in 
thn  equivalent  circuit  arc-  dolurmuied  by  the  type  oi  piezoelectric  material, 
the  dimensions  of  the  element,  and  the  material  properties  and  dimensions 
of  the  diaphragm.  If  we  assume  for  the  moment  that  the  transdm  er  output 
is  essentially  open  circuited,  the  equivalent  uncut  will  yield  a  solution  for 
the  ratio  of  E  to  P  (volts  pur  unit,  pressure). 
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Ibis  represents  the  sensitivity  for  frequencies  well  below  the  sys 
tern  resonance  frequency  and  is  independent  of  frequency.  The  system 
resonance  frequency  is 


w 


(M.D  +  MC}  C 


D  +  CC 


where  w  is  the  angular  resonance  frequency  and  equal  to  2irf 
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Figure  71,  Functional  Arrangement  of  Pressure  Transducer  Using 
the  Piezoelectric  Bender  Element 
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Figure  72.  Equivalent  Circuit  for  Figure  7  1 
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Figure  73,  Low  Frequency  Equivalent  Circuit  for  Figure  71, 
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Transducer  with  AnuiI  Elements 


The  relationships  given  have  assumed  negligible  loading  by  the 
external  electrical  circuit.  This  will  be  approximated  if  the  internal  elec¬ 
trical  impedance  is  email  compared  to  the  load  impedance.  It  ia  theoretic, 
ally  poaeible  to  achieve  any  sensitivity  desired  with  a  given  resonance  by 
properly  choosing  the  crystal  and  diaphragm  constants.  However,  a  email 
structure  with  very  high  aonsitlvity  and  high  resonance  will  result  in  an 
unreasonably  high  source  impedance,  and  the  output  will  he  subject  to  sub- 
stantial  loading  effects.  In  addition,  there  are  practical  limits  to  the  ability 
to  fabricate  very  small  ceramic  elements,  particularly  with  respect  to 
thickness. 

3.  The  Axial  Element  Transducer 

Referring  to  Fig.  71,  It  is  seen  that  the  largest  crystal  dimension 
will  be  of  the  order  of  the  diaphragm  radius  it  the  body  diameter  Is  approx. 
Innately  the  same  as  the  diaphrugm  diameter,  This  putB  a  severe  restric¬ 
tion  on  element  length  which  can  limit  sensitivity  and  source  capacitance 
This  '■an  be  overcome  by  rotating  the  element  90°  so  th&L  its  longest  dimer. 
i»lou  la  parallel  to  the  axis  of  the  cylinder.  In  order  to  properly  drive  an 
element  so  positioned  a  pivoted  drive  arm  is  required.  Finally.  In  order 
to  balance  the  reaction  back  on  the  diaphragm,  two  eloments  ate  used  In  the 
configuration  shown  In  Fig.  74.  For  a  downward  motion  of  the  diaphragm, 
both  elements  are  deflected  outward.  The  elements  are  »  onncclcd  In  series 
and  phased  So  that  this  typo  of  stressing  gives  added  voltages. 

It  tan  he  scon  that  the  lateral  force  available  for  deflecting  the 
elements  Is  noi  (he  same  as  the  force  transmitted  by  the  diaphragm,  but  Is 
some  function  ol  the  drive  angle.  Figure  75  shows  the  relationship  between 
the  force*  and  displacements  lor  the  angled  drive.  If  we  define  mm  lumbal 
compliance  as 
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for  each  of  the  two  elements 


Cy 


AY 

AFy 

2 


t-y 


AX 

tan  a 

2AFx  tan  0 
- 


C.x 

7 

tan1"  0 


It  can  be  shown  that  the  angled  drive  system  is  equivalent  to  a 
mechanical  transformer  with  a  "turns  ratio"  of  tan  0/1.  This  satisfies  the 
conditions  derived  for  force  displacement  and  impedance.  The  equivalent 
circuit  for  this  configuration  is  given  in  Fig,  76  It  is  seen  that  for  0  =  45° 
the  turns  ratio  becomes  unity  and  the  performance  becomes  identical  to  that 
of  Fig.  7l  except  for  the  fact  that  the  single  bender  element  has  been  re¬ 
placed  by  two  with  the  same  aggregate  dimensions  and  mechanical  charac¬ 
teristics.  This  cunfigui  ation  can  always  be  employed  with  a  resultant 
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reduc  ed  to  1/4  of  the  single  element  value  it  is  not  ofttn  used 


As  with  the  conventional  structure,  we  can  derive  the  low  fre¬ 
quency  sensitivity.  This  is  found  to  be 
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The  factor  "2"  appears  in  the  denominator  because  C  refers  to 

C 

each  element  of  the  two  element  system  It  can  be  seen  tliai  for  a  given  set 
01  diaphragm  and  crystal  constants  there  is  an  optimum  angle  which  gives 
maximum  output.  This  relationship  is 
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The  increase  In  sensitivity  which  can  ba  achieved  by  maldng  $ 
optimum  in  comparison  with  <p  *  45°  is  shown  in  Fig.  It,  ft  can  be  seen 
that  great  increases  are  possible  only  whore  there  Is  a  largo  disparity 
between  the  crystal  and  the  diaphragm  compliances.  At  the  same  time, 
optimum  sensitivity  might  not  always  be  doaireii  since  lower  sensitivity 
transducers  would  be  necessary  for  measurement  of  pressures  close  to 
stagnation.  The  drive  arm  angle  could  bo  a  convenient  means  for  sensi¬ 
tivity  control. 

In  addition  to  sensitivity  and  shape,  the  factor  of  acceleration 
sensitivity  is  of  consequence  when  employing  pressure  transducers.  The 
important  parameter  here  lit  the  mass  of  the  transducer  moving  system. 
Referring  to  Fig.  74,  It  can  be  seen  that  axial  acceleration  acts  on  the  dia¬ 
phragm  mass,  but  not  on  the  elements.  Acceleration  in  the  flexural  plane 
of  the  elements  produces  opposing  voltages  in  each  of  the  elements  so  that 
a  degree  of  nrceloiatiou  compensation  is  provided  in  this  structure.  • 

Drawings  for  the  initial  version  of  the  transducer  have  been  sub¬ 
mitted  to  the  Laboratory  Shop.  Its  principle  design  characteristics  arc  ns 
follows: 

Body  Diameter 
Body  Length 

(not  including  terminals) 

Sensitivity 

or 

Internal  Capacitance 
Principle  Resonance 
Diaphragm  Active  Diameter 
Optimum  Drive  Angle 
Element  Dimensions 

The  above  figures  are  based  on  an  .tpproxltiwd  t«  (.(halation  of  Ihe 
diaphragm  characteristics  for  which  an  exact  solution  in  not  possible,  and 
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are  subject  to  modification  when  expert mental  dat*  availably.  The 

design  internal  capacitance  is  sufficient  for  developmental  purposes,  but 
too  low  for  operational  transducers.  Final  units  will  use  modified  ceramic 
elements  which  will  give  four  t  imes  the  internal  capacitance  but  only  half 
the  output  voltage, 

Complete  analysis  requires  taking  into  account  a  number  of  addi¬ 
tional  factors.  Some  of  these  may  be  easily  assessed  without  adding  to  the 
equivalent  circuit  and  the  over-all  analysis.  These  include  the  mass  of  the 
drive  arm,  the  acoustic  stiffness  of  the  internal  cavity,  and  the  acoustic 
leak  required  to  equalize  for  ambient  pressure  changes.  Other  factors  may 
require  some  experimental  evidence  in  order  to  determine  their  importance 
Included  in  this  class  is  the  effectiveness  of  the  drive  arm  in  transmitting 
the  diaphragm  forces  and  unanticipated  diaphragm  modes.  Coupling  of  the 
transducer  to  tho  medium  is  usually  made  through  a  small  diameter  hole 
leading  to  a  cavity  in  front  of  the  diaphragm.  This  can  constitute  a 
Helmholtz  resonator  with  substantial  effect  on  the  transient  response. 

This  can  bent  be  avoided  by  insuring  that  the  resonator  frequency  is  well 
ubovo  the  highest  frequency  of  interest.  Future  consideration  will  be  given 
to  acoustical  or  mechanical  damping  in  order  to  improve  the  transient  re  ¬ 
sponse  of  the  transducer. 


IV.  FUTURE  PROGRAM 


During  the  next  program  period,  the  following  investigations  will  be 
undertaken: 

A.  Theoretical  Investigations 

1.  Inviscid  three-dimensional  hypersonic  flows.  Analysis  of  lnviacid 
subsonic  and  transonic  regions  in  general  three-dimensional  flows.  In  par¬ 
ticular,  the  flow  fields  about  zero  angle  of  attack,  non- axisymmotric,  and 
anglo-of-attack  axiaymmetric  blunt  bodies  in  hypersonic  flow.  Tha  analyses 
to  be  initiated  may  be  divided  into  three  separate  general  methods: 

a.  An  inverse  method  of  integration  from  an  assumed  three- 
dimensional  shock  chosen  as  a  constant  value  in  a  spheroidal  coordinate 
system. 

b.  Direct  Dorodnitsyn-Betolacrkuvskii  procedure  for  a  three- 
dimensional  body  of  a  constant  value  spheroidal  coordinate  surface, 

c.  Perturbation  analysis  along  basic  flow  streamlines  Pertur¬ 
bation  equations  will  be  solved  in  senes  form  with  coefficients  depending 
only  on  the  nxlsymmotric  or  non -axisymmetric  basic  flow  and  appropriate: 
boundary  conditions . 

/.  Ktnetu  thooiy  apprna*  b  for  a  sliaip  leading- edge  flat  plate  in 
hypersonic  flow,  A  four -stream  Maxwellian  distribution  is  employed  in¬ 
stead  of  Lees'  and  Motl-Smith's  two-stream  Maxwellian  which  is  inadequate 
for  the  flat  plate  problem. 

3.  Numerical  solution  of  the  non- equilibrium  blunt  body  equations 
(11.13)  for  selected  flight  conditions  and  vehicle  geometries. 

4  Theoretical  predictions  of  the  non- equilibrium  flow  about  selected 
geometries  in  the  presence  of  the  frozen  ambient  gas  streams  associated 
with  the  combustion-driven  shock  tunnel. 

5.  Prediction  of  shock  layer  characteristics  downstream  of  the  stag¬ 
nation  region  employing  the  Navier-Stokes  model. 
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6  Analysis  of  shock-boundary  layer  interaction  associated  with 
control  surface  deflection  in  the  presence  of  thick  boundary  layers.  Pre¬ 
diction  of  upstream  separation  and  propagation  as  a  function  of  shock 
strength,  Mach  number,  Reynolds  number,  and  leading  edge  bluntness. 

B.  Experimental  Investigations 

1.  Installation  and  assembly  of  combustion  driver  and  six-foot  nosxle, 

2.  Calibration  of  combustion-driven  shock  tunnel 

a.  Probe  static  pressure  measurements 

b.  Langmuir  probes 

c.  Total  head  and  test- wedge  moiUuJ  etiiunU 

3.  Incorporation  of  electron  beam  densitometer  into  combustion- 
driven  Bhock  tunnel. 


4.  Control  suifaco  oil  ectiv  onus*  at  high  altitudes.  Vary  ratio  of  local 
boundary  flap  dimension  between  10*^  to  I.  Vary  leading  -  edge  bluntuesu 
no  an  to  clrtA,oo»«  influence  blast  wove  phenomena.  Correlate  with 
theoretical  Investigation  (II.  A.  6). 


5.  PrcHSUin  distribution  and  heat  transfer  about  blunt  and  nhaip 
leading-edge  flat  plates  In  older  to  experimentally  detei  mine  ialluein  e  ot 
non- equilibrium  flow  phenomena.  Langmuir  ptube  measurements  in  the 
shock  layer.  Correlate  with  theoretical  investigations  (It.  A.  land  If.  A.  4). 

6.  Langmuir  probe:  base  and  wake  electron  coni  miration 
measurements . 


7.  Pressure  and  heat  transfer  measurements  on  a  slurp  leading  edge 
(Kii  >  1}  flat  plate  in  the  range  5  <  ^  <  1000.  Correlate  with  theoielual 
investigation!)  (II.  A.  2  and  II.  A.  5) 
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technique  establishes  the  partial  derivatives  along  this  new  boundary.  In 
turn,  the  partial  derivatives  normal  to  this  boundary  are  determined  from 
the  governing  flow  equations.  Successive  applications  of  this  scheme  al¬ 
lows  the  computation  to  be  extended  to  the  unknown  body  surface  which  will 
be  established  as  a  part  of  the  solutions  based  upon  the  consideration  of 
flow  continuity. 

The  equations  of  motion,  in  terms  of  the  shock-oriented  curvllln- 


ear  coordinates  may  be  written  (Fig.  24), 
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The  equation  of  state  is  written  as 

p  =  (R /)?{)  pT 

(65) 

Where  R  is  the  universal  gas  constant,  and  /W  is  the  molecular  weight, 


n  denotes  the  total  number  of  chemical  components  in  the  system.  The 
above  system  of  (n  plus  5)  equations  will  be  used  for  the  determination  of 
on  ,  T,  p,  p,  u,  and  v,  which  completely  describe  the  flow  field. 
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In  Eq.  (60)  the  mass  rate  of  production  of  the  component  by 
chemical  reactions,  ,  is  given  by 
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Mi  ^i 

A  dt 


(66) 


where  A  is  the  Avogadro's  number,  v.  is  the  number  of  particles  of  the 

t  3 

i  component  per  unit  volume,  (particle/crn  ). 

In  general,  chemical  reactions  can  proceed  in  both  forward  and 
backward  directions.  A  set  of  opposing  chemical  reactions  may  be  repre¬ 
sented  by  (Ref.  31), 
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where  and  aj'  are  stoichiometric  constants. 

In  these  representations,  any  "third  body”  chemical  species,  p.  ^ , 
which  acts  as  an  energy  source  or  sink  in  a  chemical  reaction  but  is  other¬ 
wise  unchanged,  is  included. 

The  net  rate  of  production  of  ith  species  is  then  given  by 
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At  thermodynamic  equilibrium  there  is  no  net  change  In  composition, 
therefore,  the  rate  constants  k.  and  k,  must  be  related  through  the  equi- 
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librium  constant,  K,  That  is, 
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When  multiple  chemical  reactions  .are  involved,  the  net  rate  production  of 

ith  species  can  be  obtained  by  summing  Eq.  (68)  over  all  the  chemical  re- 

actions  in  the  system.  Where  K  is  the  equilibrium  constant  in  the  unit  of 

(particle/cm  )  while  the  rate  constant,  k,  has  the  unit  of  [cm^/(sec  par- 
2 

tide  )] .  In  the  case  of  atmospheric  air  the  dominant  chemical  reactions 
with  their  respective  rate  constants  are  listed  in  Table  1. 

The  rate  constants  for  all  the  chemical  reactions  with  various 

32 

"third  bodies"  used  in  this  analysis  are  those  given  by  AVCO.  It  should 
be  understood  that  these  numerical  values  of  the  reaction  constants  are  to 
be  considered  as  recommended  estimates  only.  Modifications  can  readily 
be  introduced  when  better  experimental  data  become  available.  Related 
discussion  on  chemical  kinetics  may  be  found  in  Ref.  .13. 

Combining  Eqs.  (60),  (66),  (68),  and  (69)  on  the  basis  of  chemical 
ronrlmns  and  their  rate  constants  given  in  Table  1,  the  following  set  of 
five  equations  may  be  obtained  in  place  of  Eq.  (60): 
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In  obtaining  the  above  equations,  proper  care  has  been  taken  to 
convert  the  rate  constants  kf ,  1^,  in  the  unit  of  (cm  /particle  sec)  to  , 
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,  k|  in  the  units  given  in  Table  1.  Furthermore,  the  following 
non-dimensional  quantities  have  been  introduced: 
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where  L,  py,  py  ,  T,  r'^r  ,  uy  are  chosen  reference  quantities 

ur  =  ^  *  Pr  =  Pr  f^)  Tr  ’ 

In  terms  of  these  same  non-dimensional  quantities,  Eq.  (62)  may  be 
rewritten  as: 
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Differentiating  the  equation  of  state,  Eq.  (65),  with  respect  to  y  ,  we  have, 
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Dividing  the  partial  differential  form  of  the  energy  equation,  Eq.  (64),  by 
the  absolute  tempgrature,  we  obtain, 

8a, 


ill  Z 


T  0y  i 

Substitution  of  (81)  into  (82)  yields, 
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Eliminating  (8p/8y)  and  (8p/8y)  from  Eq.  (83)  by  using  Eq.  (63)  and  Eq.  (61) 
respectively**  we  obtain 
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The  non-dimensional  form  of  Eq.  (84)  reads, 
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The  equation  of  state  may  be  rewritten  as, 
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^0_  +  _^N  +  °NO  +  °2  + 


-1 


\%  *NO  %z  *1 


PT 


Nz/ 


(86) 


Now  the  reacting  flow  field  inside  the  shock  layer  is  governed  by  the  differ¬ 
ential  system  consisting  of  Eqs.  (70,  71,  72,  73,  74,  78,  79,  80,  85,  81, 
or  86).  There  are  ten  equations  for  the  ten  unknown  quantities 


aO’  aN  ’  a02’  QN2’  QNO’  u  -  v>  P-  P-  t 
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Where  0  denotes  any  of  the  dependent  variables  0  is  the  value 
of  0  at  the  mth  step.  The  values  of  0q,  (80/9x)q  are  initially  obtained 
from  the  oblique  shock  calculations.  The  computations  are  carried  out  by 
applying  this  scheme  successively  until  the  body  streamline  is  established 
from  continuity  consideration. 

Consider  a  streamline  crossing  the  curved  shock  at  an  arbitrary 
point,  B  (Fig.  25).  The  amount  of  fluid  enclosed  between  this  streamline 
AB  and  the  streamline  along  the  axis  of  symmetry  must  be  equal  to  the 
amount  that  crosses  the  section  BC  inside  the  shock  layer.  It  should  be 
understood  that  for  the  axisymmetric  flow,  the  control  surfaces  will  be 
obtained  by  revolving  AB  and  BC  about  the  axis  of  symmetry.  The  above 
mentioned  physical  condition  may  be  represented  by  the  following  equation: 

yw 

n  r8  Poo  Uod  =  Zn  f  T  pu  dy  •  (87) 

o 
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SYMBOLS  (Section  II.  B) 


Avogadro  number 
specific  enthalpy 

h/h  *  where  h  =  (p  Ip) 

T  w  r  'roo  rao' 

reaction  rate  constant,  cm^/[sec  (particles)^] 
dissociation  rate  constant,  [cm  /{mole  sec)],  see  Table  I 

recombination  rate  constant,  [cm^/(mole*  sec)],  see  Table  I 

3 

reaction  rate  constant,  [cm  /(mole  sec)],  see  Table  1 


curvature 


KL  =  Kno 

Equilibrium  constant,  particles /cm^ 

2  • 

Equilibrium  constant,  s  (pQ  / pQ  ) 

2 

2 

Equilibrium  constant,  5  (p^/p^  ) 

2 

Equilibrium  constant,  =  [  (pNO  Pq^pO 


Equilibrium  constant,  =  [  (pN  Pq)/(pn0 

2 


pN)] 

pN)i 


reference  length,  = 

Mach  number 
molecular  weight 
pressure 

!p/pr)  -  where  pr  =  p^ 

radial  distance  from  the  axis  of  symmetry 


SYMBOLS  (Continued) 


a 


R 

t 

T 

T 


u 

u 


V 


00 


V 

V 

X 

X 


y 

y 

Q, 

\ 

V 

P 

P 

a 

<P 

n 

n 


Universal  gas  constant 

time  ,  sec 

temperature,  °K 

=  (T/T  )  ,  where  T  =  T 

r  r  oo 

velocity  component  along  the  x  direction 

=  u/u  ,  where  u  =  *Jp  To 

r  r  roo 

free  stream  velocity 

velocity  component  along  the  y-direction 

=  v/u  ,  where  u  =  n/»  7p 

r  r  1  oo  v<x> 

coordinate  distance  along  the  shock  surface 
=  x/no 

coordinate  distance  normal  to  the  shock  surface 
=  y/no 

=  9-J  P 

=  specific  heats  ratio  for  perfect  gas 

3 

=  number  density  (particle/cm  ) 

3 

=  number  density  of  "third  body"  species  (particle/cm  ) 

=  density 
=  shock  angle 
=  arbitrary  function 
=  radius  of  curvature 

=  radius  of  curvature  of  shock  wave  at  stagnation  streamline 
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SYMBOLS  (Continued) 


Subscripts: 

i  ith  species  or  arbitrary  index  on  rate  constants 

j  jth  species  or  arbitrary  index  on  rate  constants 

0  quantities  at  x  =  c  or  sea  level  standard  conditions 

r  reference 

s  conditions  immediately  behind  the  shock 

w  wall  conditions 

oo  free  stream 


c.  Hypersonic  Wake  Flow  Studies 


1.  Introduction 

As  a  body,  such  as  a  nose  cone,  enters  the  earth's  atmosphere, 
it  leaves  a  deposit  of  energy  in  the  form  of  a  trail.  The  rate  at  which 
energy  is  deposited  is  equal  to  the  drag  energy  DU  where  D  is  the  aerody¬ 
namic  drag  and  U  is  the  flight  velocity.  Part  of  this  energy  goes  into 
heating  the  air  and  part  goes  into  giving  the  air  directed  motion  or  momen¬ 
tum  in  the  flight  direction.  It  is  the  heated  air  that  gives  rise  to  observ¬ 
ables  in  the  form  of  an  ionized  trail  and  gas  radiation. 

After  the  air  has  passed  through  the  bow  shock,  around  the  body, 
and  through  the  trailing  shock  system,  as  shown  in  Fig.  26,  it  returns  to 
ambient  pressure  at  some  distance  xq  behind  the  body.  At  this  station  a 
certain  velocity  and  temperature  profile  would  persist  downstream  if  it 
were  not  for  laminar  diffusion,  turbulent  mixing,  and  radiation. 

The  amount  of  energy  radiated  from  the  wake  is  extremely  small 
compared  to  the  other  controlling  processes  of  energy  transfer,  and,  hence, 
can  be  neglected  in  any  wake  analysis. 

Figures  26a  and  26b  are  schematics  of  what  is  qualitatively  ex¬ 
pected  in  the  trail  behind  a  blunt-nosed  body  and  a  slender -pointed  body, 
respectively.  Both  bodies  are  assumed  to  have  about  the  same  base  area 
and  volume.  The  skin  friction  or  viscous  drags  of  the  two  bodies  will  be 
roughly  the  same  at  the  same  velocity  and  altitude.  This  fact  is  reflected 
in  the  temperature  and  velocity  profiles  shown  in  Figures  26a  and  26b,  and 
is  represented  by  the  small  bumps  at  the  center  of  the  profiles.  This  is  a 
schematic  representation  only.  In  reality  the  momentum  defect  associated 
with  the  bump  should  include  not  only  skin  friction  drag,  but  nose  drag  in 
the  case  of  a  sharp  body  with  a  small  nose  radii,  and  the  drag  associated 
with  the  base  recompression  at  the  neck.  This  bump  in  the  velocity  profile 
constitutes  the  inner  core  and  its  spread  will  be  dominated  by  viscous  pro¬ 
cesses,  either  laminar  or  turbulent,  if  transition  has  occurred. 
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VISCOUS 


Figure  26a.  ScVUmatic  of  Flow  Past  a  Blunt  Body 


Figure  26b.  Schematic  of  Flow  Past  a  Sharp  Body 


The  associated  bump  in  the  temperature  or  enthalpy  profile  will 
reflect,  in  addition,  the  heat  taken  out  of  the  flow  by  the  body,  and  its  peak 
will,  therefore,  be  much  lower  than  in  the  case  of  an  insulated  body.  As 
will  become  apparent  in  the  analysis,  heat  transfer  to  the  body  plays  a 
dominant  role  in  erasing  observables  in  the  wake  of  a  sharp  body  which  are 
dependent  on  the  maintenance  of  a  high-temperature  environment.  Although 
the  heat  transfer  to  the  blunt  body  is  of  the  same  order  as  the  heat  transfer 
to  the  sharp  body  at  the  same  velocity  and  altitude,  a  high-temperature 
environment  ^s  maintained  by  the  high  temperature  air  that  has  passed 
through  the  steep  portion  of  the  bow  shock  and  surrounds  the  viscous  core 
as  shown  in  Fig.  26a,  so  that  in  the  blunt-body  case  heat  transfer  to  the 
body  does  not  play  an  important  role  in  erasing  observables. 

2.  The  Spread  of  Turbulence  in  Axisymmetric.  Wakes 

%hia  section  presents  a  simplified  analysis  of  the  spread  of  tur¬ 
bulence  %i  the  trail  wherein  certain  simplifying  assumptions  are  made. 

Thes«  assumptions  may  be  summarized: 

a.  Equilibrium  flow, 

b.  Negligible  radiation, 

c.  Negligible  laminar  diffusion  compared  to  turbulent  mixing, 
and 

d.  Turbulence  originates  in  the  viscous  core  or  on  the  body  and 

is  fully  developed  and  in  equilibrium  by  the  time  it  has  reached 
the  50  radii  station. 

Turbulence  is  thus  confined  to  a  circular  core,  the  front  of  which 
moves  out  and  "swallows  up"  mass  from  the  outer  entropy  trail.  Eventually 
the  front  will  have  moved  out  into  the  outer  cool  air  at  the  edge  of  the  en¬ 
tropy  trail,  and  the  entire  trail  begins  to  cool  by  turbulent  mixing.  As  to 
how  long  this  process  takes  (or  how  far  downstream)  before  the  entropy 
trail  produced  by  the  normal  portion  of  the  shock  of  a  blunt  body  is  cooled 
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sufficiently  to  erase  any  observables  such  as  electron  density,  will  depend 
upon  certain  initial  parameters  at  the  station  xq  .  These  parameters  are: 

a.  Total  drag  coefficient, 

b.  Ratio  of  drag  associated  with  the  turbulent  inner  core  to  the 
total  drag, 

c.  Width  of  the  entropy  trail,  and 

d.  Width  of  turbulent  inner  core. 

3.  Basic  Equations 

We  will  use  the  stream  function  41  and  x  as  independent  variables 
in  this  analysis.  The  stream  function  is  related  to  the  physical  radial  co¬ 
ordinate  r  through  the  formula 

r 

=  2  tt  J  p  u  i  dr  (88) 

o 

and  r  can  be  found  from  the  inverse  relation 


2 

r 


(89) 


Thus  the  radius  R  of  the  stream  tube  4*  in  the  free  stream  up  to  the  bow 
shock  is 


R2  = 


TT  p  U 
rCO  00 


(90) 


Let  4^  be  the  stream  tube  which  enters  the  effective  front  of  the 

inner  turbulent  core  at  x  and  let  4**  he  the  stream  tube  that  enters  at  x  , 

f  o 

the  origin  of  calculation.  Then  from  conservation  of  momentum  in  the 
x-direction  we  have 


I 


o 


d  4<  + 


(91) 
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since  there  is  no  x-momentum  transported  through  stream  tube  boundaries 
in  the  outer  inviscid  flow  and  the  pressure  is  assumed  to  be  constant. 

A  second  relation  is  founu,  the  integral  relation  between  the  kinetic 
energy  and  the  dissipation  in  the  turbulent  core,  viz. 


a  f{  1  2  .  .  1  2  d  ^  .  2  ff  2 

d x  J  2Ut  d*  =  2  Uf  “dT  ‘  4ir  J  £tr  ptUt 


2  u, 

T 


(92) 


By  substituting  the  profile 


.  dj 

u  +  (u,  -  U  )  -f- 

a  f  a  4v 


(93) 


into  the  integral  relations  (91)  and  (92)  one  obtains 


jj'f 

i(uf  +  ua)^f  =  y(uf  +ua)+f  +[  u.  d^ 

o  o  o  J 


(94) 


d  u 

F(uf  "V-ar^f  =  4Tp2/  £tr2ptut(uf-ua)2^2dv|;  (95) 

o  ^f 

Where  both  Eq.  (94)  and  its  differential  form  were  used  in  reducing  the  left 
hand  side  in  Eq.  (92)  to  that  in  Eq.  (95).  4^,  ua  -  and  uf  are  assumed  to  be 
purely  functions  of  x  with  initial  values  of  >  ua  >  a”d  uf  at  x  =  xq  . 
Equation  (95)  can  then  be  further  reduced  to 


d  u 

a 

dx' 


24  ir 


u,  -  u 
2  f  a 


/'■ 


pt  utd  + 


(96) 


In  order  to  obtain  a  simple  evaluation  of  the  dissipation  integral  on 
the  right  hand  side  of  Eq.  (96)  we  restrict  ourselves  to  the  case  where  the 
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variations  of  velocity  and  density  are  small  across  the  core.  For  this  case 
we  make  the  "incompressible"  assumption 


Pt  Ut  "  Pf  Uf 


(97) 


and  approximate  «  by 


£h  =  KPf  rf|uf  *  ua|  =  — 


n/tt 


K  \  l^{  Pf 


u,  -  u 


1  f  a 1 


(98) 


Substituting  these  into  (96)  we  obtain 
du 


— s— ~  *  1 2  K  \  /■!  —  —  I  u  c  -  u  I  (u  ,  -  u  ) 
dx  f  a  1  a 


(99) 


Now  from  the  relation  (94)  we  obtain  u  in  terms  of  u,  and  an 

a  i 


expression  for  du  /dx; 


u 


r 

1  n‘ 

•h 


(uf  +»„>+, 


+ 


o  o 


-  Ur 


(100) 


d  u 

a 

dx 


2uf  duf  ? 

vl/f  d  f  a  i  ,  Z 

.  f  f  °  °  o  -q, 


1  u‘ 


di)> 


d  d  ua  d  ^f 

dx  ~  d  iK  dx 


(101) 


Since  u.  is  assumed  to  be  known  as  a  function  of  v|; ,  it  is  a  simple  matter 
to  obtain  x  as  a  function  of  ^  by  numerical  integration  of  Eq.  (99). 

We  have 


f 

Tu 
+f 

O 


d+f 

7'  Ua  '  (uf  -  "a1 


Vf 


li[if  =  12K  n/w  (x  -  x  ) 


(102) 
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From  the  graph  of  x  vs  4^  we  can  obtain  4*£  vs  x  and,  hence,  the  flow 
quantities  and  equilibrium  electron  densities  as  a  function  of  x  and  4J  and 
ultimately  as  a  function  of  x  and  r  through  the  relation  (89). 

This  formulation  is  for  an  arbitrary  u.  as  a  function  of  i|j  .  We 
shall  look  at  the  special  case  where  m  can  be  approximated  by 

u.  (i|i)  =  UflD(l-  ae’^)  (103) 


.This  distribution,  by  the  proper  choice  of  a  and  b,  can  be  made 

to  fit  the  velocity  defect  produced  by  a  blunt  body  fairly  well.  It  gives  the 

bell-shaped  nearly  Gaussian  curve  for  the  outer  entropy  trail  shown  in  Fig.  26a 

2 

when  plotted  in  the  physical  coordinate  r  .  (Note  that  ^  ^  r  from  Eq.  (88.) 

At  the  initial  station  x  =  xq  let  the  drag  associated  with  the  inner 

turbulent  core  be  D  and  let  D  be  the  total  drag.  We  then  have 
w 

^f 

Dw  *  [  <Uoo-“t)d+  *  -ua»+£  11041 

J  o  o  o 


CO 

D-D  =  (U  -  u.  )  4 1,  +  /  (U  -  u.)  dip 

w  oo  f  Tf  J  oo  l 


o  o 


(U 


oo 


+ 


£ 

b 


(U 


00 


from  which  one  readily  obtains  the  relation 


(105) 


(106) 
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U 


D  = 


D  = 


w 


Q  = 
Q  = 


We  shall  introduce  at  this  point  the  dimensionless  quantities: 


u/u 

CO 

(107) 

p/p 
r  roo 

(108) 

t/Tqo  =  P^/P  =  1/p*  (since  p  =  const.  =  Pqo) 

(109) 

D/l  p  U  2  A 

2  ^oo  oo 

(110) 

dw4  <>«.  a 

(111) 

w 

(112) 

W+b 

(113) 

!/b^b 

(114) 

b  4*f  =  4^*/+* 

(115) 

V* 

(116) 

flight  speed 

atmospheric  density  at  flight  altitude 
absolute  temperature  at  flight  altitude 
ambient  pressure  at  flight  altitude 

2 

re-entry  body  maximum  x- sectional  area  =  tt  r^ 
maximum  body  radius 

p  U  A  a  stream  function  value  of  free-stream  tube  equal  to 
roo  oo 

body  x- sectional  area 
total  drag 

drag  contained  in  turbulent  core  at  the  start  of  calculation 
total  heat  transfer  rate  to  the  body 

Q/_.  u  ik  ili  =  dimensionless  heat  transfer  rate  coefficient. 

2  oo  Yb  T 
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Equation  (101)  can  then  be  written  in  non-dimensional  form,  viz. 


J  p* — 

•L  U,  -  u 


o 


w .  .  v  v. 
a  1  <uf  ua 


>K  , 
q  Uf  Tf  dn 


12  K  i  X  ’  o 


(117) 


where 


1  -  -  C*  }  ' _k  e  °  e  ~n 
2  D  1  +  r\0 


j  Uf*  dn  +  \  2  - 


1  -  *  1  -  k  'o' 

IcDl-nro'  , 


(118) 


1  1  -  k  -  *  ‘o  ,  ,  2  -1 

+  2  ‘  T+f£  CD  e  1  +^J  e 


(119) 


#  #  CD  ,  1  -  k  no  M  , 

u,  -  u  =  -=—  1  -  T-jTrT  e  (1  +  9)  e 
f  a  1  +  n 


(120) 


tf  = 


(121) 


*}c  v  ••  1  2  ^ 

T£  =  1+VM0O  1  -”f 


(122) 


T*  =  l+lfiM2  1  -  u*  -  Q*  J 
a  2  oo  a  H 


(123) 


u,  -  u 
f  a 
o  o 


(124) 
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It  can  be  seen  immediately  that,  if  Eqs.  (120)  and  (121)  are 
substituted  into  Eq.  (117),  a  factor  1/Cn  can  be  pulled  out.  Thus  multi 
plying  Eq.  (117)  by  we  obtain 


F  (Q‘ ,  ,  V  ,  k  ,  no  ,  ;  n)  dn  = 


12  K  C_,  x  -  x 

D  -2  =  i  (125) 


with 


F  = 


n  *  * 


l  - 


1  -  k  -1 1  fi  1  -  k  V,_ 

T  +  rT6  (1  +  r|)  e  '|  1-pq^-e  (1  +  ^1) 

o  o 


-n 


/  * 

•V1!  Uf  Tf 


(126) 


4.  Discussion  and  Numerical  Examples 

Examination  of  Eq.  (125)  and  (126)  reveals  that  the  major  param¬ 
eters  that  control  the  spreading  of  the  turbulent  core  occur  only  in  the 
combinations 


*  1  -  k 

D  ’  1  +  t,q 


n 


o 


and  that  (x  -  x^/r^  can  be  expressed  entirely  in  terms  of  these  parameters. 
It  can  further  be  shown  that  these  parameters  can  all  be  expressed  in  terms 
of 

y  ,  Mffl  ,  K ,  CQ  ,  k ,  Tf  ,  u£  -  u& 

o  o  o 


Although  the  exponential  variation  of  velocity  given  by  Eq.  (103) 
was  chosen  on  the  basis  of  a  good  fit  for  a  blunt  body  with  a  spherical  nose, 
the  approximation  can  still  be  used  for  sharp  bodies,  such  as  a  cone,  with¬ 
out  loosing  any  of  the  salient  features  in  describing  the  turbulent  core 


spreading  and  temperature  on  the  axis  as  a  function  of  (x  -  x^)/r^.  This 
is  due  to  the  fact  that,  for  a  sharp  body,  the  radial  distribution  of  velocity 
and  temperature  in  the  inviscid  flow  is  fairly  flat  and  the  gradients  at  the 
edges  are  not  important.  This  is  demonstrated  in  the  two  numerical  cal¬ 
culations  discussed  in  the  ensuing  paragraphs. 

Two  numerical  calculations  at  150,  000  ft  and  Mach  20  are  pre¬ 
sented;  one  for  a  typical  blunt  body  such  as  that  shown  in  Fig.  26a  and  one 
for  a  typical  sharp  body  shown  in  Fig.  26b.  The  drag  contained  in  the  tur¬ 
bulent  core  bucket  is  assumed  to  be  the  same,  that  is,  C  =  .06.  This, 

,Jw 

in  both  cases,  includes  the  skin  friction  and  base  drag,  phis  nose  drag  in 

the  case  of  the  sharp  body.  At  150,  000  ft,  the  skin  friction  coefficient  is 

about  0.03,  and  using  Reynolds  analogy,  the  heat  transfer  coefficient  is 

0/4-  p  U  3  A  =  0.015  .  These  values  were  assumed  for  both  cases.  Both 
2.  co  oo 

calculations  were  started  with  the  assumption  that  the  velocity  on  the  axis 
was  zero  at  the  beginning  which  would  be  the  case  at  the  effective  rear 
stagnation  point.  Each  case  will  now  be  discussed  separately. 

a.  Blunt  Body 

A  typical  pressure  drag  coefficient  for  a  blunt  body  is  0.54, 
giving  a  total  drag  coefficient  of  Cn  =  0.6  with  k  =  0.1  .  The  results  are 
shown  in  Fig.  27  where  r| ,  T^.  and  l'a  with  and  without  heat  transfer  are 
plotted  against  (x  -  x^/r^  .  A  measure  of  how  far  the  front  has  moved 
laterally  in  terms  of  Gaussion  widths  of  the  entropy  trail  is  given  by  *1 . 

It  can  be  seen  that  in  the  blunt  body  case,  the  front  must  move  many 
Gaussian  widths  before  the  center  of  the  turbulent  core  begins  to  cool. 

This  distance  is  of  the  order  of  a  hundred  body  radii.  The  solid  curve  for 

>jc  t  t 

T  is  the  temperature  ratio  on  the  axis  for  an  insulated  body,  while  the 

cl 

dashed  curve  is  for  the  assumed  value  of  the  heat  transfer  rate.  It  can  be 

seen  that  the  heat  transfer  rate  to  the  body  has  little  effect  in  depressing 

#  >!< 

the  temperature  except  at  the  very  beginning  where,  for  u&  =  0,  T&  is 
about  0.56  times  the  free  stream  stagnation  temperature  (strictly  speaking 
T>!  is  the  enthalpy  ratio  and  not  the  temperature  ratio). 
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It  can  be  concluded  that  for  a  typical  blunt  body  with  a  two-foot 
diameter  the  temperatures  will  be  sufficiently  high  to  produce  large  amounts 
of  radiation  and  sustain  high  electron  concentrations  visable  to  radar.  Fur¬ 
thermore,  the  width  of  the  turbulent  core  becomes  very  large  for  the  blunt 
body  case,  as  evidenced  by  the  large  values  of  q  in  terms  of  body  radii. 

)jc  I'f 

When  T  has  reached  roughly  half  the  value  of  T,  ,  the  initial  front  tem- 
a  *o 

perature  in  the  entrophy  trail,  ti  has  the  value  5  corresponding  to 

r^/r^  =  18.1  as  calculated  from  Eq.  (89),  and  (x  -  XQ)/r^  is  117. 
b.  Sharp  Body 

A  typical  pressure  drag  coefficient  for  a  sharp  body,  say  a  10° 

half  angle  cone  or  an  equivalent  ogive,  is  0.06  giving  a  total  drag  coefficient 

of  CD  =  0.12  with  k  =  0.5.  The  results  are  shown  in  Fig.  28  where  the 

same  quantities  are  plotted  as  in  the  blunt  body  case.  It  can  be  seen  that 

the  wake  is  cooled  in  a  much  shorter  length  than  in  the  blunt  body  case. 

$ 

The  temperature  on  the  axis  has  reached  the  value  =  9  at  a  value  of 
(x  -  xQ)/r^  =  33,  with  =  2.1  for  a  value  of  r^/r^  =  5.4. 

It  can  be  concluded  that  for  a  typical  sharp  body  the  dimensions 
of  the  high'temperature  region  in  the  wake  are  about  a  factor  of  three  less 
than  for  the  blunt  body  case.  The  heat  transfer  rate  to  the  body  has  a  more 
significant  effect  on  the  reduction  of  these  dimensions  than  in  the  blunt  body 
case. 

The  values  of  the  parameters  chosen  for  the  sharp  body  are 
conservative  and  they  tend  to  make  the  wake  larger  than  that  which  might 
actually  be  expected  at  this  altitude. 


9 1? 


Figure  28.  Wake  Temperature  Decay  for  Slender  Body  in  Hypersonic  Flow 


III.  EXPERIMENTAL  INVESTIGATIONS 


A.  Stagnation  Point  Heat  Transfer  Measurements  in 
Hypersonic  Low  Reynolds  Number  Flows 

1.  Introduction 

In  determining  the  point  at  which  protective  fairings  for  delicate 
payloads  should  be  jettisoned  during  the  launch  trajectory  of  satellites  ve¬ 
hicles,  free  molecule  flow  heat  tr  ansfer  rates  have  been  assumed.  Calcu¬ 
lations  based  on  this  assumption  provide  the  maximum  heat  transfer  rate 
predictions,  and  hence,  the  maximum  payload  protective  fairing  require¬ 
ments.  As  a  result,  the  tendency  has  been  to  carry  the  fairings  to  altitudes 
which  are  much  beyond  optimum  because  of  the  uncertainty  of  the  heat 
transfer  rates  in  the  very  high  atmosphere.  These  rates  were  not  signifi¬ 
cant  relative  to  exit  and  re-entry  heat  rates  and  these  unknowns  were 

tolerable.  Now  situations  are  arising  where  it  is  no  longer  advisable  to 
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ignore  these  uncertainties.  Also  theoretical  analysis  in  the  low  density 
regime  is,  for  the  most  part,  unsubstantiated  by  experiment..  It  is,  there¬ 
fore,  the  purpose  of  this  study  to  attempt  to  measure  the  stagnation  point 
heat  transfer  rates  to  cylinders  and  spheres  in  the  hypersonic  shock  tunnel 
at  flow  Mach  numbers  from  13  to  23.  The  Reynolds  number  range  to  be 
covered  is  from  5  to  500,  based  on  leading  edge  diameter,  corresponding 
to  the  regime  between  free  molecule  flow  and  the  continuum  regime. 

Because  of  its  high  reservoir  stagnation  temperature,  the  shock 
tunnel  is  an  instrument  peculiarly  suited  to  making  measurements  at  high 
Mach  numbers  and  low  Reynolds  numbers  typical  of  the  high  altitude  por¬ 
tions  of  the  launch  trajectory  of  satellite  vehicles.  Furthermore,  the  shock 
tunnel  can  yield  a  variety  of  stagnation  pressures,  and  hence,  Reynolds 
numbers  for  a  given  stagnation  temperature.  The  chief  drawback  inherent 
in  shock  tunnels  is  the  comparatively  short  time  interval  during  which 
measurements  can  be  made.  As  a  result,  equilibrium  conditions  between 
the  flow  and  the  measuring  instrument  usually  cannot  be  obtained  and  must 
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be  inferred  from  the  transient  response  of  the  measuring  device.  In  this 
report  heat  transfer  rate  measurements  determined  from  the  transient 
response  of  thin  film  resistance  thermometers  mounted  at  the  stagnation 
point  of  models  will  be  described  and  the  results  of  the  measurements  pre-i 
sented.  The  shock  tunnel  facility  in  which  these  tests  were  made  is  shown 
in  Fig.  29. 

2.  Theoretical  Predictions 

34 

Probstein  and  Kemp  have  considered  the  problem  of  calculating 
the  aerodynamic  heating  to  be  expected  in  hypersonic  rarified  gas  flows, 
where  the  mean  free  path  becomes  too  large  for  the  use  of  the  inviscid 
flow-boundary  layer  approach  and  too  small  for  the  free  molecule  flow 
regime  to  apply.  They  have  solved  the  Navier  stokes  equations  for  two  sets 
of  boundary  conditions;  (1)  a  "viscous  layer"  class  in  which  a  thin  shock 
wave  exists,  but  the  region  between  the  shock  and  the  blunt  body  is  fully  .< 
viscous,  and  (2)  a  "merged  layer"  class  representing  a  more  rarified  con¬ 
dition.  In  the  merged  layer  the  shock  wave  is  no  longer  thin  and  the  Navier 
stokes  equations  are  used  to  provide  a  solution  which  includes  the  shock 
structure  and  free  stream  conditions  as  outer  boundary  conditions. 

From  this  analysis  a  deviation  is  predicted  in  the  stagnation  point 
heat  transfer  rate  from  that  expected  on  the  basis  of  extrapolated  boundary 
layer  theory.  The  viscous  layer  theory  predicts  higher  heat  transfer  to 
both  spheres  and  cylinders  /normal  to  the  flow  as  the  Reynolds  number 
(based  on  nose  diameter)  is  decreased  to  the  order  of  10. 

The  "incipient  merged  layer"  theory  modifies  the  heat  transfer 
rates  to  values  which  are  less  than  in  the  viscous  layer,  but  greater  than 
in  the  boundary  layer  theory.  The  predicted  results  for  spheres  and  cyl¬ 
inders  are  shown  in  Fig.  30a  and  30b.  Although  an  incipient  merged  layer 
prediction  for  the  cylinder  i6  not  made,  it  is  easy  to  suppose  that  it  lies 
between  the  viscous  layer  and  boundary  layer  value.  These  curves  were 
computed  on  the  basis  of  a  density  ratio  across  the  shock  of  0.1  and  a  y 
for  the  gas  of  11/9. 
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Figure  29.  View  of  Shock  Tunnel  Showing  12-inch  Test  Section 


igure  30a„  Heat  Transfer  Rate  to  Cylinders  Normal  to  the  Flow  as 
Function  of  Stagnation  Point  Reynolds  Number 


'igure  30b.  Heat  Transfer  Rate  to  Spheres  as  a  Function  of 
Stagnation  Point  Reynolds  Number 


For  shock  tunnet  reservoir  conditions  of  about  240  atmospheres 
and  2840°K,  and  test  section  Mach  numbers  from  10  to  20,  the  heat  trans- 
fer  rates,  according  to  Fay  and  Riddell,  were  calculated.  These  values 
are  plotted  in  Fig.  30  for  comparison  with  the  results  of  Probstein  and 
Kemp.  A  free  molecule  flow  value  for  an  accommodation  coefficient  of 
unity  is  plotted  also. 

3.  Experiments 

In  the  shock  tunnel  for  shock  Mach  numbers  of  4.5  to  4.7,  reser¬ 
voir  test  conditions  (which  remain  steady  for  about  4  millisec)  of  between 
220  and  240  atm  and  2840  and  3000°K  are  obtained.  Using  the  twelve-inch 
diameter  test  section  and  varying  the  throat  diameter  over  a  range  of  .04 
to  .375  inches,  Mach  numbers  from  10  to  23  are  achieved.  For  each  test, 
the  shock  Mach  number  (Mg)  in  the  shock  tube,  and  the  pitot  pressure  (Pt) 
in  the  test  section,  are  measured.  From  the  shock  Mach  number,  the 
reservoir  conditions  (that  is,  the  conditions  behind  the  reflected  shock) 
are  ascertained.  These  conditions,  along  with  the  pitot  pressure  measure¬ 
ment,  determine  the  flow  Mach  number  in  the  test  section. 

Spherical  and  cylindrical  Bhells  of  Pyrex  7740  glass,  ranging  in 
diameter  from  3  to  20  mm,  were  employed  as  models.  Two  thin  film  re¬ 
sistance  thermometers  were  applied  at  the  stagnation  point  of  the  cylinder 
and  half  an  inch  on  either  side  of  the  mid-point  of  the  cylinder.  The  films 
were  applied  by  airbrushing  Hanovin  No.  05X  liquid  bright  platinum  paint 
in  a  strip  0.5  mm  wide  and  half  an  inch  long  and  then  firing  it  1250°F. 
Firing  silver  micro  paint  in  a  heavy  layer  provides  low  resistance  contacts 
with  the  measuring  circuit.  Fig.  31  is  a  photograph  of  the  cylindrical 
model.  This  model  was  mounted  in  the  test  section  normal  to  the  flow  so 
that  the  mid-point  of  the  cylinder  corresponded  with  the  axis  of  the  tunnel 
test  section.  The  cylinder  was  supported  from  the  wall  in  cantilever  fash¬ 
ion.  The  spherical  model  shown  in  Fig.  32  has  a  shorter  thin  film  at  its 
stagnation  point  and  is  supported  from  the  model  stand  as  illustrated.  The 
pitot  pressure  gauge  location  is  indicated  also  in  this  photograph. 
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‘igure  32.  Mounting  of  the  Spherical  Model  in  the  Tunnel  Showing  also 
the  Pitot  Tube  for  Measuring  Test  Section  Total  Pressure 


The  technique  of  thin  film  resistance  thermometry  is  quite  well 
known  and  very  adequately  covered  in  the  literature.  In  essence,  the  ther¬ 
mometer  measures  the  surface  temperature  history  of  the  material  on 
which  it  resides.  This  measured  temperature  history  is  converted  to  a 
heat  transfer  rate  by  the  use  of  one- dimensional  heat  conduction  theory  and 
one  basic  calibration  of  the  thermometer.  The  basic  calibration  is  gener¬ 
ally  done  by  the  pulsing  technique  in  which  a  known  heat  flux,  produced  by 
a  current  step  pulse  through  the  resistance  thermometer,  is  applied  to  the 
gauge.  The  calibration  is  used  to  determine  the  lumped  thermal  properties 
of  the  substrate  of  the  thin  film,  “s/ir  p  Ck  ,  where  p  ,  C  ,  and  k  are  the  den¬ 
sity,  specific  heat,  and  conductivity  of  the  substrate.  The  thermal  coeffi¬ 
cient  of  resistance  of  the  film  is  determined  by  measuring  the  resistance 
at  the  freezing  and  boiling  points  of  water. 

The  equation  relating  the  boat  transfer  rate  (O)  to  the  model  along 
with  the  temperature  history  of  the  model  surface  (AT/'v/t)  is  given  as 


Q  = 


n/tt  p  C  k 


‘ 


1  + 


AT' 


o 

AT7 


log 


10 


1  + 


AT 

T 


.414 


A-  (127) 
vt 


AT  is  the  temperature  change  of  the  surface  and  T  is  the  initial  tempera¬ 
ture  of  the  surface.  The  quantity  in  the  braces  ^  ^  represents  the  correc¬ 

tion  to  the  heat  transfer  rate  which  must  be  included  in  order  to  account  for 
the  change  in  the  thermal  properties  of  the  substrate  when  the  surface 

temperature  becomes  large  in  comparison  with  the  calibration  tempera- 
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ture.  b/k^  is  a  function  of  the  material,  which  for  Fyrex  is  4.73. 

Measurements  were  made  in  the  shock  tunnel  of  the  heat  transfer 
rates  to  cylinders  3mm,  10mm,  and  20mm  in  diameter  and  to  spheres  10 
and  20mm  in  diameter.  Reservoir  stagnation  conditions  ranged  from  2800 
to  3000°K  with  pressures  of  28  to  2 9  atm  and  220  to  240  atm.  Mach  num¬ 
bers  in  the  test  section  varied  from  10  to  23,  but  most  of  the  data  was  in 
the  range  of  13  to  16. 
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4.  Results 


In  calculating  the  value  of  Re^ ,  the  Reynolds  number  based  on 

leading  edge  diameter,  the  value  of  the  viscosity  is  determined  according 

to  the  results  of  the  study  of  oxygen  recently  reported  by  Hartunian  and 
36 

Marrone.  This  study  indicates  that  the  Sutherland  viscosity  values  based 

on  kinetic  theory  and  hard  spherical  molecules  are  as  much  as  50%  too  low 

for  temperatures  up  to  3000°K.  This  modified  value  of  u  also  was  used 

1  2 

to  determine  the  theoretical  value  of  Qq  ,  according  to  Fay  and  Riddell. 

The  experimental  values  of  the  reservoir  conditions,  the  flow 
Mach  numbers  and  the  heat  transfer  rates  are  tabulated  in  Table  I.  Values 
of  QQRet/-p  U'  as  a  function  of  Re^  are  plotted  in  Fig.  30a  and  30b,  along 
with  the  theoretical  curves.  The  results  for  the  cylinder  seem  to  follow 
the  boundary  layer  theory  right  on  down  to  the  free  molecule  flow  regime 
with  only  a  slight  increase  at  Reynolds  numbers  less  than  20.  The  consid¬ 
erable  scatter  in  the  data  makes  a  definitive  determination  of  the  best 
theory  impossible. 

The  experimental  data  for  the  axiaymmetric  probe  (sphere)  are 
even  more  scattered,  and  therefore,  inconclusive. 

Some  remarks  on  the  accuracy  of  these  experiments  should  be 
made.  The  shock  Mach  number  can  be  measured  to  &  1%  and  this  is  about 
the  error  expected  in  the  reservoir  conditions.  The  total  pressure  can  be 
measured  to  within  ±  10%  and  this  should  be  the  expected  error  in  p  .  As 
a  result,  ,the  theoretical  curve  which  is  plotted  from  the  experimentally 
obtained  values  for  the  tunnel  conditions  could  have  a  width  of  ±  10%.  The 
heat  transfer  measurements  have  a  possible  error  of  between  ±  5%  and 
dt  18%,  depending  on  the  smoothness  and  clarity  of  the  heat  transfer  record. 
An  example  of  a  typically  adequate  record  is  shown  in  Fig.  33  for  a  flow 
Mach  number  of  16.5  and  reservoir  conditions  of  243  atm  and  3100°K. 

No  explanation  is  immediately  available  for  the  greater  amount  of 
scatter  in  the  sphere  data.  Since  the  spheres  were  blown  from  Pyrex  tub¬ 
ing  some  irregularities  existed  in  their  radii.  However,  the  heat  transfer 
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Figure  33.  Typical  Temperature  History  for  the  Stagnation  Point  of  a 

Cylindrical  Model 

P  =  222  atm 

T  =  2980°K 

M  =16.0 
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rate  varies  inversely  as  the  square  root  of  the  body  radius,  and  hence, 
i:i  regularities  in  the  radii  should  not  give  significant  effects.  The  much 
smaller  gauge  on  the  spheres  could  limit  the  accuracy  of  the  results. 
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5 ,  Conclusions 

The  heat  transfer  rates  to  cylinders  and  spheres,  as  measured  in 

the  shock  tunnel,  are  adequately  predicted  by  the  theory  of  Fay  and  Riddell 

down  to  Re  £  20  .  At  this  point,  the  heat  transfer  rate  increases  some- 
1  34 

what,  as  predicted  by  Probstein  and  Kemp,  reaching  the  free  molecule 
flow  heat  transfer  rate  at  a  value  of  Ret  between  2  and  3. 
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Where  I  is  the  beam  intensity,  x  is  the  path  length,  p  is  the  gas  density, 
and  p  is  an  absorption  coefficient. 

The  absorption  coefficient,  p,  is  related  to  the  effective  scatter¬ 
ing  cross  section  of  the  gas.  The  cross  section  Q  ($)  is  defined  as  the 
probability  per  scattering  center  that  an  electron  will  be  deflected  by  an 
interaction  with  a  gas  molecule  into  a  cone  of  half  angle  $  which  is  greater 
than  the  effective  angle  subtended  by  the  detector  aperture  measured  from 
the  point  of  scattering.  If  there  are  n  molecules  of  gas  in  the  volume 
occupied  by  the  electron  beam,  the  total  probability  of  an  electron  inter¬ 
action  with  a  gas  molecule  is  n  Q  ,  or  over  an  increment  of  path  length  dx, 
is  nQdx,  and  the  loss  of  beam  intensity  becomes: 


^  =  nQI  (129) 

dx 


Integration  of  Eq.  (129)  yields  I  (x) 


I 

o 


(-nQx 


(130) 
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Markevitch  and  Hurlbut  have  shown  that  the  experimental  value 
of  Q  (it)  for  high  energy  electrons  in  air  is  close  to  the  theoretical  value 
based  on  elastic  collisions.  This  value  is  energy  dependent  and  has  been 
found  to  have  the  relation 


Q  (it) 


3.6  x  IQ"13 
V 


2 

cm 


(131) 


where  V  is  in  volts. 

23 

Taking  Avogadro's  number,  N  =  6.025  x  10  /mol.,  and  the 

3  3 

standard  volume,  V  =  22.4146  x  10  cm  /mol.,  then  at  a  pressure  of  one 
°  13 

micron  there  are  3.537  x  10  molecules  per  cubic  centimeter  of  a  perfect 
gas.  Then  n  =  3.537  x  10^  p.  and  pp  =  nQ  =  3.537  x  10^  pQ,  from 
which 

Q  =  2.827  x  10'14  .  (132) 

P 
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This  quantity  Q  is  the  total  collision  cross  section  of  the  gas,  and  has  the 
same  significance  as  the  sense  of  Q  ( it )  in  Eq.  (131);  that  is,  it  represents 
the  total  cross  section  which  was  determined  from  scattering  into  a  cone 
or  half  angle  it  radians. 

We  may  now  obtain  the  absorption  coefficient  p  by  combining 
Eq.  (131)  and  (132); 


pp  =  12.73  -E 
^  V 


(133) 


for  air,  where  p  is  in  microns  and  V  is  in  volts. 

The  equation  for  the  intensity  of  an  electron  beam  in  air  at  a 
constant  pressure  p  microns  is  then 


I  (x) 


_  ~  1 2 . 7  3  px 

I  r 


(134) 


where  I  (x)  is  measured  in  centimeters  from  the  point  of  origin  and  Iq  is 
the  incident  beam  current  of  an  energy  V  electron-volts. 

Equation  (134)  is  plotted  on  logarithmic  coordinates  in  Figs.  34 

and  35. 

3.  Electron  Beam  Absorption  in  an  Axially  Symmetric  Flow  Field 

The  electron  beam  densitometer  will  be  used  to  obtain  free  stream 
density  measurements  in  the  flow  field  of  the  hypersonic  shock  tunnel  and 
the  variation  of  free  stream  density  in  the  direction  normal  to  the  flow. 
Consider  a  uniform  beam  of  electrons  injected  into  an  axially  symmetric 
flow  field  where  the  gas  density  is  a  function  of  r  ,  the  radial  distance 
from  the  central  flow  axis,  as  in  Fig.  36.  It  will  be  required  to  find  the 
gas  density,  p  ,  as  a  function  of  r  .  Consider  a  uniform  beam  of  electrons 
of  intensity  I  originating  outside  the  flow  region  and  directed  in  the  posi¬ 
tive  x  direction. 
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Lambert's  law  of  absorption  holds  as  a  partial  derivative  since 
the  beam  attenuation  is  a  function  of  x  and  y  . 


61^ 

6x 


-  M-Pl 


(135) 


which  has  a  solution 


lnl(y)  -  lnIQ  = 


pdx  -•  fj-p x  (L  -  2x)  =  F(y) 


(136) 


where  p  =  f(r)  and  p^  is  a  constant. 
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From  the  geometry  of  Fig.  36,  r  =  x  +  y  ,  and  for  constant  y  , 
rdr  =  xdx  =  n/t2  -  y2  dx-  Also,  when  x  =  o  ,  r  =  y,  and  when  x  =  x, 
r  =  R  then  Eq.  (136)  becomes 


F(y) 


p(r) ...... .rd£ —  -  |ip  (L  -  2x) 

Tl 2  l 

Vr  -  y 


(137) 


A  second  change  of  variables  ib  made  by  putting  s  =  R  -  r  and 

2  2 

g(s)  =  p(r).  Now  ds  =  -2rdr,  and  at  r  =  y ,  s  =  R  -y  and  at  r  =  R, 
s=0.  Then  Eq.  (137)  becomes 


F(y) 


K-V 


ds 


g(s)  — _ 

2  2  >p  r 


-  PP1  (L  -  2x) 


(138) 


R  -  y  -  s 


2  2 

This  may  be  simplified  by  substitution  of  t  =  R  -  y  ,  then 


t 

F(y)  =  -  p  f  -  (ip,  (L  -  2x) 

J  \Jt  ~  B 


(139) 


115 


and  define 


f(y)  =  G(t) 


H  f  &!*)*? 
->  *Jl  -  8 

o 


(HO) 


then 


f(y)  =  F(y)  +  ppj  (L  -  2x) 


(141) 


Equation  (140)  is  of  the  Abel  type  with  unknown  g { s )  and  has  the 

solution 


B(a) 


zlA  f 

TTfi  ds  J 


G(t)dt 
V  8  -  t 


To  integrate  by  parts,  choose  U  =  G  (t)  and 


dV 


Equation  (142)  becomes 


g(s) 


a 

li  °(°)  +  L  f  G'  Ujdt 

«/7  "ffH  J  n/TTT 

o 


From  Eq.  (140)  G  (o)  =  o  and 


G'  (t) 


so  that 


g  (s) 


Z 


S 

_L  f  f'(y)dt 
^  J  \l  8  -  t 

o 


(142) 
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(144) 
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We  now  change  variables  by  the  substitutions  t  =  R  -  y  and 
2 

-  r  .  Note  that  when  y  =  R  ,  t  -  o  and  when  y  =  r  ,  t  =  s  ,  then 


g(s)  =  p  (r )  =  -L  f  ---fy)dy  (145) 

^  J  fz  2 

r  Vy  -  r 

By  differentiation  of  Eq.  (141), 

f '  (y)  =  4-  lnl(y)  -  2 ixp  y  (146) 

dy  O  2 

V  R  -  y 

and  substitution  into  Eq.  145, 


p(r )  = 


1T(J. 


lnl(y)  — — 


(147) 


which  leads  to 


p(r) 


dy 


lnl(y) 


dy 


+  Pi 


,  2  2 
\(y  -  v 


(148) 


I  (y)  may  be  obtained  by  measuring  segments  of  the  core  flow 
under  identical  tunnel  conditions  and  p(r)  is  evaluated  numerically.  A 
preferred  technique  is  to  incorporate  a  series  of  detectors  and  a  rapid 
sweep  on  the  electron  beam  to  measure  all  segments  during  a  single  run. 

4  Instrumentation 

The  electron  beam  densitometer  consists  of  two  major  systems: 
the  electron  source  or  gun,  with  the  associated  high  voltage  accelerator, 
heaters,  bias  and  deflection  power  supplies,  and  the  detector,  with  its 
associated  power,  bias,  and  deflection  supplies  Oil  diffusion  pumps  sup¬ 
ply  the  high  vacuum  environment  for  the  electron  gun  and  the  detector 
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helium  in  the  driver,  and  a  reasonable  degree  of  driver  gas  heating.  The 
primary  driver  gas  will  then  be  helium  with  a  stoichiometric  mixture  of 
hydrogen  and  oxygen  to  supply  the  energy  for  heating  and  to  increase  the 
driver  gas  pressure  through  a  combustion  process.  In  the  operation  of  the 
combustion  driver,  however,  the  helium  is  regarded  as  the  diluent  in  sup¬ 
pressing  detonation  during  the  combustion  process.  Thus,  a  combustion 

mixture  is  usually  referred  to  by  its  helium  content  in  percent  of  the  total 

40 

volume.  Nagamatsu  has  performed  experiments  for  hydrogen-oxygen- 
helium  mixtures  of  various  ratios.  Detonation  was  experienced  when  a 
65%  helium  diluent  was  used  with  a  stoichiometric,  mixture  of  hydrogen  and 
oxygen.  Difficulties  in  igniting  the  gas  were  encountered  when  an  80% 
helium  dilution  was  used. 

In  calculating  the  performance  of  the  shock  tunnel,  helium  is  re¬ 
garded  as  the  driver  gas.  The  product  of  combustion,  vapor  of  H.,0,  is 
then  regarded  as  the  contaminant  in  the  driver  gas.  Since  the  water  vapor 
contributes  to  an  appreciable  increase  in  molecular  weight  of  the  final  gas 
mixture  over  that  of  the  pure  helium,  its  effects  cannot  be  ignored  com¬ 
pletely  in  the  calculations.  A  detailed  Btudy  of  the  combustion  process 
would  give  an  accurate  prediction  of  the  driver  gas  conditions  with  regard 
to  temperature,  pressure,  speed  of  sound,  specific  heat  ratio,  etc.  A 
detailed  study  of  this  nature,  in  order  to  be  adequate  for  our  present  pur¬ 
poses,  should  trace  the  expansion  process  of  the  driver  gas  after  the  shock 
formation  has  been  initiated.  This  implies  the  necessity  for  the  construc¬ 
tion  of  portions  of  a  Mollier  diagram  for  each  mixture  ratio  under  consid¬ 
eration,  in  order  that  the  regions  of  interest  can  be  adequately  covered. 

For  our  present  purpose,  the  equilibrium  gas  temperature,  pressure,  and 
specific  heat  ratio  y  were  determined  at  the  time  when  the  combustion 
process  was  completed.  Isentropic.  expansion  of  the  driver  gas  was  assumed 
in  order  to  calculate  the  conditions  behind  the  interface.  Four  mixture 
ratios  were  considered  in  calculating  the  driver  gas  conditions,  and  subse¬ 
quently,  the  shock  tunnel  performance  with  regard  to  stagnation  conditions 
in  the  region  after  the  reflected  primary  shock.  It  is  on  this  stagnation 
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Calculations  for  a  shock  tunnel,  using  pure  hydrogen  as  the  driver 
gas,  is  given  in  Fig.  49  as  a  comparison.  Hydrogen  is  very  efficient  for 
use  as  a  driver  gas;  however,  its  inherent  dangerous  handling  character¬ 
istics  make  it  unattractive  for  use  in  the  present  facility. 

a.  Area  Ratio  Effect  in  Shock  Tunnels 

The  preceding  analysis  was  appropriate  for  an  ideal  shock 
tube  with  constant  area  throughout.  For  shock  tunnels  having  unequal 
driver  and  driven  tube  diameters  and  a  throat  section  at  the  diaphragm 
location,  the  performance  calculated  previously  would  be  modified  to  com¬ 
pensate  for  area  ratio  effect  and  the  restriction  at  the  throat.  A  shock  tube 

with  such  a  restriction  is  given  schematically  (see  Fig.  50).  The  effective 

42 

pressure  ratio  for  such  a  configuration  would  be 


lpl)  p3a  p3*  p3b  p3  p2  P1 


1  p2 

V4  *  1  &2 

1  -  Ms 

[l  -  — ]  g"Q 

"8  pl 

L  a4 

l  pJ8  J 

where 


'2  +  (v4-  dmL' 

2 

-2  +  (y4-  l)M3b- 

L2  +  (V4-l)M23b_ 

_2  +  (v4  -  DM3a_ 

(153) 


(154) 


is  the  gain  factor  modifying  the  actual  pressure  ratio  to  an  effective  pres¬ 
sure  ratio.  With  the  boundary  conditions  such  that  =  1  and  if 

A3  is  given,  then  M.<a  can  be  obtained  by  solving  the  equation 


A, 

3a 

2+(Y4-l)M*a" 

2(y4-1)  m3* 

A  3# 

_2  +  (y4  -  1)M^_ 

M3a 

(155) 
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SHOCK  MACH  NUMBER,  M, 


Figure  54.  Tailored  Interface  Operation  for  the  Combustion- 
Driven  Shock  Tunnel  (P^~' 
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system,  the  driven  section,  and  the  nozzle  test  section  dump  tank  is  given 
below: 

b .  Driver  Section  and  Ga s  Charging  System 

Due  to  the  possibility  of  detonation  of  the  combustion  mixture, 
the  driver  section  is  designed  for  an  ultimate  pressure  of  165,  000  psi.  The 
working  pressure,  in  general,  will  not  exceed  30,000  psi.  The  temperature 
is  expected  to  reach  between  2000°  and  3000°K,  but  it  is  transient,  as  is 
the  high  combustion  pressure. 

In  order  to  meet  these  requirements,  the  driver  section  and 
its  closely  associated  fittings  will  be  fabricated  from  forged  300M  steel  and 
heat  treated  after  rough  machining  to  maximize  the  effectiveness  of  the 
heat  treating.  The  vessel,  cap,  latching  mechanisr.i,  and  diaphragm  pack¬ 
age  are  detailed  in  Fig.  59. 

The  gas  handling  system  for  loading  the  driver  with  a  com¬ 
bustible  mixture  of  oxygen  and  hydrogen,  diluted  with  helium,  is  shown 
schematically  in  Fig.  60.  A  gas  accumulator,  with  a  volume  equal  to  that 
of  the  driver  section,  acts  as  a  fast  high-pressure  hydrogen  pump.  A  stoi¬ 
chiometric  mixture  of  hydrogen  and  oxygen  diluted  with  75  percent  helium, 
at  a  total  pressure  of  4000  psi  before  combustion,  is  used.  Oxygen  iB 
pumped  into  the  driver  section  from  the  storage  marked  to  a  pressure 
of  330  pal.  Helium  is  then  added  until  a  pressure  of  3330  psi  is  reached. 
Then  hydrogen  i.s  pumped  into  the  hydrogen  accumulator  to  a  pressure  of 
660  psi.  Helium  gas  then  passes  into  the  accumulator  behind  the  piston. 

As  the  helium  is  stor  ed  at  very  high  pressures  (15.000  -  20,000  psi),  it 
forces  the  hydrogen  into  the  driver.  This  sequence  of  loading  never  per¬ 
mits  a  mixture  of  hydrogen  and  oxygen  to  exist  in  richer  concentrations 
than  the  mixture  to  be  employed  in  the  test.  Hence,  the  probability  of  deto¬ 
nation  oc curing  during  the  filling  process  is  low. 

Proper  distribution  and  circulation  of  the  driver  charge  is 
obtained  by  using  a  long  ported  stainless  steel  tube  inside  the  driver  to 
serve  as  the  fill  tube.  The  port  area  varies  exponentially  along  the  tube, 
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for  each  of  the  two  elements 


Cy 


AY 

AFy 

2 


AX 

r  -  tan  0  _  Cx 

y  ”  2AFx  tan  0  “  Z  ^ 
- - - -  tan  <f> 


It  can  be  shown  that  the  angled  drive  system  is  equivalent  to  a 
mechanical  transformer  with  a  "turns  ratio"  of  tan  0/1.  This  satisfies  the 
conditions  derived  for  force  displacement  and  impedance.  The  equivalent 
circuit  for  this  configuration  is  given  in  Fig.  76.  It  is  seen  that  for  0  =  45°, 
the  turns  ratio  becomes  unity  and  the  performance  becomes  identical  to  that 
of  Fig.  7l  except  for  the  fact  that  the  single  bender  element  has  been  re¬ 
placed  by  two  with  the  same  aggregate  dimensions  and  mechanical  charac¬ 
teristics.  This  configuration  can  always  be  employed  with  a  resultant 
doubling  of  the  output  voltage,  but  since  the  internal  source  capacitance  is 
reduced  to  1/4  of  the  single  element  value  it  is  not  often  used. 

As  with  the  conventional  structure,  we  can  derive  the  low  fre¬ 
quency  sensitivity  This  is  found  to  be 
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The  factor  "2"  appears  in  the  denominator  because  C^.  refers  to 
each  element  of  the  two  element  system  It  can  be  seen  that  for  a  given  set 
,,f  diaphragm  and  crystal  constants  there  is  an  optimum  angle  which  gives 
maximum  output  This  relationship  is 
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